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CHAPTER  I 


INTRODUCTION 

The  concept  of  the  electrothermal  gun  was  first  developed  in  the  hopes  of  increasing 
the  muzzle  velocity  of  a  projectile  for  use  in  a  wide  variety  of  tasks.  Most  notably,  the 
use  of  this  concept  has  applications  in  the  area  of  both  space-based  hypervelocity 
"kinetic-ldll"  type  weaponry  as  well  as  more  conventional  tank  and  artillery  uses.  An 
electrothermal  gun  is  designed  to  utilize  the  interaction  of  an  electrically-generated 
plasma  or  working  fluid  to  create  the  high  temperatures  and  high  pressures  necessary  for 
projectile  motion.  "The  deposition  of  eleaiical  energy  into  the  gun  represents  a  way  of 
increasing  the  work  potential  of  the  driver  gases  as  compared  to  conventional 
propellants  " '  By  using  electricity  as  the  source  of  energy,  instead  of  the  conventional 
chemical  reactions  in  an  ordinary  propellant,  the  working  fluid  may  be  chosen  for  its  low 
molecular  weight,  high  sound  speed,  favorable  base  pressure-to-breech  pressure  ratios, 
and  other  beneficial  thermal  characteristics  without  regard  to  the  chemical  reactions 
required  to  form  the  working  fluid.  The  "tailoring"  of  the  amount  of  energy  and  type  of 
working  fluid  represents  a  significant  advancement  over  that  of  conventional  propellant 
commonly  used  in  projectile  devices. 

Typically,  electrothermal  guns  offer  a  substantial  increase  in  muzzle  velocity  over  that 
attained  through  the  use  of  conventional  propellants.  As  an  example,  most  conventional 
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t>7e  propellants  can  impart  velocities  to  projectiles  on  the  ordei’  of  2  to  3  km/sec, 
v'hereas  electrothermal  guns  have  been  demonstrateo  to  attain  projectile  muzzle  velocities 
in  excess  of  6  to  7  km/sec.  In  fact,  under  ideal  circumstances,  it  is  theoretically  possible 
to  extend  these  velocities  for  small  mass  projectiles  to  over  9  km/sec. " 

In  the  attempt  to  harness  these  advantages,  several  classes  of  electrothermal  gun  have 
been  developed  over  the  last  few  decades.  These  include  Combustion-Augmented 
Plasma,  or  "CAP",  guns,  Thermal-Electric  (ET)  guns,  and  the  Electrically  Powered  Light 
Gas  Gun  (ELGG).  Each  of  these  devices  use  slightly  different  working  characteristics, 
but  are,  in  reality,  variations  on  the  same  theme. 

This  thesis  examines  a  new  variation  on  the  ELGG  concept,  that  differs  in  both  the 
initial  configuration  of  the  device,  as  well  as  the  development  of  the  working  fluid  or 
plasma  within  the  "propellant  chamber".  The  "propellant  chamber"  is  synonymous  with 
the  breech  of  the  gun.  The  "propellant  chamber"  in  the  case  of  the  ELGG  has 
temperatures  of  up  to  several  electron  volts  and  pressures  ranging  from  850  to  1, 150 
atmospheres. 

The  extremely  high  velocities  and  subsequent  high  operating  temperatures  and 
pressures  attained  by  electrothermal  guns  suggest  that  a  great  deal  of  heat  transfer  from 
the  plasma  to  the  gun  barrel  must  take  place  and  that  this  heat  transfer  is  likely  to  be 
much  greater  than  that  in  a  conventionally  powered  device.  However,  much  of  the 
previous  research  into  electrothermal  guns  has  not  adequately  addressed  this  concern,  and 
in  fact  many  of  the  published  documents  simply  assume  that  the  heat  transfer  is  of  the 
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same  order  of  magnitude  as  that  in  conventional  guns  with  radiation  assumed  to  be 
negligible.  Little,  if  any,  further  study  on  this  subject  has  been  conducted. 

It  is  the  aim  of  this  thesis  to  examine  the  new  class  of  electrothermal  gun  and  to  study 
its  operation.  Specifically,  the  heat  transfer  mechanisms  in  the  electrothermal  gun  barrel 
environment  are  examined  in  detail.  MACH2,  a  2  1/2  dimensional 
magnetohydrodynamic  simulation  code  developed  at  the  Phillips  Laboratories  (Kirtland 
AFB,  New  Mexico)  mnning  on  an  Ohio  State  Cray  supercomputer  will  be  utilized  in 
order  to  provide  highly  accurate  and  insightfiil  computer  simulated  data  for  various  heat 
transfer  modes  and  viscous  effects  within  the  barrel  section  of  the  ELGG.  The 
combination  of  first  principle  model  development  and  MACH2  data  will  yield  a 
substantial  understanding  of  the  ELGG  operational  and  heat  transfer  characteristics. 


1.1  ~  PROBLEM  STATEMENT 


The  basic  configiaration  and  operation  of  the  EL  GG  will  be  discussed  with  the  future 
application  of  a  2,100  capacitor  bank  located  at  the  Ohio  State  University ,  known  as 
"Godzilla",  to  the  ELGG  in  mind.  Additionally,  this  thesis  will  study  the  three  main  heat 
transfer  mechanisms;  thermal  radiation,  forced  convection  and  thermal  conduction  from 
the  high  temperature  plasma  "working  fluid"  to  the  walls  of  the  electrothermal  gun  barrel. 
Since  various  theoretical  heat  transfer  radiation  models  are  available,  the  determination  of 
which  model  is  appropriate  to  the  conditions  within  the  ELGG  barrel  will  be  made,  and 
the  comparison  between  the  various  heat  transfer  modes  will  be  derived  and  their  regions 


4 


of  influence  described.  Viscous  effects  will  be  taken  into  account,  and  the  combined 
energy  loss  due  to  the  main  transfer  mechanisms  and  the  viscous  effecis  will  be  discussed. 

The  impact  of  the  heat  transfer  mechanisms  and  viscous  losses  will  be  discussed  and 
demonstrated  from  predicted  data  obtained  with  models  developed  in  this  thesis. 
Consideration  will  be  given  to  the  degree  of  ionization  and  dissociation  created  by  the 
changing  temperature  and  pressure  profiles  throughout  the  length  of  the  barrel.  Finally, 
this  thesis  vrill  discuss  the  applicability  and  useability  of  the  magnetohydrodynamic 
simulation  code,  MACH2,  to  the  problem  of  high  heat  transfer  and  viscous  effects  within 
the  barrel  region  of  an  ELGG. 

1.2  -  PRELIMINARY  DISCUSSION 


The  system  described  in  this  thesis  uses  hydrogen  gas  as  the  working  fluid,  unless 
specified  otherwise,  and  avoids  the  use  of  chemical  reactions  as  the  source  of  energy  by 
introducing  an  arc  discharge  in  the  breech  chamber.  The  region  prior  to  the  beginning  of 
the  ELGG  barrel  is  not  studied  in  this  thesis,  but  it  is  assumed  that  hydrogen  gas  at  a 
specific  stagnation  pressure  and  temperature,  as  well  as  inlet  velocity,  is  provided  as 
initial  conditions.  The  analysis  deals  with  a  hydrogen  plasma  at  initial  stagnation 
temperatures  ranging  from  one  to  three  electron  volts  (eV)  and  at  stagnation  pressures 
from  850  to  1,150  atmospheres.  These  initial  parameters,  in  addition  to  a  set  of  assumed 
velocity  profiles  at  the  i.iiet  are  consistent  with  the  flow  exiting  the  chamber.  The  initial 
conditions  pose  extreme  heat  transfer  problems  within  the  barrel  region  of  the  gun,  as 
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well  as  the  necessity  to  contain  the  force  on  the  barrel  walls  caused  by  the  high  pressures. 
The  actual  heating  of  the  working  fluid  by  an  electric  arc  is  not  discussed. 

Hydrogen  gas  was  chosen  as  the  working  fluid  for  several  important  reasons.  The  low 
molecular  weight  of  the  gas  is  one  of  the  keys  to  attaining  very  high  projectile  velocities, 
and  the  corresponding  high  sound  speed  of  hydrogen  aids  in  the  interaction  of  the  hot 
plasma  gas  and  the  projectile  by  decreasing  the  wave  transit  time  thioughout  the  working 
fluid,  further  increasing  the  possible  projectile  velocity.  Because  of  the  choice  to  use 
hydrogen  ';,2s,  the  thermal  characteristics,  opacities,  and  thermodynamic  properties  of 
hydrogen  need  to  be  well  understood. 

In  addition  to  gaining  a  good  understanding  of  the  properties  of  hydrogen  at  the 
conditions  necessary  for  electrothermal  gun  operation,  the  main  modes  of  heat  transfer 
need  to  be  well  understood.  The  nature  of  these  modes  is  vital  to  the  understanding  of 
the  overall  heat  transfer,  and  governs  the  exit  kinetic  energy  the  projectile  may  acquire 
from  interaction  with  the  plasma. 

The  importance  of  studying  the  heat  transfer  mechanisms  is  dramatically  illustrated  by 
the  results  of  an  experiment  by  Tidman  and  Massey.  ’  They  developed  a  similar  ELGG 
with  a  barrel  length  of  5.0  m,  bore  diameter  of  1 .6  cm,  and  chamber  volume  of 244  cm^ 
in  which  a  projectile  with  a  mass  of  10  grams  was  accelerated  to  speeds  of  about  7 
km/sec  with  an  efficiency  (defined  as  the  transfer  of  the  energy  contained  within  the 
propellant  gas  to  kinetic  energy  of  the  projectile)  of  32%.  This  indicates  that  fully  68%  of 
the  energy  initially  stored  in  the  plasma  is  lost  due  to  a  variety  of  sources.  Most  of  these 
are  in  the  form  of  heat  convection  out  of  the  bn.',  el  walls,  heat  transfer  to  the  projectile. 
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conversion  of  energy  into  internal  energy  modes,  energy  of  ionization,  and  dissociation 
energy.  This  represents  a  considerable  amount  of  energy  which  must  be  accounted  for.  ^ 

The  main  problem  with  describing  the  heat  transfer  mechanisms  stems  from  the  large 
pressure,  hence  density  gradients  throughout  the  length  of  the  barrel,  as  well  as  the  large 
temperature  gradients  between  the  plasma  and  the  surrounding  barrel  walls,  especially 
near  the  wall  surface  within  a  thermal  boundary  layer.  In  addition,  the  effects  of  each  of 
the  heat  transfer  modes  is  difficult,  if  not  impossible  to  uncouple  from  each  other.  As  the 
plasma  travels  down  the  barrel,  the  movement  of  the  projectile  changes  the  volume  the 
plasma  must  occupy,  thus  further  decreasing  the  pressure  and  temperature  in  the  plasma 
itself,  resulting  in  larger  gradients  in  the  direction  of  the  flow  and  less  plasma  energy 
available  to  be  imparted  to  the  projectile. 

Electrothermal  guns  operate  in  regions  of  temperature  and  pressure  where  more 
conventional  artillery  guns  do  not  function  As  a  consequence  of  this,  very  little  data  for 
thermodynamic  and  opacity  properties  of  hydrogen  has  been  collected  v»ith  guns  in  mind. 
However,  much  of  the  data  used  in  this  thesis  was  obtained  from  studies  of  stellar 
surfaces,  which  in  many  ways,  has  similar  environmental  factors  to  those  in  an 
electrothermal  gun  barrel. 

In  addition  to  losses  due  to  heat  transfer,  viscous  effects  are  also  important  and  are 
studied.  Problems  associated  with  viscous  effects  include  loss  of  flow  enerf^  and 
possible  choking  of  the  flow  within  the  barrel  region  of  the  ELGG.  Regimes  of  operation 
in  which  choking  becomes  a  problem  will  be  discussed  and  will  be  correlated  with  heat 


transfer  studies  as  well. 


CHAPTER  II 


BACKGROUND 


In  order  to  have  a  good  understanding  of  what  the  system  physically  looks  like,  and 
what  regions  of  the  gun  the  considerations  of  heat  transport  are  important,  it  is  necessary 
to  understand  the  operational  characteristics  of  an  electrothermal  gun.  There  have  been 
several  concepts  of  electrothermal  guns,  and  what  differentiates  one  design  from  another 
is  also  important  to  understand,  l  o  begin  v.ath,  an  electrothermal  gun  is  an  advanced 
concept  which  involves  addition  of  electrical  energy  via  a  plasma  working  fluid  to  a 
projectile.  A  working  fluid  is  merely  a  gas  or  plasma  which  imparts  energy,  in  this  case, 
to  a  projectile. 

The  working  fluid  maintains  the  high  temperatures  and  pressures  necessary  for  the 
acceleration  of  a  projectile  dovvn  a  gun  barrel.  The  separation  of  the  tasks  of  power 
generation  and  gas  generation  (both  accomplished  by  one  propellant  source  in  a 
conventional  gun  approach)  is  inherently  more  flexible,  allowing  the  tailoring  of  the 
working  fluid  to  products  with  low  molecular  weight.  In  the  case  of  the  ELGG  studied  in 
this  thesis,  hydrogen  is  the  propellant  selected  due  to  it's  very  low  molecular  weight  of 
2.016x10^  Kg/Kmol.  The  use  of  a  low  molecular  weight  propellant  helps  to  minimize  the 
pressure  losses  along  the  gun  barrel. 
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A  major  advantage  of  the  ELGG  is  that  by  separating  the  sources  of  thermal  energy 
and  working  fluid,  it  is  possible  to  select  a  working  fluid  based  on  its  physical  and 
thermodynamic  characteristics  without  worry  as  to  the  energy-producing  chemical 
reactions  taking  place.  In  addition,  if  the  assumption  that  the  rate  of  production  of  the 
working  fluid  is  "directly  related  to  the  rate  of  deposition  of  electrical  energy  into  the 
system,  pressure  time  traces  can  be  much  more  readily  optimized  in  terms  of  producing 
maximum  muzzle  velocity  for  a  fixed  maximum  allowable  breech  pressure."  “ 

2.1  --  mSTORY  AND  CONnGURATIONS  OF  ELECTROTHERMAL  GUNS 

"  1  hermal-electric,  or  ET  propulsion  as  it  is  known  today,  has  its  origins  in  the  1 950's 
and  60's  when  researchers  attempted  various  means  of  utilizing  electrical  energy  in  the 
interior  of  a  gun  tube.  Primarily,  these  experiments  involved  the  conversion  of  electrical 
energy  to  thermal  energy  of  the  propelling  gases  in  the  breech  or  barrel."  During  the 
1970's,  electromagnetic  propulsion  (EM)  devices  were  studied  and  developed  .  Several 
key  technologies  were  developed  which  were  then  adapted  for  use  in  ET  devices,  most 
notably  the  advanced  power  supplies  required  by  both  types  of  systems. 

The  1980's  saw  increased  interest  in  ET  devices,  and  the  development  of  the 
Combustion-Augmented  Plasma  (CAP)  gun  provided  renewed  study  of  this  technology. 
The  CAP  gun  is  "typically  a  hybrid  system  combining  the  attributes  of  electrothermal  and 
liquid  propellant  systems.  Electrical  energy  stored  in  a  capacitor  bank  is  discharged 
through  a  breech  mounted  cartridge.  A  capillary  section  within  the  cartridge  consists  of 


9 


the  anode/cathode  terminals,  consumable  liner,  and  a  solid  fuel  source.  The  electrical 
discharge  strikes  a  plasma  arc  within  the  capillary  which  ablates  the  liner  material  and 
vaporizes  the  fuel  source.  A  plasma  jet  is  formed  and  accelerated  through  a  nozzle  and 
into  the  combustion  chamber.  Through  the  interaction  of  the  rich  plasma  jet  and  liquid 
liner,  thermal  energy  is  transferred."  ^  A  schematic  of  a  typical  CAP  gun  demonstrates  its 
basic  configuration; 


Figure  2.  i  CAP  Gun  basic  configuration 


The  above  description  demonstrates  that  this  type  of  electrothermal  gun,  although 
using  a  plasma  as  part  of  the  gas  used  for  projectile  motion,  is  really  a  different  type  of 
chemical  propellant  device  and  is  not  too  far  removed  from  conventional  artillery  guns. 
Indeed,  it  has  been  discussed  that  the  main  fraction  of  energy  within  the  accelerating  gas 
is  formed  from  chemical  reactions,  and  the  plasma  has  little  if  any  effect  on  the  energy 
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content  of  the  gas.  Juhsz  notes  that,  "one  strong  specific  concern  [of  the  electrothermal 
community]  was  that  with  a  large  percentage  of  the  propulsion  energy  coming  from  the 
chemical  source  (a  common  figure  was  a  75/25  chemical  to  electrical  energy  ratio)  the 
electrothermal  [CAP]  gun  could  reduce  to  a  modem  version  of  the  ill-fated  bulk  loaded 
liquid  propellant  gun  concept."  ^ 

The  system  described  in  detail  in  this  thesis  uses  hydrogen  gas  as  the  sole  working 
fluid,  unless  noted  otherwise,  and  does  not  rely  on  any  chemical  reactions  for  the 
de/elopmciit  of  the  plasma.  Typically  the  hydrogen  is  at  stagiiation  pressures  of  850  to 
1,150  atmospheres  and  stagnation  temperatures  initially  betv/een  one  eV  and  three  eV. 
This  system  is  generally  configured  as  a  chamber  separated  from  a  tube  (the  muzzle)  by  a 
diaphragm  designed  to  break  at  a  given  pressure. 

Hydrogen  gas  is  introduced  into  the  "combustion  chan'.ber"  (or  capillary,  or  breech), 
at  some  initial  pressure  and  temperature,  and  then  hydrogen  pla.sma  is  created  by  the 
interaction  of  the  initial  hydrogen  "charge"  with  an  electric  arc  created  between  an 
annular  anode  casing,  or  nozzle,  and  a  cathode  probe  A  wire  connects  the  anode  and 


cathode  to  facilitate  tiie  stait  of  llic  arc  discharge.  This  wire  explodes  with  the 
application  of  electrical  current  (usually  within  the  first  few  microseconds),  but  its  short 
life  span  is  sufficierit  to  create  the  path  for  the  arc  to  develop  It  is  assumed  in  this  thesis 
that  the  mass  of  the  exploding  wire  atoms  does  not  significantly  impact  the  propellant 
mass,  and  therefore  does  not  alter  the  molecular  weight  of  the  working  fluid,  nor  does  it 
impact  the  thermal  characteristics  and  the  ratio  of  specific  lieats  of  the  hydrogen  gas 
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Once  initiated,  the  arc  is  self  sustaining  as  long  as  electrical  energy  is  provided.  "The 
actual  thermal  conversion  takes  place  when  electrons  which  are  driven  from  one  end  of 
the  capillary  [or  breech]  to  the  other  collide  with  ions  or  atoms  in  the  plasma.  Titis 
increases  the  energ>'  of  the  plasma  constituents  which  is  directly  related  to  the  plasma 
temperature  The  collision  rate  and  the  number  of  electrons  available  define  the  electrical 
resistance  of  the  disciiarge  so  that  the  usual  relation  for  ohmic  heating  holds  " '  It  is 
assumed  that  the  electrical  energy  is  uniformly  distributed  throughout  the  working  fluid. 

A-;  a  specified  pressure  (usually  between  850  and  1, 150  atmospheres),  the  diaphragm 
breaks,  allowing  the  energetic  plasma  to  expand  either  directly  into  the  ELGG  barrel  or 
through  a  conveigent-diverge.ii  iiozz’e  and  into  the  barrel  where  the  fluid  interacts  with 
the  ptojectile.  This  inrerai^tion  takes  place  as  long  as  the  projectile  remains  in  the  muzzle, 
and  the  effects  of  muzcle  wall  interactions,  as  well  as  the  increasing  volume  due  to 
projectile  motion  and  heat  transfer,  must  be  taken  into  account.  The  main  components  of 
the  ELGG  art'  the  electric  power  source  (capacitor  bufik),  an  electrical  switch,  the 
"combustion  or  arc  chamber"  with  anode  casing  and  cathode  probe,  diaphragm,  muzzle 


or  barrel,  and  tho  projectile.  The  followirig  figure  graplucalliy  lllusiiaies  a  schematic  a 


typical  ELGG.  Note  that  tiiis  figijre  does  not  illustrate  a  convergent-divergent  nozzel 
section  prior  to  the  ban ci  inlet. 


12 


Figure  2.2  Electrothermal  Light  Gas  Gun  configuration 


The  kinetic  energy  transferred  to  the  projectile  is  a  function  of  the  velocity  attained 
throughout  the  muzzle.  "It  is  generally  accepted  that  the  practical  velocity  limit  of 
chemical  propulsion  is  on  the  order  of  two  kilometers  per  second.  It  is  also  thought  [by 
the  requirements  for  higher  velocity  weaponry]  that  this  may  not  be  adequate.  Three 
faaors,  propeiiant  energy  density,  progressivity,  and  the  partitioning  of  kinetic  energy 
between  the  propellant  gases  and  projectile  are  seen  as  potential  contributors  to  this 
velocity  limit." ' 

"I  he  first  parameter,  propellant  energy,  is  determined  by  the  heats  of  formation  and 
the  relative  amounts  of  propellant  constituents.  Given  the  mature  state  of  the  art  in 
energetic  materials,  extending  the  energy  limits  of  new  formulations  significantly  may  not 
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energetic  materials,  extending  the  energy  limits  of  new  formulations  significantly  may  not 
be  practical.  While  formulations  with  higher  energies  are  certainly  possible,  their 
potentially  increased  hazard  properties  may  not  be  acceptable." ' 

The  second  parameter,  progressivity,  "i.e.  the  increase  in  gas  generation  rate  as 
burning  progresses  concerns  the  ceiling  on  the  amount  of  propellant  which  may  be 
profitably  burned  in  the  gun." '  Since  progressivity  is  governed  by  the  chemistry  or  the 
granulation  of  the  propellant,  there  is  an  optimum  limit  to  which  a  chemical  propellant  can 
impart  kinetic  energy  to  the  projectile.  The  energy  limitations  inherent  in  the  propellant 
chemistry  are  thereby  avoided  by  the  use  of  the  electrothermal  gun. 

"The  third  parameter,  partitioning  of  kinetic  energy  between  the  propellant  gases  and 
the  projectile  becomes  impcrtant  due  to  the  high  propellant  mass  (charge)  to  projectile 
mass  ratios  (C/M)  required  to  attain  high  velocities.  In  interior  ballistics  it  is  generally 
estimated  that  for  C/M  approxin  ately  1,  the  kinetic  energy  of  the  propellant  gases  is 
about  one  third  the  projectile  kinetic  energy.  Above  a  C/M  of  three,  the  kinetic  energy  of 
the  propellant  gases  begins  to  exceed  the  kinetic  energy  of  the  projectile,  making  high 
velociC'  systems  increasingly  less  efficient." ' 

The  following  figure  graphically  illustrates  the  effects  of  changing  any  combination  of 
the  propellant  molecular  mass,  C,  the  sound  speed  of  the  working  fluid.  A,  the  projectile 
mass  ,M ,  amd  the  projectile  muzzle  velocity.  V,  has  on  the  gas  expansion  ratio,  E.  ^ 
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Figure  2.3  Correlation  of  main  parameters ' 

"The  degree  of  expansion,  E,  is  a  paramster.  Sound  speed,  A,  in  electric-discharge 
gun  appears  to  be  much  more  than  powder  [conventional  propellant]  gun  mainly  because 
of  mass  of  light  gases  is  considerably  less,  and  to  be  considerably  higher  than  in  piston 
light-gas  gun  due  to  rather  high  temperature.  Therefore,  at  other  equal  (sic.)  conditions, 
the  higher  muzsde  velocity  is  obtained."  ^  The  following  table  illustrates  the  advantages 
of  the  ELGG  over  conventional  powder  (i.e.  chemical  reaction  dependent)  guns  and 
piston  light  gas  guns. 
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Table  2  .1  Comparison  of  three  different  types  of  guns.  (Note 
conventional  propellant  has  a  relative  number  of  1  for  comparison.) ' 


Parameter 

Conventional 
powder  gun 

Piston  light 
gas  gun 

ELGG 

Minimum 
molecular 
jmass  of  gas 

1 

0.06 

0.06 

Maximum 
temperature  of 
gas 

1 

CO 

O 

1 

O 

1.5 -2.0 

Chamber 
speed  of 
sound, A 

1 

2 

4-5 

Relative  mass 
of  projectile, 
M/C 

1 

2-3 

5-6 

Maximum 

projectile 

muzzle 

velocity  at  the 
same 
expansion 
ratio,  E 

1 

1.5-2 

2-3 

Weight  of 
barrel 

1 

2 

1 

n..  ■  ■■■■  j=as=:=;  > 

The  electrothermal  gun  concept  is  a  way  of  avoiding  the  disadvantages  of  conventional 
propellant  and  offer  the  ability  to  achieve  higher  velocity  projectiles.  Since  electrothermal 
guns  use  gas  pressure  to  impart  motion  to  the  projectile,  it  has  been  speculated  that 
conventional  barrel  technologies  and  materials  could  be  used.  However,  this  is  a 
discussion  which  needs  to  be  analyzed  by  considering  the  heat  transfer  from  the  hot 
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plasma,  which  will  be  accomplished  in  chapters  three  and  four.  Typical  electrothermal 
plasma  temperatures  can  range  from  lO.OOOK  (about  1  eV)  to  34,000K  (about  3  eV), 
whereas  typical  conventional  gun  flame  temperatures  are  in  the  4,000K  tc  5,000K  regime. 

"For  a  conventional  solid  propellant  gun  the  constraining  factors  on  performance  in 
terms  of  muzzle  velocity  are  gun  strength,  chamber  volume  and  tube  length  as  well  as  the 
combustion  characteristics  of  the  propellant.  The  gun  chamber  volume  limits  the 
pressures  at  which  the  gun  can  operate.  Burning  characteristics  of  the  propellant 
determine  the  pressure  distribution  throughout  the  gun.  An  electrothermal  gun  is  also 
limited  by  similar  constraints.  Gun  strength  again  limits  operational  pressures.  The  rate 
and  magnitude  of  the  electrical  input  in  combination  with  properties  of  the  working  fluid 
determine  the  pressure  profiles  in  the  gun.  The  [ELGGj  differs  from  the  solid  propellant 
gun,  however,  since  there  is  no  limitation,  at  least  in  theory,  on  the  total  available  energy 
for  the  system."  * 

"Although  the  amount  of  electrical  energy  input  is  theoretically  unlimited,  from  the 
practical  point  of  view,  the  maximum  operating  gun  temperature  limits  the  amount  of 
which  can  be  introduced  into  the  system.  In  fact,  thermal  management  for  an  [ELGG] 
gun  is  much  more  difficult  than  for  conventional  systems.  In  the  conventional  gun,  the 
overall  temperature  is  essentially  limited  by  propellant  flame  temperature.  Even  if  the  rate 
of  energy  input  increases  by  increasing  the  burning  rate  of  the  propellant,  the  gun 
temperature  is  still  limited  by  the  flame  temperature  which  is  independent  of  the  burning 
rate.  For  the  [ELGG]  gun,  on  the  other  hand,  there  is  no  such  upper  limit.  The 
temperature  of  the  gases  resulting  from  the  ...  electrically  heated  plasma  ...  is  an 
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increasing  function  of  the  amount  of  electrical  energy  being  transmitted  to  a  unit  mass  of 
working  fluid.  Thus,  temperature  limitation  results  only  through  controlling  the 
magnitude  and  rate  of  electrical  energy  input  in  combination  with  the  mass  and  properties 
of  the  working  fluid."  *  From  this  discussion,  it  is  apparent  that  the  magnitude  and  the 
rate  of  electrical  energy  input,  and  the  properties  of  the  working  fluid  are  the  controlling 
factors  in  the  gun  pressure  profiles  which  detemhne  projectile  velocity  and  overall  gun 
temperature. 


2.2  -  HYDROGEN  PLASMA  AS  A  WORKING  FLUID 


^l_  If  -1  *  «*/^*f* 

me  qucSu«.ui  Oi  Vviiy  *;iiuusc  iiyoiugeii  g<is  as  me  mam  wofKing  liuiQ  IS  easuy 
answered  by  considering  what  the  ideal  working  fluid  would  be.  "It  is  particularly 
desirable  that  gun  propellants  have  low  molecular  weight  since  such  products  result  in 
minimization  of  the  amount  of  product  mass  which  has  to  be  accelerated  through  the  gun 
barrel  along  with  the  projectile.  Accordingly,  low  molecular  weight  products  result  in 
reduction  of  the  amount  of  energy  "wasted"  as  propellant  product  kinetic  energy  for  a 
given  imparted  p'^ojectile  kinetic  energy.  Looked  at  another  way,  low  molecular  weight 
products  result  in  reduced  pressure  differential  from  the  gun  breech  to  the  base  of  the 
projectile  in  the  gun  banel.  Since  maximum  allowable  breech  pressure  is  an  important 
gun  design  parameter,  this  reduction  in  pressure  gradient  results  in  increased  application 
of  force  to  the  projectile  base  for  a  given  gun  design."  “ 
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The  superiority  to  conventional  chemical  propulsion  goes  beyond  simply  reducing  the 
molecular  weight  and  pressure  gradients  of  the  plasma.  Higher  propellant  gas  energies, 
improved  ballistic  pressure-time  profile  tailoring,  higher  sound  speeds  and  lower  gas 
temperatures  are  all  benefits  of  using  hydrogen.  A  way  to  measure  the  effectiveness  of 
projectile  acceleration  with  various  propellants  is  to  compute  the  impetus.  The  impetus  is 
defined  as  "an  index  generally  used  by  ballisticians  to  describe  the  potential  of  a 
propellant  to  do  work." '  It  is  defined  as: 


I^RT/m  (2  1) 

R  is  the  universal  gas  constant,  T  is  the  adiabatic  flame  temperature  and  m  is  the 
molecular  weight  of  the  propellant  gas.  By  considering  this  equation,  it  is  apparent  that 
to  maximize  the  impetus,  either  the  temperature  can  be  increased,  or  the  molecular  weight 
can  be  decreased.  "For  conventional  propellants,  T  and  M  are  determined  by  the 
propellant  chemistry.  In  the  ELGG,  however,  the  temperature  can  be  increased  by 
introducing  more  electrical  energy  into  the  working  fluid.  This  is  an  enormous  potential 
advantage," '  The  following  table  demonstrates  the  differences  in  the  impetus  for  JA2  (an 
Army  artillery  chemical  gun  propellant),  hydiogen  gas,  water,  and  carbon  dioxide.  All  are 
calculated  at  a  temperature  of  3,450K. ' 
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Table  2.2  Impetus  calculations  for  various  propellants ' 


Material  (Species) 

Molecular  Weight 
of  Species 

Impetus  (J/g) 

JA2 

24.9 

1,142 

2 

14,342 

ILO 

18 

1,594 

CO, 

44 

652 

— - =lli 

At  the  flame  temperature  of  3, 4501s  the  impetus  of  pure  hydvogen  would  be  nine 
times  greater  than  pure  water,  and  more  interestingly,  it  is  ovsr  twelve  times  g’^eater  than 
that  available  from  JA  2.  Clearly,  lowering  of  the  molecular  weight  of  propulsion  gases 
has  a  significant  beneficial  effect  on  impetus.  "To  tne  evtent  trial  ELGG  cait  take 
advantage  of  lower  molecular  weight  propelling  gases,  therefore,  it  c'ruld  offer  sig'iifican; 
increases  in  ballistic  performance. '  The  the  ELGG  could  favorably  increase  the  kinetic 
energy  distribution  between  the  projectile  and  the  propellant.  This  is  accomplished  by 
increasing  the  impetus  with  a  low  molecular  weight  working  fluid.  T.he  kinetic  eiiergy 
budget  in  a  gun  may  be  described  as; 

KE„,  =  +  KE,„  (2.2) 

KE„  =  I/2(M  V (I  +  C/3M)  (2.3) 

M  is  the  projectile  mass,  V  is  the  projectile  velocity,  and  C  is  the  propellant  charge 
mass.  "With  higher  impetus  gasses  (working  fluids)  the  value  of  C  decreases  for 
equivalent  performance,  altering  the  kinetic  energy  distribution  in  favor  of  the  projectile. 


20 


For  example,  in  conventional  propellants,  C/M  ratios  of  approximately  3  are  required  to 
obtain  velocities  of  2  km/sec.  This  means  that  in  conventional  ballistics  there  is  more 
kinetic  energy  in  the  propellant  gases  than  the  projectile  once  a  velocity  of  2  km/sec  is 
exceeded.  Lower  molecular  weight  ELGG  working  fluids  leading  to  reduced  C, 
therefore,  could  present  a  significant  advantage  at  high  muzzle  velocities." '  Clearly  from 
this  discussion,  hydrogen  is  the  propellant  of  choice  for  these  considerations. 

Unfortunately,  it  is  not  a  simple  matter  of  simply  choosing  a  particular  gas  to  use  as 
the  working  fluid  based  solely  on  a  high  sound  speed  and  low  molecular  mass.  The 
interior  of  the  ELGG  barrel  experiences  rapid  changes  in  working  fluid  temperature, 
pressure  aiid  density  due  to  the  various  heat  loss  mechanisms  and  the  viscous  effects. 

TIus  ciiaJiging  eir/ironmeni  causes  drastic  changes  in  the  thermodynamic  properties  of  the 
working  fluid,  as  well  as  the  optical  thickness  associated  with  radiation  heat  transport. 

For  example,  the  ratio  of  specific  heats  for  hydrogen  can  vary  by  an  order  of  magnitude 
simply  from  a  temperature  change  from  10,000  K  to  7,OOOK.  This  magnitude  (or  larger) 
of  tempeiature  decrease  can  be  expected  throughout  the  length  of  the  ELGG  barrel, 

Additionally,  the  changes  in  the  thermodynamic  properties  of  the  working  fluid  greatly 
affect  the  flow  clmracteristics  within  the  ban'el  region,  allowing  for  the  possibility  of 
shock  development  and  thus  choked  flow  prior  to  the  exit  of  the  barrel.  The  boundary 
layer  growth  and  development  within  the  tube  is  highly  dependent  on  the  initial  conditions 
of  the  fluid  entering  the  barrel,  as  well  as  the  conditions  existing  at  any  point  in  the  flow 
up  to  the  barrel  exit. 
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2.3  -  POWER  SOURCE  FOR  THE  ELECTROTHERMAL  LIGHT  GAS  GUN 


At  the  Aeronautical  and  Astronautical  Research  Laboratory  of  The  Ohio  State 
University  a  gigawatt-level,  quasi-steady,  high  speed  flow  capacitor  bank  and  facility 
commonly  referred  to  as  "Godzilla"  has  been  constructed  .  "A  flexible  arrangement  has 
been  adopted  involving  2,100,  43  microfarad  capacitors,  at  6  kilovolts,  that  allows 
giga-watt  power  levels  (into  a  matched  load)  at  333  kilo-amperes  and  3,000  volts,  for 
1.63  milliseconds,  and  various  lowei  current,  longer  pulsetime  combinations,  (e.g.,  1 1 1 
kilo-amperes  for  4.9  milliseconds)."  ® 

Aluminum  support  racks  were  obtained  and  assembled  (each  rack  or  bank  is  28  feet 


long,  1 4  feet  high,  and  4  feet  deed)  into  three  equal  sections.  "CousU  uctiun  of  tlic  racks 


allowed  seven  equal  horizontal  rows  on  each  side  of  a  rack.  These  rows  are  in  turn 
separated  into  five  sections  supporting  ten  capacitors  each."  ’ 

The  capacitance  per  section  is  calculated  to  be  approximately  3  micro-faradays.  The 
charging  voltage  of  6  kilovolts  and  power  input  using  a  matched  load  of  10’  Watts,  the 
total  transmission  line  impedance  is  9  milliohms,  ’  "The  two-way  transit  time  in  the  line  is 
1.625  milliseconds,  which  corresponds  to  the  release  of  1 .625  Mega  Joules  (MJ)  at  1 
Giga-watt.  The  total  current  is  333  Kilo-amperes  or  1 1 1  Kilo-amperes  per  line."  ’ 

By  disconnecting  two  banks  from  the  system,  1 1 1  Kilo-amperes  can  be  provided  into 
a  27  milliohm  load  for  4.875  milliseconds.  With  two  banks  connected,  an  intermediate 
222  Kilo-amperes  can  be  provided  for  3.25  milliseconds.  ’  The  technical  specifications  of 
"Godzilla"  are  well  within  the  energy  requirements  for  the  electrothermal  gun,  and  make 
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it  an  obvious  choice  for  electrical  energy,  at  least  for  short  duration  testing.  For 
development  of  the  heat  transfer  studies  in  following  chapters,  the  technical  specifications 
of  "Godzilla"  will  be  used. 


CHAPTER  HI 


APPROACH 


3.1  -  DESCRffTION  OF  PLASMA  STATE 

Before  continuing  the  study  of  the  plasma  contained  within  the  electrothermal  gun, 
the  static  and  equilibrium  properties  of  the  plasma  must  be  established.  For  example,  the 
population  of  ions  and  electrons  that  exist  at  any  specified  total  pressure  and  temperature 
must  be  calculated  This  can  be  accomplished  by  beginning  with  the  assumption  that  the 
gas  is  in  thermodynamic  equilibrium,  and  then  move  on  to  conditions  that  do  not  satisfy 
this  assumption. 

Using  the  Saha  equation  which  determines  the  fraction  of  particles  that  have  been 
ionized  in  tl'.ermal  equilibrium,  we  can  determine  the  degree  of  ionization  of  the  plasma. 
If  the  ionization  is  high,  then  a  significant  amv»unt  of  energy  ij  used  for  this  process  and 
represents  a  loss  that  could  otherwise  have  bei  ‘  used  for  projectile  motion.  The 
following  is  a  variation  of  the  familiar  Saha  equation:'® 
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(^7)2exp[™]  (3.2) 


where: 


u. 


A)j 


k 

h 

T 

X 


=  Partition  mnction  for  ionized  hydrogen  ( -  I) 

=  Partition  function  for  ground  state  hydrogen  ( -•  2) 
=  Ionization  energy  for  hydrogen  (  =  13  .6  eV) 

-  Total  pressure  of  the  ( in  PsicaJ) 

Rcjt  mass  of  elearon  ( =  9. 11  x  10‘^'  Kg) 

-  Planck's  Constant  (  =  1.38  x  10  “  J/K) 

=  Boltznuui  Constant  ( =  6.63  x  10  ^  J  sec) 

=  Total  single  temperature  ( in  K) 

-  Degree  of  ionization 


"In  order  to  compute  the  degree  of  ionization,  the  partition  function  ha.s  to  be  known. 
For  this  we  need  the  statistical  weights  of  the  different  states  of  excitation,  which  are 
given  by  quantum  mechanics.  Since  the  highei-  states  contribute  little  to  the  partition 
function,  we  may  approximate  it  by  the  weight  of  the  ground  state,  go  o  =  2,  while  for 

iofized  hydrogen,  U;  "  I  " Since  11^  incr^ftses  with  temperature  and  decreases  with 


pi 


#4  va^v  uaiM  vaav  \/a  iViitxAAUVii  aiiva 


*  w 


with  temperature  and  decrcTiSes  with  the  gas  pressure.  "This  can  be  easily  understood: 
with  increasing  Umpeiatuie  the  collisions  become  more  violent,  and  the  process  of 
'kicking  off  the  elsctrons  from  the  atoms  more  frequent.  If,  on  the  other  hand,  the 
temperarorc  is  kept  constant  but  the  pressure  increases,  then  the  probability  grows  that 
tlie  ion  .nceis  an  electron  and  recombines." '® 


3.2  -  DESCRIFflON  OF  ENERGY  LOSS  MECHANISMS 


As  a  starting  point  to  understanding  the  heat  transfer  process  within  the  electrothermal 
gun,  it  is  helpfijl  to  model  the  gun  as  a  control  volume  and  illustrate  the  first  law  of 
thermodynamics.  The  following  is  an  illustration  of  the  model: 


The  study  of  heat  traatfcr  concerns  the  rate  term  in  the  above  figure.  There  ai  e 
three  general  types  of  heat  transfer  mechanisms  operating  within  the  elearothermal  gun 
barrel  region.  The.se  are  1)  them'.al  radiation  from  the  plasma  to  the  walls  and  projectile, 
2)  thermal  conduction  from  the  plasma  to  the  walls  and  projectile,  and  3)  forced 
convection  due  to  the  motion  of  the  plasma. 

Conduction  is  caused  by  the  fact  that  molecules  in  higher  temperature  regions  have 
higher  khietic  energies  than  molecules  in  adjacent,  but  cooler  temperature.  These  higher 
energy  molecides  transmit  some  of  their  energy  to  the  cooler  adjacent  molecules  through 
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collisions  Conduction  can  occur  both  internally  (among  particles  making  up  the  fluid),  or 
by  contact  with  another  substance,  such  as  a  ELGG  barrel  wall. 

In  convection,  heat  is  carried  from  a  higher  to  a  lower  temperature  region  within  a 
moving  fluid.  This  is  a  macroscopic  effect,  and  occurs  in  one  of  two  main  ways.  The 
first  is  natural  convection,  which  is  caused  simply  because  of  density  differences  caused 
by  temperature  gradients.  The  second  is  forced  convection,  which  is  caused  by  some 
external  mechanism,  such  as  an  electric  arc  heating  a  fluid  to  high  temperatures,  then 
expanding  rapidly  down  the  barrel. 

In  contrast  to  conduction  and  convection  which  rely  on  collisions  or  bulk  movement 
of  the  fluid,  radiation  allows  heat  energy  to  travel  over  distances  without  the  aid  of 
contact.  Radiation  operates  by  the  expelling  and  recapturing  photons. 

In  addition  to  the  loss  of  energy  due  to  the  various  heat  transfer  mechanisms 
detailed  above,  viscous  effects  represent  another  form  of  energy  loss.  The  forces  of 
friction  act  as  a  shear  stress  between  the  walls  and  the  moving  fluid,  which  in  essence 
drags  the  fluid  velocity  near  the  wall  from  the  freestream  value  in  the  mid-line  of  the  tube 
to  zero  at  the  face  of  the  wall.  This  velocity  gradient  forms  a  boundary  layer  in  which 
heat  transfer  by  conduction  takes  place,  as  well  as  serves  to  decrease  the  available  area 
for  the  freestream  to  travel  along  the  interior  of  the  tube. 


3.3  -  DEVELOPMENT  OF  MODEL 


Development  of  a  model  for  the  heat  loss  in  a  fluid  traversing  a  gun  barrel  requires  the 
adoption  of  several  assumptions  in  order  to  make  the  equations  of  motion  and  equation  of 
state  tractable.  The  purpose  of  this  model  is,  in  effect,  to  provide  insight  into  the  actual 
processes  occurring  in  the  barrel  region  of  an  ELGG.  These  equations  and  the  model  to 
be  developed  will  be  used  in  order  to  verify  and  compare  results  given  in  MACH2,  a  2 
1/2  dimensional  magnetohydrodynamic  simulation  code  for  problems  with  complex 
shapes  and  areas  of  operation. 

To  develop  this  model,  the  method  of  a  one  dimensional  compressible  flow  through  a 
constant  area  duct  will  be  used.  Thus,  it  will  be  assumed  that  the  flow  variables  such  as 
pressure  or  temperature  vary  in  the  flow  direaion  only.  Consider  the  flow  through  a  one 
dimensional  control  volume  as  shown  in  the  following  figure: 


Figure  3.2  Flow  through  control  volume.  " 
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Consider  this  region  as  one  of  heat  addition  or,  in  this  case,  heat  extraction.  What  this 
is  suggesting  is  that  flow  properties  change  as  a  function  of  the  position  aJong  the  flow 
direction.  Since  the  flow  is  assumed  to  be  one  dimensional,  the  flow  properties  are 
constant  on  either  side  of  the  region  of  heat  loss  bounded  by  the  control  volume.  The 
areas  are  equal  since  we  are  dealing  with  a  constant  area  barrel.  It  v.ill  be  assumed  that 
the  flow  is  steady,  requiring  all  time  derivatives  to  be  zero. 

With  these  assumptions,  the  following  equations  of  motion  re.sult." 

CONTINUITY;  piAiU}  =  P2A2U2  (3.3) 

TA.f-  P  ,  a.  —  P.»  X  /-i  A\ 

ENERGY;  hi  +  j -h  g  =  h2  +  J  (3.5) 


kv  Vliv  VV{U<*VlVUO  \Jl  lllWilUII, 


II  wtu  uw  ooduiii^u  iiiai  me;  iiuiud  udcu  cue; 


calorically  perfect  and  are  represented  by  the  following  equation  of  state;" 


EQUATION  OF  STATE;  P  =  pPT  (3.6) 

In  order  to  complete  the  development  of  the  one  dimensional  model,  the  first  and 
second  laws  of  thermodynamics  must  be  discussed.  Consider  a  system  which  has  a  fixed 
amount  of  heat  added  for  an  incremental  amount  of  work  done  on  the  system  by  the 
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surroundings.  It  is  assumed  that  the  processes  within  this  model  are  isentropic  unless 
otherwise  specified.  The  change  in  specific  internal  energy  is  thus  written  as,  using  v,  to 
be  the  specific  volume:" 

de-dq-pdv  (3.7) 

The  second  law  of  thermodynamics  is  important  in  that  it  describes  the  direction  a 
process  will  proceed.  This  directed  movement  of  a  process  is  termed  entropy,  and  is 
defined  as; 


ds  —  y,  +  ds^irrev)  (^-8) 

From  this  equation,  5q  is  an  incremental  amount  of  heat  added  reversibly  to  the  system, 
and  dS(^^  is  a  dissipative  loss  which  is  not  reversible,  such  as  fiiction.  This  suggests  that 
any  dissipative  process  increases  the  entropy,  and  we  can  see  that  at  best,  entropy  merely 
remains  constant.  For  an  adiabatic  process,  5q  is  equal  to  zero. 

Using  relations  for  thennally  perfect  gasses  in  place  of  de  and  dh,  we  find:" 

S2~S\=  Cp\n{^)  -  R  \n(~)  (3.9) 

or 

52  -  Si  =  CvlnC^)  +  R ln(-^) 


(3  10) 
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These  two  equations  allow  the  calculation  of  the  change  in  entropy  between  two  states 
of  a  calorically  perfect  gas  in  terms  of  either  pressure  or  volume  and  temperature.  For  an 
isentropic  process  (by  definition,  ds=0),  the  entropy  is  constant  throughout  the  system. 
From  this,  important  and  useful  relations  are  developed  and  are:" 


Pi 


(3.11) 


This  set  of  relations  is  valuable  in  that  it  relates  the  flow  properties  of  pressure, 
temperature,  and  density  at  different  locations  in  the  flow.  Using  the  stagnation 
conditions,  we  can  relate  the  stagnation  enthalpy  to  the  static  flow  properties:" 


ho  =  CpTo 

(3.12) 

CpT+^  =  CpTo 

(3.13) 

In  this  equation,  T^  is  the  stagnation  temperature  and  would  exist  if  we  were  able  to 
adiabatically  bring  the  fluid  element  to  rest.  The  following  useful  relations  can  be 
obtained  by  using  combinations  of  previous  equations:" 


(3.14) 


(3.15) 
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^  =  (3,16) 

These  relations  do  not  require  that  the  flow  be  adiabatic  or  isentropic,  which  allows 
them  to  be  used  throughout  the  muzzle  section  of  the  ELGG  and  for  various  conditions 
of  the  flow.  However,  if  the  general  flow  field  is  isentropic  throughout,  then  T„,  P„,  and 
Po  are  constant  at  all  locations  in  the  flow.  The  change  in  entropy  is  related  to  the  change 
in  stagnation  pressure  from  the  following  equation  " 

S2~Si  =-R\n(^)  (3.17) 

The  change  in  specific  energy'  is  given  by  the  Kugonoit  equation  as  follows;" 

e2-ei  =— (3,18) 

From  these  equations  of  motion,  the  equation  of  state,  thermodynamic  relations,  and 
derived  fluid  flow  relations,  the  model  is  now  capable  of  being  tailored  for  use  in 
calculation  of  the  various  heat  transfer  modes,  and  can  be  adapted  for  viscous  effects. 


3.4  -  VISCOUS  EFFECTS  MODEL 


In  order  to  develop  a  model  to  account  for  the  frictional  effects  of  flow  through  a  gun 
barrel,  assumptions  need  to  be  made  which  will  make  the  equations  analytically  solvable. 
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Consider  a  one  dimensional  flow  of  a  compressible  fluid  in  a  pipe.  It  is  assumed  that  the 
flow  is  steady,  adiabatic,  and  shockless.  The  assumption  of  shockless  flow  is  important 
and  will  be  discussed  shortly.  If  this  fluid  is  viscous,  the  friction  acting  between  the 
moving  fluid  and  barrel  will  cause  the  flow  properties  to  vary  as  they  travel  down  the 
tube. 

If  this  frictional  effect  is  modeled  as  shear  stress  at  the  wall  acting  on  the  fluid  with  an 
assumed  uniform  set  of  properties  across  the  cross  section,  then  the  following  set  of 
equations  of  motion  may  be  used;" 

CONTINUITY:  PlUiA  i  =  P2U2A2  (3.19) 

MOMENTUM;  Pi  -h  piU^i  =  P2  +  P2uj  +  ^  f  Ty,dx  (3.20) 

where:  —  |  X^dx  is  evaluated  from  0  to  L  and  accounts  for  frictional  losses, 

2  2 

ENERGY;  /li  +^  =  h2+~  (3,21) 

EQUATION  OF  STATE;  P  =  pRT  (3.22) 

The  modification  of  the  momentum  equation  is  due  to  the  shear  stress,  ,  acting  on 
the  surface  of  the  gun  barrel,  thus  imparting  a  surface  force  in  the  integral  formation  of 
the  momentum  equation.  The  frictional  effect  added  to  the  momentum  equation  can  be 
rewritten  by  considering  the  total  pressure  drop  over  a  distance  dx: 


dp  =  -Aj^dx-pVdV 


(3.23) 


33 


where:  dp  =  total  pressure  change 

f  ”  friction  factor  (dimensionless) 

p-  density 

V  ~  velocity 

D  ==  hydraulic  diameter 

dx  -  length  of  tube 

dV  -•  incremental  velocity  change 


By  definition,  the  fiictiou  factor,  f,  is;" 


(3.24) 


Using  the  constant  area  duct  assumption,  and  the  perfect  equation,  it  can  be 
shown  that:" 


_  f 

D  ^3  25) 


The  fiictional  faaor  depends  greatly  cn  the  flov^'  Reynolds  number  as  well  as  the 
Mach  number  and  surface  roughness,  among  others.  The  effect  the  Reynolds  number  has 
is  to  determine  whether  the  flow  is  laminar  or  turbulent.  The  sttrface  rouglmcss  is  a 
measure  of  the  hydraulic  smootlmess  of  the  interior  of  the  gun  oaiTei.  For  this  thesis,  it 
will  be  assumed  that  the  surface  is  hydraulically  smooth,  meaning  that  the  surface 
roughness  is  less  than  .001  times  the  interior  diamctei  of  the  tube." 
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If  the  inlet  flow  is  supersonic,  friction  will  tend  to  drive  the  Mach  number  to  I .  If  the 
flow  is  subsonic,  the  friction  will  tend  to  bring  the  mach  number  up  to  1 .  "Whether  or  not 
[the  Mach  number]  will  become  unity  depends  on  the  extent  of  the  friction,  which  is 
actually  dependent  on  the  length  of  the  pipe."'* 

If  the  actual  pipe  or  barrel  length  exceeds  the  theoretical  length  necessary  for  the  inlet 
Mach  number  to  be  driven  to  Mach  1,  choking  will  occur,  and  if  the  inlet  flow  is 
supersonic,  then  shock  phenomena  will  develop  and  may  progress  upstream  within  the 
barrel.  When  choking  occurs,  the  flow  must  readjust  itself  since  the  mass  flow  has 
reached  a  maximum  and  cannot  be  increased  further.  There  are  two  main  ways  in  which 
this  can  be  eliminated.  For  supersonic  flow,  the  shock  wave  formation  can  be  accepted, 
and  allow  it  to  "correct"  the  flow  mlel  conditions,  sufiering  a  loss  in  eflTiciency,  or  tlie 
inlet  conditions  can  be  adjusted  in  advance,  allowing  the  Mach  number  to  at  most  become 
sonic  at  the  exit  of  the  barrel.  This  would  represent  the  maximum  mass  flow  possible 
\vithout  shock  formation  (supersonic  flow)  and  thus  provide  the  maximum  gas  enthalpy 
which  would  theoretically  be  able  to  provide  the  best  possible  projectile  acceleration. 

Note  that  Mach  one  can  be  approached  from  either  supersonic  or  subsonic  inlet 
conditions. 

The  following  figure  shows  a  graphical  analysis  of  the  effect  of  various  frictional 
factors  on  the  theoretical  maximum  length  of  pipe  for  various  inlet  Mach  numbers.  The 
full  theoretical  length  analysis  and  calculated  data  is  located  in  appendbt  J. 
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A  ^  ^  V*  Via  VllW'l  IWll^kll  va  ^l^V 


The  final  consideration  is  the  selection  of  the  average  fnction  factor  which  depends  on 
whether  the  flow  is  laminar  or  turbulent.  "In  almost  all  practical  cases,  the  flow  is 
turbulent,  and  the  variation  of  f  must  be  obtained  empirically.  Anderson  further  suggests 


that  for  our  purposes,  it  is  reasonable  to  assume  an  approximate  constant  value  of  .005"'^ 
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Experiments  by  Keenan  and  Neumann  were  undertaken  in  the  hopes  of  determining 
the  turbulent  friction  factor  for  flow  in  pipes.  "They  report  that  for  supersonic  flow  of  air 
in  a  tube,  the  friction  coefficient  is  influenced  as  usual  by  the  Reynolds  number  and  a 
dimensionless  parameter  L/D,  in  which  L  is  the  length  of  the  test  section  and  D  is  the 
inside  diameter  of  the  tube."'*  They  also  discovered  that  at  high  velocities,  the  Mach 
number  is  also  an  influencing  factor.  "Keenan  and  Neumann  studied  friaional  phenomena 
in  smooth  pipes  for  Mach  numbers  from  0.22  to  3.87  and  for  Reynolds  numbers  from 
1x10^  to  8.7xl0^  They  found  that  at  supersonic  speeds,  the  results  were  influenced 
greatly  by  the  presence  of  shock  waves.'"' 

"They  concluded  further  that,  when  the  ratio  L/D  was  greater  than  50.0,  the  friction 
factor  was  approximately  equal  to  the  friaion  factor  for  incompressible  flow  for  the  same 
Reynolds  number.  They  found  also  that  supersonic  flow  can  rarely  be  maintained  if  the 
L/D  ratio  is  greater  than  50.0."'^  Thus  it  is  determined  that  Mach  number  acts  as  a 
limitation  to  the  possible  range  of  L/D  values  . 

Keenan  and  Neumann  conclude  that  the  coefficient  of  friction  is  not  well  understood 
for  turbulent  flow,  mostly  due  to  the  fact  that  the  friction  effects  can  vary  from  place  to 
place  in  turbulent  flow.  This  necessitates  the  consideration  of  an  apparent  or  average 
friction  factor  rather  than  a  local  one.  They  suggest  that  for  most  supersonic  flow,  a 
value  for  the  friction  factor  between  .002  and  .003  in  hydraulically  smooth  pipes  is 
consistent  with  their  experimental  data.'^  In  fact,  values  of  the  friction  factor  for  similar 
systems  have  been  shown  to  range  from  .002  to  .005,  so  this  is  the  range  of  values  that 


will  be  used  for  this  thesis. 


37 


3.5  -  CONDUCTIVE  AND  RADIATIVE  HEAT  TRANSFER  MODEL 


These  two  different  heat  transfer  modes  have  been  grouped  together  for  an  important 
reason.  The  basic  model  to  be  discussed  is  based  on  a  bulk  fluid  flow  assumption,  which 
is  equally  valid  for  conduction  and  radiation,  except  where  noted.  The  heat  dissipation 
term,  q,  to  be  introduced  in  to  the  energy  equation  can  be  either  from  conduction, 
radiation,  or  a  combination  of  both. 

In  a  similar  manner  to  that  used  in  the  development  of  the  viscous  effects  model,  a 
model  is  developed  for  the  loss  of  heat  energy  from  the  bulk  working  fluid.  For  heat 
transfer  calculations,  note  that  inside  the  control  volume  illustrated  in  figure  3.2, 
something  is  happening  which  causes  the  flow  properties  in  region  2  to  be  different  to 
those  in  region  1 .  This  is  caused  by  the  loss  of  heat  energy  due  to  the  high  temperature  in 
this  region  in  relation  to  the  barrel  wall,  which  is  assumed  to  be  at  a  constant  temperature 
of  300K.  By  losing  heat  energy  by  conduction  or  radiation  (or  a  combination  of  both), 
the  enthalpy  and  specific  internal  energy  are  decreased  and  an  increase  in  entropy  results. 

The  following  one  dimensional  flow  equations  are  similar  to  those  for  fiiction,  except 
there  is  not  an  extra  term  in  the  momentum  equation  (viscous  effects  are  neglected)  and 
the  energy  equation  now  contains  the  quantity,  q.  This  is  the  heat  loss  per  unit  mass. 

CONTINUTTY  Pi  1  =  P2W2^2  (3.26) 


MOMENTUM 


P\  +  plWj  =-P2  +  P2W2 


(3.27) 


ENERGY 


/?  1  4-  Y  +  ^  -  /l2  Y 


(3.28) 


For  £  calorically  perfect  gas,  it  is  possible  to  achieve  an  anal>lical  solution  vviiich  in 
many  ways  is  similar  to  that  described  for  friction.  However,  there  are  some  importanv 
difrerences.  Unlike  fiiaion,  heat  can  be  lost  or  gained,  and  the  process  can  be  reversed 
whereas  friction  only  results  in  a  loss  of  enei'gy.  For  heat  loss,  which  will  be  the  nature  of 
the  flow  through  the  ban  el,  the  Mach  number  diverges  away  from  a  Mach  =  1  condition. 

1  has,  unlike  the  eftccts  caused  by  friction,  choking  of  the  flow  will  not  occur  because  of 
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3.6  -  CONDUCTIVE  AND  CONVECTIVE  HEAT  TRANSFER 

The  heat  exchange  between  a  wall  and  a  fluid  whose  motion  has  been  caused  by 
external  forces  is  called  forced  convection  heat  transfer.  The  resistance  to  heat  transfer  is 
concentrated  primarily  in  a  thin  layer  immediately  adjacent  to  the  surfece  of  the  boundary. 
The  heat  transfer  is  determined  by  the  intensity  of  the  conductive  and  convective 
transport  within  this  boundary  layer.  Thus,  the  heat  transfer  coefficient  is  determined 
essentially  by  the  thickness  of  the  boundary  layer.  The  boundary  layer  thickness  is 
dependent  on  the  type  of  fluid  flow  which  occurs  along  the  surface.  With  forced 
convection,  the  differential  equations  which  describes  the  fluid  and  heat  flow  processes 
yield  the  result  that  the  flow  field  and  hence  me  development  of  the  boundary  layer  are 
independent  of  the  heat  transfer  if  the  properties  of  the  fluid  are  temperature  independent. 
However,  in  general,  all  properties  of  the  fluid  are  temperature  dependent,  but  as  long  as 
the  internal  temperature  differences  are  not  extreme,  the  temperature  dependence  of  the 
property  values  can  be  neglected.  For  this  investigation,  it  will  be  assumed  that  the 
temperature  dependence  can  be  accounted  for  by  use  of  tabulated  thermodynamic  data 
provided  in  appendbc  D 

For  the  development  of  a  simplified  model  to  describe  the  convective  and  conductive 
heat  transfer  within  the  ELGG,  many  fluid  properties,  such  as  mass,  viscosity,  and 
specific  heats,  have  been  taken  into  account  and  actual  experimental  data  is  used 
v/henever  possible.  Density  is  also  included  as  a  property  in  these  calculations.  For  a 
gas,  density  is  highly  pressure  dependent,  especially  at  high  velocities.  Generally,  the 
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limit  for  which  density  becomes  highly  pressure  dependent  is  approximately  one  third  of 
the  sonic  velocity. 


3.6.1  -  BOUNDARY  LAYER  DISCUSSION 


There  are  two  concepts  in  fluid  flow  theory  which  are  essential  for  understanding  heat 
transfer.  These  are  the  concept  of  the  boundary  layer  and  the  concept  of  turbulence.  Due 
to  viscosity  effects,  fluids  suffer  ffictiotial  forces  while  in  motion.  However,  real 
situations  can  be  approximated  by  assuming  a  frictionless  fluid.  Later,  the  effects  of 
viscosity  will  be  added. 

There  are  two  kinds  of  forces  that  will  be  considered  in  the  development  of  the  model. 
These  are  forces  of  inertia  and  forces  of  pressure.  Gravity,  centrifugal  force  and 
electrical  forces  will  be  neglected.  Additional  forces  exist  in  real  fluids  which  are  caused 
by  viscosity.  These  are  shear  stresses  between  individual  stream  lines  which  move  at 

HiflPprpnt  v^lr>r*iti#»c  Th#*  ctr#»cc  t  r*an  ooiiation- 


The  factor,  p  ,  is  the  dynamic  viscosity,  p  is  the  fluid  density,  v  is  the  kinematic 
viscosity,  and  the  shear  stress  is  seen  to  be  directly  proportional  to  tlie  velocity  gradient 
normai  to  the  wall.  The  kinematic  viscosity  of  gasses  is,  by  the  ideal  gas  law,  inversely 
proportional  to  pressure.  Since  numerical  values  of  the  viscosity  are  comparatively  small, 
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high  shear  stresses  in  the  flow,  according  to  Newton's  equation,  will  exist  only  where 
large  velocity  gradients  exist. 

There  a^  e  basically  two  different  forms  of  flow  within  the  ELGG  barrel  region;  laminar 
and  turbulent.  In  laminar  flow,  the  individual  streamlines  are  orderly  afid  p?rallel,  while 
the  turbulent  case  sees  the  ^neamlines  intcnningied  in  an  inegular  way.  Convective  heat 
exchange  frotii  the  working  fluid  to  the  barrel  walls  is  enhanced  by  the  fluctuation  and 
turmoil  in  the  turbulent  region.  This,  of  course,  makes  the  convective  heat  transfer  in  a 
turbulent  flow  considerably  higher  than  it  would  otherwise  be  in  a  laminar  flow. 

Turbulcnl  flow  may  also  exist  wthiii  the  boundary  layer. 

If  fluids  of  different  viscosity's  are  considered ,  the  laminar-turbulent  transition  occurs 
at  U.Xj^v,  which  has  been  known  for  a  long  time  as  tlie  Reynolds  number,  R^.  Typically, 
it  can  be  expected  that  the  boundary  layer  will  become  turbulent  when  ti«e  Reynolds 
number  exceeds  approximately  500,000. 

3.6.2  -  DOIINDARY  LAYER  DEVELOPMENT  IN  A  PIPE 

The  flow  conditions  in  a  pipe  in  the  vicinii '  of  the  inlet  are  similar  to  those  on  a  flat 
plate  in  parallel  flow.  The  thickness  at  the  inlet  is  i;cro  and  increases  rapidly  35  the  flow 
travels  down  the  pipe  (or  barrel).  It  will  be  assumed  that  the  inlet  geometry  is  such  that 
the  flow  into  the  ELGG  barrel  for  the  arc  chamber  or  convergent-divergent  nozzle 
proceeds  without  separation. 
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Figure  3 .4  Growth  of  boundary  layers  in  a  pipe 


At  a  certain  distance,  L„  from  the  inlet,  the  boundary  layers  expand  to  fill  the  entire 
pipe.  The  velocity  profile  downstream  from  this  point  has  the  form  of  a  parabola  (laminar 
flow)  or  an  arched  curve  (turbulent  flow).  Further  downstream  the  velocity  profile  does 
not  change  shape.  The  developed  flow  region  will  be  turbulent  if  the  Reynolds  number  is 
greater  than  a  critical  value.  By  defining  the  Reynolds  number  as  the  mean  velocity  over 
the  cross  section  of  the  pipe,  U„,  and  the  tube  diameter,  d,  the  critical  value  under  normal 
conditions  is:'^ 

3,000  (3.30) 

"In  gen  Tal,  transition  point  calculations  is  the  weakest  link  in  the  convective  heat 
transfer  calculations."” 

3.6.3  -  PIPE  FLOW 


Consider  flow  in  a  pipe.  For  two  dimensional  flow  in  the  entrance  of  a  channel  or  a 


tube,  with  circular  cross  section,  the  equation  for  axisymmetric  flow  are  used.  The 
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following  figure  shows  how  flow  develops  in  the  pipe  firom  inlet  down  stream: 


Figure  3 , 5  Developed  region  of  flow 

The  flow  enters  from  the  left  and  progresses  to  the  right.  The  velocity  in  the  core  of 
the  flow  outside  the  boundary  layers  increases  with  increasing  distance  from  the  entrance 
because  the  fluid  mass  flowing  through  any  cross  section  must  be  constant  and  the 
boundary  layer  thickens  as  it  moves  downstream.  For  laminar  flow,  the  ratio  of  the 
entrance  length,  L„  to  the  diameter,  d,  of  the  tube  is  a  function  of  the  Reynolds  number: 

^  =  Q.02$SRed  (3.31) 

is  calculated  by  using  the  mean  vetccity,  U„  ,  of  the  flow  in  the  tube  and  the  tube 
diair’iCter,  d.  Tiic  laminar/ turbulent  transition  point  movci  toward  the  entrance  of  the  tube 
’  as  the  Reynolds  number  increases.  Since  the  turbulent  boundary  layer  giows  more 
quickly  than' the  laminar  one,  the  entrance  length,  L,,  becomes  shorter  as  the  transition 
point  moves  upstream.  "At  the  critical  Reynolds  number  (R^  approximavely  3000)  it  can 
be  assumed  that  the  flow  in  the  entrance  region  is  laminar. The  equation  above  yields 


44 


at  this  Reynolds  number  an  entrance  length  of  approximately  lOOd.  If  the  Reynolds 
number  increases,  then  the  transition  occurs  in  the  entrance  region.  The  results  of 
numerical  analysis  show  that  L,  at  first  decreases  to  between  lOd  and  40d,  depending  on 
the  Reynolds  number,  and  then  starts  to  slowly  increase. 

The  velocity  profile  in  the  fully  developed  turbulent  flow,  up  to  Reynolds  numbers  of 
about  100,000  can  be  represented  by  Prandtl's  equation  (previously  discussed)  with  the 
boundary  layer  thickness  replaced  by  the  radius  of  the  tube.  The  velocity,  Us,  is  now  the 
velocity  of  the  flow  at  the  axis.  "Integration  over  the  tube  cross  sectional  area  yields 
U„=0.82  U,."’  Thus  the  followdng  equations  may  be  used  to  describe  tube  flow;’ 


Xw  = 


(3.32) 


Urn 


2.44 


(3.33) 


63.5 


(3.34) 


Prandtl  has  proposed,  for  Reynolds  numbers  greater  than  100,000,  a  general 
equation  for  the  fiiction  factor,  f,  which  sh')uld  be  used.’  This  equation  appears  as 
follows; 


-^  =  2,0log(^)  +  5.5 


(3.35) 
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r  =  (3-36) 

Here,  V*,  is  called  the  shear  stress  or  fiiction  velocity.  Due  to  the  complicated  nature 
ot  this  equation,  solving  for  the  friction  factor  has  been  discussed  previously,  and  based 
on  the  work  of  Keenan  and  Neumann,  a  friction  factor  between  .002  and  .005  should  be 
used. 

The  previous  development  of  a  turbulent  model  in  a  tube  is  for  smooth  walls  only. 

For  this  thesis,  the  tube  is  considered  to  be  hydraulically  smooth. 

3.6.4  -  HEAT  FLOW  EQUATION 

When  a  body  submerged  in  a  flow  is  heated  or  cooled,  a  temperature  field  is  built  up 
in  the  surrounding  medium  and  a  temperature  boundary  layer  is  formed.  Within  this 
layer,  the  temperature  varies  from  the  value  at  the  wall  to  the  free  stream  temperature. 

At  at  distance  6,  from  the  wall,  the  temperature  reaches  the  free  stream  temperature,  t,. 

|J(r.-0(/ar  =  a(|).  (3.37) 

This  is  the  heat  flow  equation  for  the  thermal  boundary  layer.  By  using  tabulated  and 
experimental  data,  included  in  the  appendices,  the  temperature  and  pre..,sure  dependence 
of  the  specific  heat,  Cp,  the  thermal  conductivity,  k,  and  the  viscosity,  p,  is  accounted  for 
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■at  each  region  of  differing  pressure.  "The  equation  of  state  may  be  used  to  determine  the 
pressure  values  for  gasses."’ 

The  ratio  v/a  is  a  dimensionless  property  which  is  used  a  lot  in  heat  transfer 
equations.  This  is  termed  the  Prandtl  number,  Pr: 

Pr  =  ^-‘■f  (3.38) 


"In  general,  the  Prandtl  number  is  a  function  of  temperature  only.  The  temperature 
dependence  for  gasses  is  rather  small."’  Using  this  fact  and  solving  for  L  yields;” 


( 1  \  75  >  t 

r  X  J  1  '  •  *  ,  f- 

^  -  '  ■n.i. wnere  L~oj  6^ 

1.026{P,}5 


(3.39) 


"Gases  typically  have  Prandtl  numbers  somewhat  smaller  than  1".’  In  this  case,  the 
assumption  that  L  <  1  made  previously  would  be  in  error.  "However,  since  the  smallest 
value  for  gaseous  Prandtl  number  is  about  0.6,  the  largest  value  of  L  is  equal  to  1.18  The 
error  introduced  in  using  the  above  equation  for  such  values  of  L  is  very  small."” 

From  this  discussion,  flow  through  a  tube  can  now  be  considered.  If  the  wall  of  a 
tube  through  which  a  fluid  flows  is  heated,  relative  to  the  fluid  flowing  within  the  tube,  a 
thermal  boundary  layer  builds  up  along  the  walls.  This  layer  grows  in  the  downstream 
direction  and  meets  at  the  axis  of  the  tube  at  some  distance  fl'om  the  entrance.  This  is 
where  the  thermal  starting  region  ends.  The  heat  flux  is  assumed  to  be  independent  of  r. 
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"The  heat  transfer  coefficient  for  the  flow  in  a  tube  is  usually  calculated  on  the  basis 
of  the  difference  between  the  average  temperature  of  the  fluid  and  the  wall  temperature  of 
the  tube."’’  The  average  temperature  is,  by  convention,  defined  as  that  temperature 
which  the  fluid  would  assume  if  it  were  instantly  and  adiabatically  mixed  after  leaving  the 
cross  section  of  the  tube  and  is  called  the  bulk  temperature,  t^. 

Note  that  the  bulk  temperature  change  At^/L  per  unit  length  is  constant.  This 
means  that  the  bulk  temperature  changes  linearly  in  the  axial  direction.  Since  the  heat 
transfer  coefficient  is  constant  in  the  thermally  developed  region,  it  follows  that  the  tube 
wall  temperature  also  varies  linearly  in  this  region. 

Heat  transfer  to  the  barrel  walls  for  thermally  developing  conditions  can  now  be 
discussed  when  the  velocity  profile  is  already  fully  developed.  Graetz  and  Nusselt 
developed  the  calculations  for  the  developing  conditions.  Their  findings  show  two 
important  facts.  First,  the  temperature  profile  is  rectangular  in  shape  in  the  section  where 
heat  transfer  starts,  changes  in  the  flow  direction  because  of  thermal  boundary  layer 
formation  and  growth  along  the  wall.  The  point  where  the  boundary  layers  meet  at  the 
axis  of  the  tube  defines  the  end  of  the  starting  region.  From  this  point  on,  the 
temperature  profile  does  not  change,  it  only  decreases  in  size  gradually  in  the  direction  of 
flow.  Second,  the  Nusselt  number  is  infinitely  large  at  the  beginning  of  the  tube  and 
decreases  until  it  reaches  the  value  3.65.  The  thermal  starting  length  is  calculated  from; 


^  =  O.OSR^Pr 


(3,40) 
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This  equation  gives,  with  good  approximation,  the  entrance  length  for  the  case  of 
simultaneous  flow  and  thermal  development.  "H.  Hausen  has  developed  a  formula  which 
represents  the  results  of  the  theory  of  both  Graet2  and  Nusselt  This  is  as  follows: 


NUd  =  3.65  + 


0.06g(g)/?«xPr 

l+0.04{(|)ft,A)^ 


(3.41) 


The  influence  of  variable  viscosity  is  taken  into  account  by  multiplying  the  right  hand 
side  of  the  equation  by  the  ratio  (Pt,/  where  is  the  viscosity  evaluated  at  the  bulk 
temperature  and  p^  is  the  viscosity  at  the  wall  temperature.  The  properties  in  the 
dimensionless  parameters  of  the  equation  are  then  based  on  bulk  temperature. 


3.6.5  -  FORCED  CONVECTION  IN  TURBULENT  FLOW 


In  turbulent  flow,  heat  is  transferred  by  the  turbulent  mixing  motions  in  addition  to  the 
heat  flow  by  conduction  and  by  bulk  convection  due  to  the  motion  of  the  fluid.  The 
turbulent  fluctuations  become  smaller  as  a  solid  wall  is  approached  and  the  main  features 
of  heat  transport  from  the  flow  to  the  wall  can  be  specified  by  a  simple  model,  assuming  a 
laminar  sublayer  which  exists  immediately  next  to  the  wall.  The  her "  ^  must  pass 
through  this  layer  before  it  is  transported  to  the  wall  by  turbulent  mixing.  The  interaction 
between  the  laminar  sublayer  and  the  turbulent  core  determines  the  heat  transfer  process. 
By  using  relationships  proposed  by  Reynolds  and  Prandtl,  it  is  possible  to  derive  formulas 
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for  heat  transfer  which  require  only  hydrodynamic  iiiformation.  Generally,  the  fluid  layers 
near  the  wall  are  the  main  component  influencing  heat  transfer.  The  velocities  in  these 
layers  are  parallel  to  the  wall  while  the  heat  flow  is  normal  to  the  wall.  It  is  assumed  that 
the  gradient  in  velocity  is  in  the  y  direction  (normtJ  to  the  wall)  only.  Thus  the  heat  flux 
is  in  the  y  direction,  and  the  temperature  change  is  mainly  in  the  y  direction  as  well. 

The  turbulent  shear  stress  may  be  written  as;"'® 

Xt  =  m'{U'  -U)  (3.42) 

and  the  relationship  between  the  heat  transfer  and  shear  stress  in  turbulent  flow 
becomes;'® 


^ f  ~  (3 .43) 

"This  relation  was  first  derived  by  Reynolds  in  1 874  and  is  therefore  called  the 
Re5molds  analogy."'^  The  ratio  of  the  heat  transfer  to  the  shear  stress  has  the  same 
relationship  in  laminar  or  turbulent  flow  when  k/p  -  c^  or  when  Cpp/k  =  =  1 .  A  fluid 

having  a  Prandtl  number  equal  to  1  exhibits  the  same  heat  transfer  relationships  in  laminar 
or  in  turbulent  flow.  For  fully  developed  flow  in  a  tube,  the  free  stream  values  are 
replaced  with  the  mean  velocity,  U„,  and  the  bulk  temperature,  t^. 

Qyy  =  =  R^itb  -  ty.) 


(3.44) 
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R  is  the  resistance  (R=x„A)  and  is  the  heat  flow  at  the  wall.  Note  that  R  is  only 
valid  for  Prandtl  number  very  close  to  1 .  "Since  most  gasses  have  Prandtl  numbers  which 
deviate  from  1  only  slightly,  this  equation  is  very  usefiil  in  obtaining  a  first  approximation 
to  the  heat  transfer  where  the  resistance  is  known. The  resistance  in  a  tube  with 
pressure  drop  Ap  is  R  =  Apd^7t/4  and  the  volume  flowing  through  the  tube  cross  section  is 
V  =  U„dV4.  Thus  the  power  needed  to  maintain  the  flow  is: 


W=V^p^RU„,  (3.45) 

The  ratio  of  the  heat  flow  to  power  is  an  important  relationship  and  demonstrates  that 
the  power  needed  for  a  certain  heat  transfer  shrinks  when  the  mean  velocity,  U^,  is 
decreased.  T!iis  is  given  by: 

%  =  jk(h-tn)  (3.46) 

3.6.6  -  TURBULENT  FLOW  IN  A  TUBE 

In  fully  developed  flow,  the  boundary  layers  have  met  at  the  tube  axis  with  velocity  U,. 
It  becomes  apparent  that  the  heat  transfer  coefficient  described  in  the  last  equation  should 
be  based  on  the  temperature,  t„  along  the  axis.  However,  as  stated  earlier,  the  heat 
transfer  coefficient  in  a  tube  is  based  on  the  bulk  temperature,  t^,  and  the  mean  velocity, 
U„.  "This  is  not  quite  correct  because  of  the  different  definitions  of  t^  and  U^,  but  the 
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error  introduced  in  this  way  is  quite  small."”  Using  this  information,  the  final  formula  for 
heat  transfer  to  the  walls  of  a  tube  with  turbulent  flow  is: 


St  = 


zl 

o.omR*j 

n  5" 

l+1.7R^(P,-I)P/ 


(3.47) 


Note  that  for  the  heat  transfer  equation  as  written,  the  term  on  the  left  hand  side  is 
equivalent  to  NUy(R^  P,)  which  is  the  Stanton  number  (S).  The  shape  of  the 
temperature  profile  for  different  Prandtl  numbers  can  be  derived  from  these  equations. 
Note  that  for  a  fluid  with  Prandtl  number  equal  to  1 ,  the  ratio  of  the  temperature 
differences  is  equal  to  the  ratio  of  velocity  differences.  As  Prandtl  number  increases,  the 
temperature  drop  in  the  laminar  sublayer  begins  to  dominate.  The  convection  turbulent 
heat  flow  to  the  wall  of  a  tube  of  length  L  is:” 


Q  =  hdnL(tw  -  tb)  (3  48) 

The  methodologies  discussed  throughout  this  section  are  important  for  the 
understanding  of  the  actual  flow  problems  inherent  in  an  ELGG  barrel,  and  are  necessary 
to  properly  evaluate  the  data  and  plots  obtained  through  the  MACH2  calculations. 
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3.7  -  RADIATIVE  HEAT  TRANSFER  PROCESSES  AND  HYDROGEN 
OPACITY 

As  discussed  previously,  there  are  three  types  (generally)  of  heat  transfer.  These  are 
heat  convection,  heat  conduction  and  thermal  radiation.  "If  the  temperature  of  a  gas  is 
not  too  high,  and  the  density  is  not  too  low,  the  transfer  of  heat  by  thermal  radiation  is 
usually  negligibly  small  in  comparison  with  that  by  conduction  and  convection."” 
However,  we  are  dealing  with  very  high  stagnation  temperatures  and  the  densities  are 
relatively  high,  such  that  thermal  radiation  must  be  considered.  This  can  be  seen 
graphically  in  the  following  figure;'* 
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Figure  3  6  Regimes  of  energy  transport 

There  are  tL  ee  thermal  radiation  effects  on  the  flow  field  of  a  high  temperature  gas; 
1)  radiation  stresses,  2)  radiation  energy  density,  E,^,  and  3)  heat  flux  of  radiation,  Q^. 
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Usually,  radiation  energy  density'  and  the  radiation  pressure  are  ot  the  same  order  of 
magnitude.  The  following  is  a  definition  of  the  radiation  energy  densiry;'" 


This  radiation  energy  density  should  be  compared  with  the  internal  energy,  U„  =  c^T 
of  the  fluid.  Since  radiation  energy  density  is  proportional  to  T*,  while  internal  energy  is 
proportional  to  T,  a  reduction  in  the  temperature  by  one  order  of  magnitude  will  result  in 
the  radiation  energy  density  becoming  three  orders  of  magnitude  smaller  than  that  of  the 
internal  energy.  By  comparison  with  radiative  energy  density  and  internal  energy,  it  can 
be  determined  when  radiation  energy  density  becomes  important. 

For  the  radiation  equilibrium  condition,  the  radiative  flux  is: 

Qr  =  (3.50) 

where  ct  is  the  Stelfkn-Boltzmann  constant.  The  heat  flux  by  convection  is: 

Qv~UCvT  w  =  the  typical  flow  velocity  (3.51) 


The  radiation  flux  should  be  compared  to  the  heat  flux  by  convection.  Since  is 


proportional  to  T*,  while  Q„  is  proportional  to  T,  at  high  temperatures  such  as  those 
encountered  in  the  ELGG,  may  be  of  the  same  order  of  magnitude  or  higher  than  Q„. 
If,  for  example,  the  density  is  roughly  1/1000  of  air  density  and  temperature  is  lO.OOOK 
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or  higher,  we  should  consider  the  radiation  flux.  At  IGO.OOOK  or  more,  both  radiation 

energy  density  and  radiation  flux  should  be  considered. 

The  thermal  radiation  may  be  expressed  in  terms  of  a  specific  intensity,  1„ 

The  amount  of  radiative  energy  flowing  through  an  area  5a,  in  fi^equency  range  v  to  v 
+  5v  in  direction  of  the  ray  at  an  angle  0  with  normal  of  60  within  solid  angle  6co  in  time 
interval  6t  is  5E„.  The  total  amount  of  energy  radiated  over  the  whole  spectrum  is: 

dE  =  \(^)d'o  =  /  COS  Bdadwdt  (3.53) 

I  I\idx>  this  is  the  integrated  intensity 

If  the  specific  intensity,  I„,  is  kno  wn,  the  flux  of  heat  energy  by  thermal  radiation,  the 
radiation  energy  density  for  the  whole  spectrum,  and  the  radiation  stress  tensor  can  be 
determined. 

At  this  point,  the  process  to  detemiine  the  specular  absorption  coefficient,  K  of 
radiation  is  quite  difficult  and  requires  either  computer  simulation  or  experimentation. 
Fortunately,  this  information  is  available  and  Avill  be  used  extensively  throughout  this 
thesis.  Note  that  the  >'alue  of  the  absorption  coefficient  determines,  among  other  things, 
the  loss  of  specific  intensity  along  the  ray  over  a  distance  in  a  medium,  the  optical 
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thickness,  T,  and  the  heat  flux  due  to  radiation.  The  absorption  coefficient  is  a  function  of 
both  temperature  and  density  of  the  medium,  as  well  as  the  fi'equency  of  radiation. 

The  optical  thickness  indicates  efteclive  length  in  the  absorption  of  radiation.  For  a 
given  length  (s  -  s^),  if  the  optical  thickness  is  large,  the  medium  is  said  to  be  optically 
thick.  If  the  optical  thickness  is  small,  then  we  say  the  medium  is  optically  thin. 

For  most  radiation  heat  flux  considerations,  the  assumption  of  Local  Thermodynamic 
Equilibrium,  or  l,TE,  will  be  used.  This  means  that  at  least  locally,  the  various  processes 
occurring  in  the  fluid  have  h?.d  enough  time  to  relax  to  a  equilibrium  state,  and  it  can  be 
said  that  the  particles  satisfy  the  Boltzmann  and  Maxwell  relations.  Again,  this  is  an 
assumption  for  local  conditions  only.  The  system  itself  can  be  in  a  non-equilibrium  state. 
In  conjunction  with  the  assumption  of  LTE,  the  optical  thickness  of  a  medium  allows  us 
to  choose  different  ways  in  handling  the  thermal  radiation  equation,  which  is:''* 

for  LTE  T^  +  f^=PKC(fiu-/»)  (3  54) 

k(,  =Ku[l-exp^]  (3.55) 

kJ  is  the  spectral  absorption  coefficient  corrected  for  the  effects  of  bound-bound, 
bound-free,  and  free-free  transitions. 

Bound-Bound  transitions  are  simply  a  radiative  change  in  an  atom  or  molecule  from 
one  quantized  energy  level  to  another.  "Since  the  entrgy  difference  between  two 
quantized  energy  levels  is  establisiied  by  internal  molecular  stascttire,  the  emission  and 


absorption  coefficients  for  bound-bound  transitions  aie  sfiarply  peaked  functions  of  the 
photon  energy  hv.  For  this  reason,  bound-bound  transitions  are  said  to  give  rise  to 
emission  and  absorption  lines.  In  a  monatomic  gas,  such  transitions  are  between 
electronic  states  of  the  atom 

Bound-Bound  transitions  resulting  in  absorption  and  emission  lines  have  natural  line 
widths  due  to  the  distribution  of  emitted  and  absorbed  energy  over  a  narrow  range  of  hv. 
The  naturd  line,  width  has  an  associated  uncertainty  v-nich  is  further  coniplicated  tiy  line 
broadening  due  to  Doppler  broadening  and  collision  broadening  Doppler  broadening  is 
the  result  of  thermal  motion  of  the  atoms  or  molecules,  and  collision  broadening  results 
from  a  perturbation  of  the  radiating  atom  by  encounters  with  ether  atoms. 

Bound-Free  absorption  occurs  with  the  interaction  of  radiation  with  an  atom  or 
molecule  causing  photoionization  or  photodissociation  of  the  atom  or  molecule.  The 
reverse  process  is  called  radiative  r  ecombination.  This  results  in  the  removal  of  an 
electron  from  it's  bound,  quantized  state  by  a  photon  of  precisely  the  correct  energy.  This 
ionizes  the  atom,  resulting  in  a  free  electron.  At  this  point,  the  electron  can  take  on  any 
value  of  kinetic  energy,  which  makes  the  absorption  coefficiem  for  photoionization  a 
continuous  non-zero  fu.iction  of  photon  energy  for  all  values  of  photon  energy  that 
exceed  t'.ie  ionization  potential  of  the  atom.  (Vincienti  and  Kruger)  The  process  of 
photodissociation  and  radiative  recombination  is  similar. 

Free-Free  transitions  result  also  in  continuous  emission  and  absorption  which  are  both 
associated  with  changes  in  energy  of  the  free  electrons  in  a  gas.  This  process  is  also 
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called  Bremsstrahlung,  and  occurs  when  a  free  electron  interacts  with  the  field  of  a 
positive  ion. 

Getting  back  to  the  method  of  solving  the  radiative  heat  equation,  the  idea  of 
opt'cal  thickness  can  beused  to  simplify  the  solution  process.  I’bst,  frit  opacity  of  the 
medium  needs  to  be  calculated,  which  will  consist  of  several  different  species  piesent  in 
the  medium.  For  hydrogen,  R.  W.  Patch  has  identified  seven  specii^s  which  occur  at 
different  conditions  which  contribute  significantly  to  the  opacity  of  tlie  medium.  He 
considered  a  system  from  JOOK  to  40,000K  and  pressures  (lela  ing  to  demsities)  from  1  to 
1,000  atmospheres  The  specie.s  he  considered  v/ere;“ 


K  -  hydrogen  atouts 

FT  -  protons 

e'  -  free  electrons 

Hj  -  hydrogen  molecules 

If  -  negative  hydrogen  ions 

H,*  -  hydrogen  diatomic  molecular  ions 

H3*  -  hydrogen  triatomic  molecular  ons 


By  treating  the  plasma  as  an  ideal  gas  except  for  coulomb  interactions  between  free 
charged  particles,  considering  only  singly  excited  electronic  states  of  Hj,  non-degenerate 
free  electrons,  and  assuming  the  photon  pressure  was  negligible,  among  other 
assumptions.  Patch  was  able  to  create  several  tables  and  figures  shov/ing  the  number 
densities  of  the  species  as  a  function  of  temperature  and  pressure.  His  work  includes 
fifteen  different  chemical  reactions,  takes  into  account  bound-bound,  bound-free,  and 
free-free  transitions,  and  corrects  for  varying  pressures  and  tevriperatures.  This 


information  led  to  the  calculation  of  the  Planck  and  Rosseland  opacities  which  are  used  in 
the  solution  process  of  the  radiative  heat  transfer  equation.  The  data  presented  by  Patch 
is  located  in  appendix  D. 

Vincienti  and  Kruger  show  that  the  time  derivative  of  the  radiative  heat  transfer 
equation  car>  be  neglected  by  considering  the  characteristic  speed  and  distances  within  a 
system  to  that  of  the  speed  of  light.  They  show  that  for  most  systems,  this  ELGG 
included,  the  time  derivative  can  be  ne^ected  without  loss  of  information.  There  are  two 
main  ways  in  which  the :  atiiative  heat  transfer  equation  can  be  made  relatively  easy  to 
solve.  These  are  the  optically  thin,  governed  by  Planck's  law,  and  the  optically  thick 
cases,  governed  by  the  Rosseland  difiusion  equation. 

The  optically  thin  case  assumes  that  the  spcc  fjc  intensity,  is  much  less  than  the 
Planck  function,  or  otherwise  known  .ri  the  equilifcu  o  va!t’.e  \(T).  When  this  case  is 
met,  the  gta  is  said  to  be  emission  dominated,  la  other  vvx,.'ds,  vory  little  absorption 
occurs.  If  this  assumption  holds  true  isotropically  a  ■»!.’  for  all  fteqi  cucies,  the  heat  flux 
will  depend  only  or.  local  conditions  and  will  be  equal  to  the  total  x  olun'ctric  rate  of 
spontaneous  emission.  In  this  case,  the  Planck  absorption  coefficient  ic  u  -ed,  or  i.^oie 
accurately,  the  Planck  mean  absorption  coefficient,  Op,  and  it's  ii  verse,  the  .VTvnck  mean 
free  path,  The  Planck  function  and  associated  heat  flux  are  is  eVS  follows; 

Planck  Function:  ^v(7)  =  II  ~  ^  -'5} 

Planck  heat  flux:  ^  =  4ap<jT^  (3.57) 


59 


The  optically  thick  case  occurs  when  the  specific  intensity,  I„,  is  of  the  same  order  of 
magnitude  as  the  Planck  function,  or  otherwise  known  as  the  equilibrium  value  B^(T).  In 
this  case,  the  medium  is  considered  to  be  optically  thick,  and  is  absorption  dominated. 
Vincienti  &  Kruger  demonstrate  that  in  the  case  of  optically  thick  medium,  isotropic 
radiation,  and  for  the  case  that  all  fi'equencies  are  accounted  for,  the  following  Rosseland 
diffusion  approximation  can  be  obtained:'® 


dT 

~  2*aft  dx, 


(3.58) 


where  is  the  Rosseland  mean  absorption  coefficient,  and  it's  inverse  is  the 
Rosseland  mean  free  path. 

At  this  point,  it  is  important  to  be  able  to  understand  which  regime  of  operation, 
Rosseland  diffusion  or  Planck  would  be  expected  to  dominate  in  the  ELGG  barrel 
section.  The  following  is  an  excerpt  of  simulations  accomplished  using  Mathcad  and 
shows  that  the  Rosseland  diffusion  approximation  for  the  temperatures  and  pressures 
considered  is  not  a  valid  assumption.  The  medium,  in  this  case  hydrogen,  turned  out  to 
be  optically  thin,  so  the  Planck  method  should  be  exi  scted  to  be  the  dominant  player  in 
order  to  calculate  the  heat  transfer.  However,  model  development  will  be  conducted 
using  both  methods  for  comparison  purposes  to  the  results  obtained  from  MACH2.  The 
results  of  these  calculations  are  in  appendices  C  and  h.  The  following  is  the  basic 
methodology  used  to  run  the  Mathcad  simulations  presented  in  the  appendix; 
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Figure  3.7  Test  case  for  Rosseland  and  Planck  evaluations 

In  this  figure,  A,  is  the  increment  considered  along  the  x  axis.  The  scale  size  of  the 
system  that  considered  is  10  millimeters,  thus  .lA  would  be  I  mm,  etc...  The  temperature 
gradient,  assumed  to  be  linear,  turned  out  to  be  900K/nim.  The  following  is  a  graph 
calculated  with  Mathcad  which  shows  the  Rosseland  and  Planck  mean  fi'ee  paths  at  1,000 
and  100  atmospheres.  Note  that  the  Rosseland  mean  free  path  is  larger  than  that  of  the 
Planck  mean  fi'ee  path,  suggesting  that  the  radiation  coming  out  of  the  working  fluid  may 
be  flux  limited,  and  that  the  Rosseland  approximation  may  not  be  a  valid  choice  for  the 
C(  nditions  present  in  the  barrel  region  of  the  ELGG. 
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Figure  3.8  Results  of  Mathcad  Rosseland  and  Planck  mean  free  path  versus 
temperature  calculations  (l.OOOK  to  lO.OOOK) 


Figure  3.9  Results  of  Mathcad  Rosseland  and  Planck  mean  free  path  verses 
temperature  calculations  (lO.OOOK  to  35,OOOK) 


63 


3,8  -  MACH2  OPERATION  AND  TEST 


"MACH2  is  a  2  1/2  dimensional  MHD  computer  simulation  tool  for  problems  in 
complex  geometry's"’^  This  code  uses  the  technique  of  domain  decomposition  with  the 
use  of  logically  rectangular  blocks  and  uses  an  arbitrary  Lagrangian-Eulerian 
hydrodynamic  algorithms.  This  program  was  written  by  the  combined  efforts  of  several 
military  and  civilian  scientists  at  the  Phillips  Laboratories  at  Kirtland  AFB,  New  Mexico. 

MACH2  utilizes  an  input  deck,  labeled  im2,  for  problem  identification  and  to  inform 
the  MACH2  code  which  s^t  of  equations  to  use.  For  example,  it  is  possible  to  have  the 
code  numerically  solve  the  flow  problems  within  the  barrel  region  of  the  ELGG  with 
viscous  effects  "turned  ofip’,  and  have  radiation  effects  "turned  on".  This  greatly  helps  in 
the  analysis  of  the  dynamics  of  the  system,  and  allows  for  selective  studies  of  various  flow 
parameters. 

MACH2  affords  a  gi  eat  deal  of  flexibility  in  designing  a  particular  problem,  and  has  a 
wide  variety  of  output  data.  In  order  to  test  the  code  and  verify  it's  proper  operation,  an 
ira2  was  setup  which  had  all  effects  "turned  off",  and  the  code  was  run  to  calculate  the 
exit  conditions  of  the  ELGG.  After  33,000  cycles  using  a  time  step  of  1 .0  e-8  seconds, 
the  code  reached  a  steady  state  solution.  The  total  duration  of  the  flow  through  the  1 
meter  long,  30  mm  inner  diameter  tube  was  found  to  be  2.9891  e-03seconds.  The  exit 
conditions  predicted  by  MACH2  reached  steady  state  operation  and  had  exactly  the  same 
values  as  the  inlet  conditions.  Samples  of  the  im2  and  output  data  run  for  this  test  are 
included  in  appendix  K.  Due  to  limitations  in  the  availability  of  hydrogen  data  within  the 
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MACH2  code  itself,  this  test,  and  subsequent  runs  will  use  argon  as  the  working  fluid. 
Although  argon  has  vastly  different  thermal  and  physical  characteristics  from  hydrogen, 
the  trends  and  observations  made  with  the  hydrogen  models  discussed  in  this  thesis  will 
be  applicable  to  the  results  obtained  with  MACH2. 

The  test  results  of  the  MACH2  coding  show  that  the  code  is  working  correctly  and  is 
providing  data  which  should  be  reliable,  It  was  important  to  run  this  test  for  two  reasons. 
The  first  is  the  need  for  MACH2  to  be  tested  at  various  different  geometry's  and 
conditions,  and  the  second  was  to  foster  familiarity  with  the  code,  the  im2,  and  the  output 
produced. 


CHAPTER  IV 


ANALYSIS  OF  HEAT  TRANSFER 

4.1  ~  DETERMINATION  OF  DEGREE  OF  IONIZATION  AND  WORKING 
FLUID  COMPOSITION 

In  order  to  ensure  the  applicability  of  the  models  developed  in  this  thesis,  it  is  prudent 
to  consider  how  the  working  fluid  is  affected  by  the  conditions  it  will  encounter  within 
the  ELGG  barrel  region.  The  degree  of  ionization  and  the  species  created  by  these 
conditions  are  important  because  these  effects  an  greatly  affect  the  way  in  which  the 
working  fluid  behaves  and  what  order  of  magnitude  can  be  expected  for  heat  change. 

The  ionization  of  the  working  fluid  essentially  robs  some  of  the  flow  energy  within  the 
fluid  and  if  the  degree  of  ionization  is  high,  less  flow  energy  will  be  available  for  projectile 
acceleration.  Additionally,  recombination  effects,  if  significant,  will  change  the  flow 
energy  as  well. 

Ionization  is  a  process  which  removes  an  electron  fi'om  an  atom,  leaving  the  atom  with 
a  net  positive  charge.  The  first  ionization  potential,  V,  (measured  in  eV),  of  the  atom  is 
the  energy  required  to  completely  remove  an  outer  electron  from  it's  normal  state  in  a 
neutral  atom  to  a  distance  beyond  the  sphere  of  influence  of  the  nucleus. 

In  a  plasma,  the  temperature  is  high,  and  the  ionization  is  essentially  due  to  thermal 
ionization.  Ionization  by  thermal  collision  will  eventually  result  in  a  maxwellian  velocity 
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distribution  and  involves  an  energy  transfer  from  the  fast  moving,  energetic  electrons  and 
the  heavy  particles.  The  degree  in  which  the  electron  temperature  decreases  dep<  nds  on 
the  rate  of  recombination,  transfer  of  e  gy  to  vibrational,  rotational,  and  translational 

energy  states,  and  radiation  effects.  It  also  depends  on  the  individual  specie?,  the  volume 
and  the  initial  temperature  of  the  electrons,  ions  and  neutrals.  Additionally,  it  depends  on 
the  physical  boundaries  of  the  system.  By  considering  dissociation,  ionization,  and 
recombination  processes  separately,  it  is  possible  to  get  an  idea  of  what  is  happening  in 
this  system  and  develop  expressions  which  allow  a  first  approximation. 

The  ionizatio,!  level  and  conductivity  are  primarily  determined  by  electron  temperature 
rather  than  the  heavy  particle  temperature.  This  state  is  present  in  the  inlet  rev  ion  of  tiw 
ELGG  due  to  the  use  of  a''c  heating  of  the  working  fluid.  According  to  Jahn,  there 
appears  to  be  some  justification  for  computing  the  ionization  level  from  a  Saha  equarion, 
using  the  electron  temperature.  Additionally,  Rosa  states  that  the  use  of  the  equilibrium 
constant  for  calculation  of  species  is  justified  because  it  is  simply  a  relation  between  the 
forward  and  reverse  rate  constants  and  therefore  holds  for  non-equilibrium  conditions  as 
well. 

Using  the  Saha  equation  discussed  in  section  3.1,  the  following  graphs  and  data  were 
calculated  at  stagnation  temperati-ies  ranging  from  1,000K  to  35,OOOK  and  at  stagnation 
pressures  ranging  from  100  to  1 ,000  atmospheres.  The  range  of  these  temperatures  and 
pressures  are  expected  to  be  present  witiiin  the  barrel  region  of  the  ELGG. 


Figjwe  4. 1  loiii2ation  at  100  aim 


Fiiiiure4  .2  Ionisation  at  1,000  atm 


From  tnese  charts  and  the  data  cuiitained  in  appendix  A,  it  is  clear  vliat  ionkation 
Vrithin  the  band  secuou  of  the:  LLGG  will  noi  be  an  in^rtant  consideration 
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until  the  inlet  temperature  exceeds  approximately  20,000K  at  1,000  atmospheres.  The 
trends  of  these  charts  show  that  as  the  temperature  increases,  the  ionization  increases  as 
well,  but  with  increased  pressure,  the  ionization  decreases.  At  the  extremes  of  tlie 
conditions  which  cart  be  expected  within  the  ELGG  barrel  region,  the  ionization  rate  at 
inlet  stagnation  temperatures  of  35,OOOK  and  pressures  of  1,000  atmospheres  is 
approximately  50  percent  of  the  working  fluid  composition. 

These  results  are  confirmed  by  considering  the  species  number  densities  provided  by 
R.  W.  Patch  and  included  in  appendix  D.  From  this  data,  it  can  be  seen  that  the  number 
density  of  electrons  and  ions  are  9.291  xlO^  at  a  pressure  of  1,000  atmospheres  and 
30,000K.  By  col^parisor^  the  number  density  at  these  conditions  for  H  is  7.020  xl0“. 
From  tliis  empirical  data,  tlte  hydrogen  is  neuly  70  percent  ionized,  and  this  supports  the 
results  obtained  by  the  Saha  equation. 

As  the  working  fluid  travels  through  the  barrel  region,  the  loss  of  heat  energy  due  to 
the  transport  processes  and  viscous  effects  will  ca^rse  chemical  relaxation  towards 
molecular  hydrogen.  To  determine  the  chemical  relaxation  rate,  it  must  be  understood 
that  all  chemical  processes  take  place  by  molecular  collisions  or  radiative  interactions.^ 
For  purposes  of  this  discussion,  radiative  interactions  will  be  neglected  and  collision 
processes  are  assumed  to  dominate.  According  to  Anderson,  nearly  200,000  collisions 
are  required  for  either  dissociation  of  the  molecule,  or  molecular  recombination  to 
occur.”  Obviously  this  takes  a  finite  amount  of  time,  during  which  non  -equilibrium 
effects  are  present.  Non-equilibrium  can  be  defined  as  the  effects  of  finite  reaction 
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p’ocesi'es  which  require  i.  definite  amount  of  time  in  order  to  react  to  a  rapid  change  in 
conditions  in  which  a  gas  is  subjected. 

The  time  required  for  relaxation  to  occur  is  governed  by  l  iite  constants.  When  the 
working  Iluid  is  heated  by  the  arc  prior  to  the  barrel  ixilet^  a  certain  percenttige  of  the 
reaction  will  be  in  the  form  of  cliss'  ciation.  In  general,  molectjlar  dissodaiion  energies 
are  equal  to  or  less  than  the  ionization  energies  of  the  constituent  atoms.  Due  to  the  vety 
higli  energies  and  electron  temperatures  involved,  it  will  be  assumed  that  regardle-is  of 
whether  or  not  a  system  is  in  equilibrium,  the  rate  at  which  a  reaction  proceeds  in  the 
forward  direction  is  equal  to  the  rate  in  the  reverse  direction  at  a  given  temperarure.'* 
This  assumes  that  other  non-equiiibrium  effects  do  not  afiFect  the  rate  of  reactions 
Therefore,  the  forward  and  rev<a  se  reaction  rates  can  be  related  via  tlje  equilibrium 
constant. 

To  determine  whether  oi  not  the  chentical  processes  are  firot^ti,  in  equilibrium,  or 
non- equilibrium  in  th‘R  barrel  section,  the  chemical  attie  of  recombination- dissociation 
nroj'St  be  compared  to  tlie  characteristic  time  the  plasnia  takes  to  traverse  the  length  of  the 
bwTcil,  If  die  recombination  rate  is  less  than  tl;e  characteristic  transit  time  of  the  working 
fluid,  the  flo.y  is  said  lo  be  f'Xizen. 

.!n  addition  to  chemical  relaxation,  ir-tcmal  modes  of  energy  also  seek  out  equilibrium 
conditions.  /MUio.jgli  chesiiical  rt-ictions  require  somewhere  around  200,000  collisions  in 
ordiir  to  rcar  h  a  ricw  equilibriic??  state,  vibrational  energy  is  exc’ianged  and  equilibrium 
att'vin-t.d  m  a^(pioH.i?aately  an  ordei'  of  liiagnitude  less,  typically  around  20,000  collisions. ” 
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Vibrational  energy  is  quantized  into  different  specific  energy  levels,  and  it  is  tlirough 
collisions  tliat  a  change  fi-om  one  energy  level  to  another  is  accomplished.  If  the  system 
consists  only  of  hannonic  oscillators  at  equilibrium,  a  Boltzman  distribution  will  exist, 
assuming  constant  temperature.  During  the  time  the  distribution  of  particles  is  not 
described  by  the  Boltzman  equation,  the  system  is  in  a  non-equilibrium  state. 

When  a  particle  moves  from  one  energy  level  to  another,  it  undergoes  a  transition 
which  requires  collisions  to  occur.  The  transition  probability  is  governed  by  the 
collisional  fi'equency,  Z,  and  the  number  of  particles  in  the  system,  N.  The  following  are 
used  to  determine  the  collisional  frequency. 


Z  =  nncpJ^ 

(4.1) 

(4,2) 

X  —  1 -  (mean  free  path) 

V  2iuf-n 

(4.3) 

CJ  =  (collision  cross  section) 

(4.4) 

The  transfer  of  energy  through  collisions  from  vibrational  tc  translational  energy 
creates  a  temperature  rise  in  the  system.  While  vibrational  energy  is  experiencing  the 
molecular  collisions,  the  working  fluid  is  in  a  non-equilibrium  s^ate  and  energy  is  "locked" 
inside  the  vibrational  energy  levels,  thus  some  of  the  energy  of  the  Quid  is  not  available 


71 


for  projectile  motion.  According  to  Anderson,  the  vibralional  relaxation  time  is  a 
function  of  both  temperature  and  pressure.  This  comes  form  the  fact  that  the  transition 
probability  and  the  collision  frequency  are  also  fiinctions  of  pressure  (:\r  indication  of 
how  probable  a  collision  will  be)  and  temperature  (indicating  the  relative  kinetic  energy  of 
the  particles). 

The  following  is  a  table  of  calculated  relaxation  times  for  botli  ciiemicai  reactions 
(recombination)  and  vibrational  energy. 


Table  4. 1  Calculation  of  relaxation  t.\Ties  a?  variovis  velocities,  tempci  atures 
and  constant  pressure  of  lOOC  atm  (**diita  ffon;  Brown‘d) 


11,605  4.216x10’ 


11,605  3.713x10-" 


coll'ssioi-i  arer 

««  IS 

(.mciers)  Is 


5.157x10"' 


4.641x10-’  5,386x10*  i  2.424x10''  I  9.749x10 


,7.424x.!0'’  1 

1  9.7490110- 

1 

[  - -  - 

1.34lx;0"' 


8  228x10’ 


2.377x10“ 


6  938x10“ 


6  244x10-’ 


1  804x10“ 


l.SWxlU*  I  2.’i?.4xl0'‘  I  9  749x10- 


2.735x10"  l,367x)0'‘  j  S.178)t!0-'‘’ 


rOSi.tlO'"  1.367x10'*  1  g.^/axio'" 

3.603x10"  !.802xf0'*l  7.2‘,  !.<10 


3.603x.’0''  1.802x10'*!  7.271x10''' 

I  ! 


3.603x10"  1.802x10''  7  2/1x10"' 


Coa'-ideriEg  the  results  of  the  calculations  for  both  the  chemical  and  vibrational 
relaxation  timet;,  it  is  clear  that  even  the  slowest  chemical  relaxation  time  occurring  at  an 
inlet  temperature  of  1  eV  and  velocity  of  13,000  m/sec,  the  flow  is  able  to  reach 
^uilibfium  wthin  .033  diameters.  This  suggests  that  non-equilibrium  will  not  be  a  factor 
in  the  barrel  region,  and  that  the  particles  can  be  coiisidered  in  a  Boltzman  distribution 
throughout  the  length  of  the  barrel.  Vibrational  modes  are  also  in  equilibrium. 
Additionally,  f  ozen  flow  losses  are  negligible  for  the  entire  operating  range  of  inlet 
temperalui  es.  As  the  flow  progresses,  the  electron  and  heavy  particle  ternperamres 
should  rapidly  achieve  equilibrium,  and  associated  with  this,  a  rapid  loss  of  ionization  will 
result,  leaving  the  working  fluid  composed  of  only  molecular  hydrogen  and  atomic 
hydrogen  at  the  exit  of  the  barrel.  The  effects  of  recombination  will  serve  to  increase  the 
amount  of  energy  within  the  working  fluid,  tending  to  counter  the  efiects  of  transponi; 
los  ies,  thus  providing  more  energy  for  the  purpose  of  projectile  motion. 

The  following  figure  demonstrates  the  species  number  densities  within  the  ELGG 
barrel  region  aj^d  corroborates  the  findings  of  the  above  calculations.  At  the  inlet 
stag/iatJOi'  ',emperatu.re  and  pressure  of  30,000K  and  1,000  atmosplieres,  figure  4.4 
sugjtiests  that  the  working  fluid  is  mainly  composed  of  ionized  atomic  hydrogen,  with 
molecular  hydrogen  making  up  only  s  fraction  of  the  remainder  of  the  fluid.  Thus  it  can 
be  assumed  that  the  working  fluid  is  completely  dissociated,  and  the  number  of  fi'ee 
electrons  is  of  the  same  order  of  magahude  as  that  of  atomic  hydi  ogen.  Note  that  the 
stagnation  tempeiature  and  pressure  in  the  inlet  region  of  the  barrel  rapidly  decreases  as 
the  clTecs  of  conduction,  raviiation  and  viscous  effects  deplete  t!ie  energy  from  the  fluid. 
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This  causes  a  n^id  recombination  of  atomic  hydrogen,  H,  into  molecular  hydrogen, 
Patch  describes  the  number  density  of  H  decreasuig  from  6.757  xlO“  at  lO.OOOK  to  8.97 
xlO“  at  4,OOOK  while  Hj  increases  from  5.695  xlO“  at  lO.OOOK  to  1.745  xlO”  at 
4,000K.  Temperature  drops  of  this  magnitude  typically  can  be  expected  to  occur  within 
the  first  10  diameters  into  the  barrel  region,  which  equates  to  approximately  3.34  x  10* 
seconds  at  an  inlet  velocity  of  3,000  m/sec.  The  following  graph  from  Patch 
demonstrates  the  number  densities  of  species  at  1,000  atmospheres. 


Figure  4.3  Number  densiiic;-.  of  Hydrogen  at  1,000  atm 


From  this  discussion  and  the  resuils  o.f  the  Sah;'*  equaviori  and  data  provided  by  Patch, 
the  use  of  the  molecular  hydrogen  and  the  va.vious  thermodynamic  and  opacity  p^ropenies 
associated  with  H,  witliin  the  models  developed  is  jjsiiiiable,  v/ith  only  siitall  variances  at 


inlet  temperatures  of  up  to  lO.OOOK  and  pressures  of  1,000  atmospheres.  The  working 
fluid  within  the  barrel  region  of  the  ELGG  for  these  conditions  can  best  be  described  as  a 
weakly  ionized  plasma  at  best,  and  is  actually  a  high  temperature  flow  of  slightly 
dissociated  gas  with  rapid  recombination  within  the  first  few  diameters  into  the  pipe. 
However,  at  inlet  stagnation  temperatures  in  excess  of  lO.OOOK,  the  dissociation  begins 
to  dominate,  and  the  thermodynamic  properties  of  the  working  fluid  begin  to  experience 
rapid  fluctuations.  Fortunately,  by  using  the  empirical  data  from  Patch,  the  change  in 
species  is  included  in  the  thermodynamic  data  used  for  the  heat  transfer  analysis  within 
this  thesis.  This  also  applies  to  the  increasing  ionization  of  the  working  fluid  as  inlet 
temperatures  exceed  approximately  20,000K. 

From  experimental  data  gathered  by  Zube,  the  translational  relaxation  occurs  so 
quickly  that  non-equilibrium  in  this  energy  mode  can  be  observed  to  be  negligible  for  the 
regions  of  interest  in  the  ELGO.  Consider  the  following  table: 


Table  4  2  Translational  Relaxation  times 


Process 

-  -  1 

relaxation  time  (seconds) 

elecuon  self  relaxation 

10  "  to  10"’ 

hydrogen-hydrogen  kinetic  relaxation 

10’ 

jelectron-proton  kinetic  relaxation 

10^ 

From  this  information,  and  the  fluid  velocities  expected  at  the  inlet  of  the  ELGG 
barrel,  the  slowest  translational  relaxation  takes  approximately  .0005  diameters  to  reach 
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equilibrium.  This  demonstrates  that  the  translational  energy  is  in  equilibrium  immediately 
past  the  inlet,  and  coniirms  the  assumption  that  translational  non-equilibrium  can  be 
neglected.  This  holds  true  for  rotational  non-equilibrium  as  well. 

4.2  -  CONDUCTIVE  HEAT  TRANSFER  ANALYSIS 

In  the  development  of  the  conduction  one  dimensional  model  for  the  ELGG  barrel 
region,  the  assumptions  of  a  calorically  perfect,  inviscid  gas  with  steady  flow  have  been 
used.  From  the  discussions  of  the  species  number  densities  and  the  degree  of  ionization, 
these  assumptions  are  acceptable  for  development  of  trends  and  general  conductive  heat 
transfer  calculations.  These  results  will  give  insight  into  the  processes  taking  place  as 
well  as  data  in  which  to  validate  the  MACH2  coding  and  operation  for  problems  of  this 
kind. 

As  discussed,  conduction  occurs  in  two  main  ways.  Conduction  acts  internally,  which 
tends  to  transfer  heat  energy  within  the  fluid  itself,  and  by  contact  between  the  fluid  and 
the  walls  of  the  barrel.  Internal  conduction  along  the  flow  direaion  has  been  assumed  to 
be  negligible  compared  to  the  conduction  to  the  walls . 

To  begin  the  discussion  of  the  effects  of  conduction  only,  neglecting  viscous  effects, 
radiation,  and  convection  for  the  moment,  a  set  of  data  for  conduction  has  been 
computed  beginning  with  inlet  stagnation  temperatures  ranging  from  8,000K  to  34,OOOK 
(.68  eV  to  2.93  eV)  and  pressures  from  85^^  to  1,1 50  atmospheres.  By  using  the 
tabulated  data  on  the  thermal  conductivity  of  hydrogen  throughout  a  wide  range  of 
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temperatures  and  pressures,  the  following  graph  has  been  compiled  which  shows  the  heat 
transfer  per  unit  mass  at  various  temperatures  and  pressures. 


(10000Kand850atm) 
(lOOOOK  and  1150  atm) 
(11605Kand850atin) 
—  (23000Kandll50atm) 


Figure  4.4  Heat  transfer  per  unit  mass  for  conduction  only 


This  graph,  and  the  associated  data  contained  in  appendix  F,  indicate  several  key 
findings  about  the  conduction  within  the  barrel  of  the  ELGG.  The  first  is  the  effect  of 
velocity  on  the  heat  transfer.  For  this  model,  the  effects  of  velocity  changes  is  noticed  for 
a  variety  of  conflicting  reasons.  The  duration  of  time  the  fluid  remains  in  the  barrel 
region  is  an  important  consideration.  Conduction  is  highly  dependent  on  the  physical 
contact  between  the  fluid  and  the  wails,  and  by  changing  the  duration  that  this  occurs  is 
reflected  in  the  heat  transfer.  Velocity  also  affects  the  static  temperatures  and  pressures 
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and  these  effects  are  felt  by  conduction  through  the  amount  of  mass  present  within  the 
flow  and  the  temperature  gradients  created  between  the  walls,  assumed  to  be  at  a 
constant  300K,  and  the  fluid.  Finally,  the  change  in  the  temperature  and  pressure  serves 
to  change  the  thermal  conductivity,  which  correspondingly  changes  the  rate  of  heat 
transfer. 

The  figure  graphi  ly  depicts  these  various  effects,  and  the  following  trends  are 
noted.  First,  if  the  inlet  temperature  is  held  fixed,  and  the  inlet  pressure  (both  are 
stagnation)  is  varied,  the  heat  transfer  actually  decreases  with  increasing  pressures.  This 
can  be  explained  by  considering  the  conduction  heat  flow  equation.  The  heat  flow  is 
given  in  terms  of  Joules/Kg,  which  equates  the  heat  flow  to  the  unit  mass.  As  pressure 
increases,  the  amount  of  mass  in  the  fluid  increases.  This  effect  causes  the  heat  flow  per 
unit  mass  to  decrease,  which  is  what  the  chart  depicts.  What  tliis  means  for  the  £LGO,  is 
that  for  an  increase  in  base  pressure,  there  will  be  a  higher  mass  flow  rate  (1  j  a  point),  and 
with  it,  more  heat  should  be  expected  to  be  conducted  to  the  barrel  walls. 

However,  it  can  be  seen  that  by  holding  the  inlet  pressure  constant,  but  increasing  the 
temperature,  the  conduction  heat  transfer  per  unit  mass  increases.  This  is  explained  by 
the  fact  that  higher  temperature  gradients  are  experienced  by  the  flow,  and 
correspondingly,  higher  conduction  heat  transfer  will  result.  Additionally,  with  increasing 
temperature  at  the  range  of  constant  pressures  to  be  expected  in  the  ELGG,  the  thermal 
conductivity  increases,  also  increasing  the  amount  of  heat  transfer  that  should  be 


observed. 
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Considering  the  range  of  heat  transfer  based  on  inlet  conditions,  it  is  apparent  that  a 
wide  range  of  values  can  be  expected.  The  following  table  is  a  section  of  the  data 
calculated  and  included  in  the  appendices.  It  illustrates  the  heat  transfer  and  percent  of 
inlet  power  remaining  for  various  inlet  conditions. 

Table  4.3  Conduction  data  for  heat  transfer  and  percent  power 
available  at  exit 


inlet  velocity 

inlet 

temperature 

inlet  pressure 

exit  pressure 

beet  transfer 

%  power  at 

exit 

1,500 

10,000 

1,000 

1,003 

-3.956x10’ 

72.09 

13,000 

■ 

10,000 

1,000 

2,260 

-1.588x10’ 

88.8 

1,500 

23,000 

1,000 

1,011 

-3.125x10' 

40.14 

13,000 

23,000 

1,000 

1,435 

-5.501x10’ 

83.13 

1,500 

11,605 

1,000 

1,008 

-5.432x10’ 

67.53 

13,000 

11,605 

1,000 

1,930 

-1.818x10’ 

88.96 

1,500 

10,000 

r' 

850 

857.6 

4.654x10’ 

67,16 

13,000 

10,000 

850 

2,402 

-1.868x10’ 

86.82 

1,500 

10,000 

1,150 

1,158 

-3.44x10’ 

75,73 

13,000 

10,000 

1.150 

2,288 

-1,381x10’ 

90.08 

From  this  table,  it  is  evident  that  conduction  effects  are  most  pronounced  at  low  inlet 


velocities  and  high  inlet  temperatures.  The  power  remaining  at  the  exit  is  highly 
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dependent  on  the  velocity,  as  can  be  seen.  The  highest  velocities  correspond  the  shortest 
duration  of  conduction  within  the  barrel  itself 

4.3  -  ANALYSIS  OF  VISCOUS  EFFECTS 

Viscous  or  frictional  effects  on  a  fluid  flowing  through  a  pipe  has  been  examined  with 
a  one  dimensional  calorically  perfect  model  and  the  results  follow  in  this  section.  As 
mentioned,  viscous  effects  teiid  to  drive  the  inlet  Mach  number  toward  unity,  thus  the 
possibility  of  choked  flow  is  present  whenever  viscosity  is  taken  into  account.  Also  note 
that  the  efiects  of  viscous  action  are  opposite  to  those  of  the  transport  processes  and  this 
difrerence  will  be  explored.  The  results  in  this  section  are  for  viscous  effects  only. 

In  order  to  address  the  effect  of  friction  on  the  flow  of  gases  down  the  barrel  region  of 
the  ELGG,  the  choice  of  friction  factor  must  be  discussed.  There  are  analytical 
expressions  for  the  calculation  of  friction  flictors  for  laminar  flow  in  a  pipe,  but  there  is 
not  a  definitive  expression  for  turbulent  flow.  As  previously  discussed,  fnction  factor  in 
turbulent  flow  is  a  function  of  Reynolds  number,  Mach  number  and  the  ratio  IVD  The 
ratio  L/D  for  the  ELGG  discussed  in  this  thesis  turns  out  to  be  33.3,  wliich  according  to 
Keenan  and  Neumann  is  within  acceptable  limits  for  supersonic  flow.  Additionally,  the 
range  of  Reynolds  numbers  are  between  those  discussed  for  Keenan  and  Neumann,  which 
is  suggestive  that  their  research  and  conclusions  are  applicable  here. 

The  following  are  tables  of  results  based  on  numerical  simulations  ran  using  the 
procedures  detailed  in  the  third  chapter.  These  tables  show  the  impacts  of  various 
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friction  factors,  inlet  stagnation  temperatures  and  pressures  at  a  specific  inlet  velocity  on 
the  viscous  effects.  The  method  of  evaluation,  and  the  full  table  of  calculations  and 
results  is  located  in  the  appendices. 


Table  4  4  Viscous  effects  with  varying  input  data  and  fiiction 
factors  with  constant  velocity  of  4,500  m/sec 


— 

fiiction  factor 

inlet  temp 

inlet  press 

exit  press 

then  length 

press  ratio 

S.OxlO'^ 

10,000 

1,000 

988.83 

37.11 

98.89 

30x10* 

10,000 

850 

840.59 

37.112 

98.89 

3.0x10’ 

10,000 

1,150 

1,137 

37.11 

98.87 

3.0x10’ 

8,000 

1.000 

986.1 

28.43 

98.61 

■»  ftvin  ’ 

1 

11,605 

t  An/\ 

1  ,vn;u 

990.48 

44.14 

yy.U3 

2.0x10’ 

10,000 

1,000 

992.62 

55.67 

99.26 

4.0x10’ 

10,000 

1,000 

985.2 

61.48 

98.52 

Table  4.5  Viscous  effects  with  varying  input  data  and  fiiction 
factors  with  constant  velocity  of  13,000  m/sec 


. . .  ■■ 

fiiction  tactor 

inlet  temp 

inlet  press 

exit  press 

tbeo  length 

pren  ratio 

3.0x10’ 

10,000 

1.27 

28.06 

3.0x10’ 

10,000 

850 

238.67 

1.27 

28.08 

3.0x10’ 

10,000 

1,150 

323.07 

1.27 

28.09 

3.0x10’ 

8,000 

1,000 

96.01 

1.63 

9.6 

3.0x10’ 

11,605 

1,000 

402.54 

1.04 

40.25 

2  0x10’ 

10,000 

1,000 

358  91 

1.9 

35.89 

10,000 

.. 

1,000 

183.78 

1.04 

18.38 
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The  data  used  to  determine  the  viscosity  of  the  working  fluid  was  obtained  from  Patch 
and  from  INSPI  and  is  included  in  appendix  D.  This  data  takes  into  account  the  differing 
species  due  to  the  various  temperatures  and  pressures  inherent  in  the  barrel  region  of  the 
ELGG. 

From  the  data  tables,  and  the  information  included  in  the  appendices,  it  is  clear  that 
the  faction  factor,  as  well  as  the  inlet  conditions  largely  effect  the  flow  characteristics. 
Additionally,  the  effects  of  the  theoretical  length  are  apparent.  First,  note  that  at  low  inlet 
velocities,  the  effect  of  pressure  changes  is  nearly  negligible  on  the  viscous  losses. 
However,  inlet  stagnation  temperature  changes  has  a  more  pronounced  effect  since  the 
temperature  directly  effects  the  Mach  number.  At  the  relatively  low  inlet  velocity  of 
1,500  m/sec,  the  Mach  number  ranges  from  .  18  to  .22.  Stagnation  losses  in  this 
environment  would  be  expected  to  be  minimal,  as  is  borne  out  by  table  4.2. 

In  addition,  the  impact  of  changing  the  friction  factor  can  be  seen  in  this  table  as  well. 
Recall  that  the  effect  of  friction  is  to  drive  the  inlet  mach  number  towards  unity.  As  the 
friaion  factor  increases,  more  frictional  effect  >  are  noticed  by  the  flow,  and  the  higher  the 
stagnation  losses.  Hence  by  going  from  a  friction  factor  of  .003  to  .004,  there  is  a 
corresponding  drop  in  the  stagnation  pressure  at  the  exit.  The  results  of  the  changing 
friction  factor  are  clearly  evident  in  the  preceding  tables. 

The  changing  of  the  initial  conditions  and  friction  factors  are  even  more  dramatic  at 
the  higher  speeds  indicated  in  table  4.3.  This  is  to  be  expected  since  the  shear  stress  at 
the  wall  is  a  function  of  flow  velocity  and  the  Reynolds  number.  Note  again  that  inlet 
pressure  is  not  really  a  factor  due  to  the  relatively  small  density  changes,  but  inlet 
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temperature  dramatically  changes  the  flow  characteristics.  As  temperature  increases,  the 
Mach  number  drops,  the  Reynolds  number  increases  and  the  shear  stress  at  the  wall 
decreases,  thus  explaining  the  increase  in  the  stagnation  pressure  ratio  at  the  exit.  In 
effect,  higher  temperatures  tend  to  decrease  the  frictional  effects,  at  least  in  the  realm  of 
the  ELGG  operation  studied.  The  following  figures  graphically  demonstrate  the  effects 
of  changing  inlet  conditions  at  constant  fiiction  factors,  and  constant  inlet  conditions  at 
different  friction  factors. 


P\0.6 

-Q- 
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VX^ 

(lOOOOKand  lOOOatm) 
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■O'  (SOOOKand  1000  atm) 


Figure  4.5  Ratio  of  stagnation  pressures  verses  velocity  with  a  fiiction 
factor  of  .003 


Figure  4.6  Ratio  of  stagnation  pressures  verses  velocity  at  different 
friction  factors 

As  discussed  in  the  model  development  earlier,  frictional  effects  tend  to  drive  the 
Mach  number  towards  unity,  thus  bringing  the  distinct  possibility  of  choked  flow.  By 
looking  at  the  data  in  the  previous  tables,  and  appendices,  it  is  apparent  that  theoretical 
length  changes  not  only  due  to  changes  in  frictional  factor,  but  also  by  inlet  conditions. 
Again,  figure  4.5  demonstrates  that  pressure  changes  at  the  inlet  have  negligible  effect  on 
the  theoretical  length,  but  temperature  changes  have  a  dramatic  effect.  This  is  due  to  the 
strong  influence  Mach  number  has  on  the  calculation  of  the  theoretical  length.  For  the 
cases  where  temperature  was  at  10,000K  and  a  friction  factor  of  .003,  the  entire  range  of 
inlet  velocities  from  4,500  m/sec  to  12,500  m/sec  are  not  acceptable  due  to  choking  in  the 
pipe.  With  an  inlet  temperature  of  8,00OK,  this  range  dropped  to  include  velocities  as 
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low  as  4,050  m/sec  and  as  the  temperature  rose  to  2.5  eV,  velocities  needed  to  be  higher 
than  20,500  m/sec  in  order  to  preclude  choking.  These  inlet  conditions  must  be  avoided. 

It  is  highly  likely  that  in  the  actual  barrel  of  the  ELGG,  the  specific  heats  will  not 
remain  constant,  and  that  their  ratio  will  change.  This  effect  on  the  theoretical  length  of 
the  tube  has  been  calculated ,  and  the  following  trends  were  discovered.  As  the  ratio  of 
specific  heats  decreases,  the  theoretical  length  tends  to  increase  at  the  same  Mach 
numbers.  Thus,  a  decreasing  y  is  actually  beneficial  in  terms  of  helping  to  prevent 
choking  of  the  flow.  The  opposite  is  true  for  an  increase  in  y.  This  data  is  located  in  the 
appendices,  and  in  section  3.4. 

From  the  results  of  this  analysis,  it  will  be  important  to  compare  the  effects  of 
radiation,  conduaion,  and  convection  within  the  range  of  inlet  conditions  which  prevent 
choked  flow. 

4.4  -  ANALYSIS  OF  CONVECTION 

Convection  is  determined  in  the  ELGG  by  the  pressure  differentials  resulting  in  forced 
flow  within  the  gun  barrel.  As  discussed,  the  flow  entering  the  barrel  forms  boundary 
layers  which  grow  to  meet  at  the  tube  axis  However,  the  large  values  of  Reynolds 
numbers  t^ori  the  order  of  10* )  precludes  boundary  layer  growth  analysis  since  fiilly 
developed  conditions  persist  immediately  past  the  tube  inlet.  According  to  the  work  of 
Siegel,  Piacesi,  and  Bixler,  "this  conclusion  is  justified  and  the  turbulent  boundary  layer 
completely  fills  the  barrel  throughout  its  length."  Additionally,  in  his  studies  of  gun 


barrel  dynamics,  "Williams  assumed  the  flow  to  be  ftilly  developed  turbulent  flow 
imme  lely  because  of  the  large  Reynolds  numbers  attained  by  the  driver  gas  almost 
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immediately." 

With  this  in  mind,  the  thermal  starting  region  was  found  to  begin  immediately  at  the 
entrance  region  of  the  barrel  as  well,  and  both  the  fluid  boundary  layer  and  the  thermal 
boundary  layers  extend  to  the  axis  of  the  tube.  Thus  the  downstream  x'elocity  profile 
does  not  change  from  that  at  the  inlet.  The  following  table  represents  flow  at  inlet 
conditions  of  10,OOOK  and  1,000  atmospheres  at  various  inlet  velocities.  The  full  range 
of  data  is  included  in  the  appendices. 


Table  4.6  Effect  of  fiiction  factor  on  convection  heat  transfer  to  barrel  walls 


inlet 

velocity 

inlet 

temperature 

inlet 

pressure 

exit 

pressuit* 

heat  transfer 

"q"  to  walls 

%  power 

Press  Ratio 

3.0x10’ 

1,500 

10,000 

1,000 

1,004 

■B 

1004x10’ 

3.0x10’ 

4,500 

10,000 

1,000 

1,028 

-1.581x10’ 

88.85 

3.0x10’ 

13,500 

10,000 

1,000 

1,503 

-ZSSxlO" 

94.47 

2.0x10-’ 

1,500 

10,000 

1,000 

1,004 

-1  39x10’ 

84.83 

1004x10’ 

2.0x10’ 

4,500 

10.000 

1.000 

1,028 

-1.581x10’ 

88.85 

1028x10’ 

2.0x10’ 

13,500 

10,000 

1.000 

1,503 

-783.d0‘ 

94.47 

1503x10’ 

4.0x10’ 

1,000 

1,004 

-2.139x10’ 

84.83 

1004x10’ 

4.0x10’ 

4,500 

10,000 

1,000 

1,028 

-1.581x10’ 

88.85 

1028x10’ 

4.0x10’ 

13,500 

10,000 

1,000 

1,503 

-7.83x10’ 

B 

1503x10’ 
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The  first  interesting  thing  to  note  from  the  data  presented  in  the  table  is  that,  other 
than  limiting  the  range  of  inlet  conditions  necessary  to  avoid  choked  flow  as  discussed  for 
viscous  effects,  the  friction  factor  does  not  impact  the  convection  results.  This  is  due  to 
the  fact  that  the  flow  is  immediately  fully  developed  because  of  the  large  Reynolds 
numbers  involved,  and  boundary  layer  effects  are  not  important  in  this  regime  of 
operation  since  the  entire  flow  is  essentially  contained  in  the  boundary  layer. 

What  does  dramatically  impact  the  amount  of  heat  loss  due  to  convection  is  the  inlet 
temperature  and  velocity.  Pressure  again,  as  for  viscous  effects,  does  not  play  a  key  role. 
Consider  the  following  table  containing  calculated  data  from  the  convection  model. 


Table  4.7  Effect  on  convection  heat  transfer  to  walls  using  various  inlet  conditions 


inlet  velocitj 

inlet  temp 

inlet  press 

exit  press 

heal  trans 

%  power 

press  ratio 

1,500 

10,000 

1,000 

1,004 

-2.139x10’ 

84  83 

1004x10’ 

1,500 

10,000 

850 

853.65 

-2.17x10’ 

8469 

1043x10^ 

1,500 

10,000 

1,150 

1,155 

-2.114x10’ 

85.08 

1043x10'' 

1,005 

-1.552x10’ 

86  31 

1005x10’ 

1,000 

1,004 

-2  633x10’ 

84  01 

1004x10’ 

10,000 

1,000 

1,503 

-7.83x10'' 

94.47 

1503x10’ 

As  can  be  seen,  inlet  stagnation  pressure  has  a  minimal  effect  on  changes  in  exit 
stagnation  pressure  and  on  the  amount  of  heat  per  unit  mass  transferred  by  convection. 


87 


Pressure  only  has  an  efl'ect  by  changing  the  density,  which  in  turn  is  the  only  mechanism 
pressure  has  which  can  influence  the  heat  transfer  by  convection. 

Temperature,  however,  has  a  large  effect  on  the  heat  transfer  by  convection. 
Temperature  is  a  main  component  in  the  determination  of  the  heat  transfer,  and  in  fact  the 
heat  transfer  is  directly  proportional  to  temperature  changes.  Velocity  plays  two  key 
roles  in  the  heat  transfer  by  convection.  First,  by  the  fact  that  higher  velocities  simply 
reduce  the  transit  time  of  the  fluid  through  the  tube,  less  time  is  allowed  for  the 
convection  heat  transfer  to  the  barrel  wal’s  to  take  place.  This  limits  the  amount  of  heat 
lost  and  explains  the  larger  percentage  of  inlet  power  remaining  at  the  exit  for  higher 
velocities.  Secondly,  the  velocity  directly  impacts  the  Reynolds  number  and  Stanton 
number.  These  two  factors,  along  with  the  Prandtl  number  dramatically  effects  the  heat 
transfer  coefficient. 

In  01  der  for  the  Reynolds  analogy  assumptions  to  be  valid,  the  Prandtl  number  should 
be  close  to  1 .0,  as  previously  mentioned  in  chapter  3.6.  For  the  entire  range  of  data 
considered,  the  Prandtl  number  ranged  from  no  lower  than  .929  and  to  no  higher  than 
1 .098.  This  represents  approximately  a  9  percent  variation  from  the  optimum  of  1.0,  but 
is  well  within  the  acceptable  limits  set  forth  by  Eckert.  The  following  table  illustrates 
the  results  of  calculations  at  various  inlet  stagnation  temperatures  and  pressures  and  irilet 


velocities. 


88 


Table  4.8  Convection  parameters  at  varying  inlet  conditions 


inlet  velocity 

inlet  temp 

inlet  press 

Reynolds# 

Pruidtlii 

Sunlon  tt 

Coeff  ollrict 

1,500 

10.000 

1.000 

1.273x10* 

1006x10’ 

IxlO’ 

2x10’ 

4.500 

10,000 

1.000 

3.251x10* 

1006x10’ 

.931x10’ 

2x10' 

U.500 

10.000 

1.000 

9.48x10’ 

1006x10’ 

1x10’ 

3x10’ 

1.500 

10,000 

850 

1.082x10* 

1098x10'’ 

1x10’ 

2x10’ 

4.500 

10,000 

850 

2.763x10* 

1098x10’ 

.989x10’ 

2x10’ 

13.500 

10.000 

850 

8.061x10’ 

1098x10’ 

1x10’ 

2x10’ 

1.500 

1,150 

1.464x10* 

929x10’ 

1x10’ 

2x10’ 

4.500 

10.000 

1.150 

3  739x10* 

929x10’ 

.959x10’ 

2x10’ 

13.500 

10.000 

1.150 

1  091x10* 

929x10’ 

1x10’ 

3x10-’ 

1,500 

' 

!.0C0 

!.S!2x!0* 

miiiiQ 

AAIU 

4.050 

8,000 

1.000 

4.173x10* 

1006x10’ 

.896x10’ 

2x10’ 

13,500 

8,000 

1.000 

3.154x10’ 

1006x10’ 

2x10’ 

3x10’ 

1.500 

1 1,605 

1.000 

1.044x10* 

994x10’ 

1x10’ 

2x10’ 

1,000 

2.258x10* 

994x10’ 

976x10’ 

1,000 

1.333x10* 

994x10’ 

2x10’ 

From  this  data  and  the  relationships  of  pressure,  temperature  and  flow  velocity,  we  can 
consider  the  amount  of  heat  transfer  due  to  convection.  The  following  graph  illustrates 
the  heat  flow  per  unit  mass  for  convection  at  various  inlet  velocities.  Note  from  this 
graphical  analysis  that  the  heat  transfer  by  convection  is  only  affected  appreciably  by 
temperature  and  velocity  changes,  and  pressure  changes  are  negligible. 


Figure  4.7  Keai  transfer  per  unit  mass  for  convection  at  a  friction 
factor  of  .003 


4.5  ANALYSIS  OF  RADIATION 


Radiation  is  the  iranspon  of  energy  out  of  the  system  mainly  by  virtue  of  the 
temperature  within  the  fluid.  However,  how  the  thermal  energy  is  transferred  depends  on 
the  properties  of  the  fluid,  most  notably  the  opacity.  To  determine  the  optical  depth  of 
the  working  fluid,  the  models  described  in  section  3.7  are  used.  This  methodology 
initially  indicated  that  at  lO.OOOK  and  lOO  atmospheres,  the  mean  free  path  for  the 
working  fluid  based  on  Rosseiana  ( '.  -‘usion  is  approximately  .45  to  .55  meters  as  opposed 
to  Planck  mean  free  paths  '  f  approximately  .2  to  .3  meters.  This  would  suggest  that  the 
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radiation  coining  from  the  fluid  may  be  flux  limited,  or  in  other  words,  is  limited  by  the 
equilibrium  specific  intensity  or  Planck  fiinction. 

However,  this  data  has  been  shown  to  not  completely  describe  the  system  studied  in 
this  thesis.  The  following  figures  show  the  calculation  of  the  heat  transfer  per  unit  mass 
for  Rosseland  and  Plank  radiation. 


(lOOOOKand  lOOOatm) 
+  (8000K  and  1000  Btm) 
•o-  (ll605Kand  1000  aim) 
(lOOOOK  and  SSOalin) 
(lOOOOKand  llSOatm) 


Figure  4.8  Rosseland  radiation  heat  transfer  per  unit  mass 
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Figure  4. 10  Comparison  between  Planck  and  Rosseland  radiation 
heat  transfer  at  1,000  atm  and  1 1,605K 
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A3  it  can  be  seen  from  these  graphs,  at  various  points  either  method  is  usable,  with 
some  points  flux  limited  by  Planck  radiation  and  others  absorption  limited  by  Rosseland 
radiation.  However,  by  looking  at  conditions  of  23,000K  and  1,000  atmospheres,  it  is 
clearly  seen  that  Rosseland  difilision  is  appropriate  for  the  model  of  radiation  heat  loss 
This  effect  become?  more  pronounced  as  temperature  increases.  At  nearly  3  eV  and 
1,000  atmospheres,  the  Rosseland  difilision  heat  transfer  is  3  orders  of  magnitude  less 
than  Planck,  suggesting  a  highly  absorption  dominated  working  fluid.  As  pressure  drops, 
this  was  also  found  to  be  the  case.  It  seems  that  since  both  Rosseland  and  Planck  have 
similar  mean  free  paths  down  to  approximately  5,000K,  the  methodology  to  use  is 
Rosseland  difilision,  and  not  Planck  as  previously  expiected  .  This  would  tend  to  indicate 
that  the  ELGG  system  studied  is  mainly  absorption  dominated  throughout  its  operating 
range  rather  than  emission  dominated. 

The  difference  between  the  initial  results  and  those  determined  by  the  models 
developed  in  sections  3.5  and  3.7,  is  that  initially  purely  blackbody  radiation,  oT*,  was 
assumed  to  be  sufficient  for  Planck  calculations,  but  this  is  not  completely  true.  As 
pointed  out  in  section  3.7,  there  should  have  been  included  in  the  preliminary  thesis  work 
the  factor  4ap,  which  makes  the  difierence  and  accounts  for  the  Planck  mean  absorption 
coefficient. 

Tlie  variation  in  radiation  heat  transfer  over  velocity  changes  is  explained  by  the  fact 
that  higher  flow  speeds  decreases  the  working  fluid  transition  time  through  the  tube,  thus 
limiting  the  time  for  the  fluid  to  radiate  energy.  This  is  similar  to  the  effects  encountered 


in  the  convection  and  conduction  results  as  well. 
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As  discussed  in  section  3.7,  the  radiation  energy  density  can  be  computed  and 
compared  to  th'^  convection  heat  6ux  and  the  internal  energy  to  determine  when  radiation 
should  be  considered.  The  radiation  energy  density  initially  found  in  the  tube  can  be 
calculated  from  the  inlet  conditions.  It  should  be  noted  that  as  the  temperature  changes 
with  the  flow  down  the  tube,  this  radiation  energy  density  drops  as  the  fourth  power  of 
the  temperature,  thus  radiation  will  rapidly  become  negligible  towards  the  end  of  the 
barrel  region.  The  combined  effects  of  radiation,  convection  and  conduction  serve  io 
force  the  flow  temperature  down. 

The  radiation  energy  density  for  lO.OOOK  and  1,000  atmospheres  inlet  conditions  with 
inlet  velocity  3,500  m/sec  is  found  to  be  5.67x10’  Watts/meter  whereas  the  internal 
energy  is  1.005x10®  Watts/meter,  the  convection  heat  flux  is  approximately  1 .75x10’,  and 
the  conduaion  heat  flux  is  approximately  3.05x10’  Watts/meter.  Clearly  the  effects  of 
thermal  radiation  is  of  the  same  order  of  magnitude  as  that  of  conduction  and  convection, 
and  must  be  included  in  heat  transfer  studies.  I  he  assumptions  used  in  models  used  by 
other  researchers,  namely  Siegel,  Stephenson,  and  the  Weapons  Laboratory  at  Rock 
Island,  that  radiation  is  not  a  factor  is  incorreci  and  by  avoiding  the  inciusioii  of  radiation 
effects  must  raise  questions  as  to  the  accuracy  of  their  results. 

4.6  -  COMPARATIVE  ANALYSIS 

In  this  section,  the  various  transport  mechanisms  and  the  viscous  effects  will  be 
compai  id  and  determination  of  the  regimes  of  operation  where  these  heat  loss  processes 
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dominate  will  be  made.  The  regimes  of  operation  show  where  these  mechanisms  are  of 
equal  order  of  nmgnitudes,  and  trends  discovered  from  the  one  dimensional  models  will 
be  discussed  in  order  to  fully  analyze  the  MACH2  results.  The  power  remaining  in  the 
working  fluid  is  examined  as  is  the  ratio  of  exit  to  inlet  stagnation  pressures.  Internal 
energy  changes  associated  with  the  various  heat  loss  mechanisms  will  also  be  looked  at. 

It  has  been  mentioned  that  the  cooling  effects  of  conduction,  radiation,  and 
convection  on  the  working  fluid  tend  to  drive  the  Mach  number  away  from  unity,  thus 
avoiding  the  prospects  of  choked  flow,  shockwave  development  and  the  corresponding 
loss  in  efficiency  and  the  inevitable  "correction"  of  the  flow  inlet  conditions.  However, 
viscous  effects  act  to  drive  the  Mach  number  towards  unity,  suggesting  that  there  may  be 
a  trade  off  of  effects. 

Since  all  these  processes  are,  in  reality,  closely  coupled  together,  it  is  not  possible  to 
analytically  describe  all  these  effects  simultaneously.  However,  the  data  presented  up  to 
this  point  can  give  some  understanding  of  what  can  be  expected  for  the  flow  of  a  high 
temperature  working  fluid  down  the  barrel  of  the  ELGG. 

Consider  the  ratio  of  stagnation  pressures  from  the  exit  to  the  inlet.  This  is  an 
indication  of  the  overall  effects  each  of  the  dissipative  modes  has  on  the  fluid.  For 
conduction,  it  has  been  determined  that  at  inlet  velocities  below  approximately  Mach  .5, 
the  ratio  of  stagnation  pressures  remains  nearly  constant  for  inlet  pressures  of  1 ,000 
atmospheres,  and  temperatures  up  to  1.8  eV  Above  this  temperature  range,  the  ratio 
increases  slightly,  from  1.003  at  1  eV  to  1.07  at  2.92  eV.  This  shows,  in  effect,  that  the 
temperature  at  constant  pressure  of  1,000  atmospheres  has  negligible  effect  on 
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conduction.  This  was  also  determined  to  hold  true  for  pressure  ranges  from  850 
atmospheres  to  !,150  atmospheres. 

As  the  velocity  was  increased  beyond  Mach  .5,  the  results  changed,  si>owing  a  higl?. 
degree  of  variation  in  the  stagt^atior.  pressui  e  ratio  with  boih  temperatiue  and  ptessaic. 
The  change  from  low  to  high  velocity  ranges  is  manifested  in  the  duration  of  time  the 
fluid  remains  in  the  ELGG  barrel  region,  and  the  variations  noted  indicate  the  dependence 
of  the  thermal  conducti’/ity  to  both  temperature  and  pressure.  According  to  data 
provided  by  iNSPI,  both  temperature  increases  and  pressure  decreases  show  marked 
increases  in  the  thermal  conductivity  of  the  hydrogen  working  fluid. 

For  radiation,  it  has  been  computed  that  the  velocity  changes  do  not  have  the  large 
variai  ions  in  stagnation  pressure  ratios  exhibited  by  conduction,  and  that  changing 
pressure  has  little  efiect  on  the  ratio  as  well.  There  was  nearly  a  linear  correspondence  to 
increasing  temperature  and  increasing  pressure  ratio.  It  is  interesting  to  note  that  both  the 
Rosseland  and  Planck  radiation  data  indicated  pressure  ratios  which  were  to  within  1%  of 
each  other,  although  the  heat  transfer  was  far  greater  for  Planck. 

Viscous  effects  tended  to  decrease  the  exit  stagnation  pressures,  whereas  all  the  other 
heat  loss  mechanisms  tended  to  increase  the  exit  stagnation  pressure.  Tliis  is  consistent 
with  the  particular  models  used,  and  dramatically  shows  the  difference  between  heat 
transport  mechanisms  and  viscous  eflfects.  Viscous  effects  showed  little  change  based  on 
temperature  or  pressure  variations,  but  were  dramatically  affected  by  velocity  changes. 
This  is  consistent  with  expectations  since  viscous  effects  are  driven  by  gradients  in 
velocity  throughout  the  flow  region,  and  higher  velocities  tend  to  create  large  velocity 
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gradients,  most  notably  near  the  tube  v/alls.  Pressure  ratios  as  low  as  .0926  were 
calculated  at  speeds  in  excess  of  Mach  3.2  at  a  temperature  of  8,000K  and  pressure  of 
1,000  atmospheres. 

Com  ection  showetl  suxiilar  results  as  conduction,  bat  the  pressure  ratios  were  never  as 
large  as  those  found  iu  conduction.  Converticn  pn'ssure  ratios  were  negligibly  aifected 
by  pressure  variations,  but  were  strongly  aftected  by  inlet  temperature  increases  and  flow 
velocity.  Pressure  ratios  ranged  from  1.806  at  8,OOOK,  1,000  atmospheres  and  flow 
velocity  of  13,000  m/sec,  to  1 .3 10  at  34,000K,  1,000  atmospheres  and  13,000  m/sec. 

Although  is  cannot  simply  be  a  matter  of  looking  at  these  results  and  combining  the 
information  to  predict  the  final  outcome,  it  is  possible  to  note  that  increases  in  stagnation 
pressure  ratios  arc  associated  witlr  increased  heat  dissipation  from  the  working  fluid,  and 
mat  viscous  effects  tend  to  decrease  this  ratio  As  a  first  approximation,  it  is  noted  that 
as  the  usiet  temperature  of  the  working  increases,  regardless  of  the  pressure  at  the  inlet 
(within  the  scope  of  this  thesis)  the  combined  increase  of  the  pressure  ratio  due  to 
conduction,  radiation,  and  convection  are  practically  offset  by  the  decrease  caused  by 
viscous  dissipation.  However,  as  the  inlet  velocity  increases  to  well  over  Mach  1 .0,  the 
situation  changes,  and  it  can  be  expected  that  \iscous  effects  will  tend  do  dominate  over 
the  other  three.  Again,  this  is  only  a  first  approximation,  and  recall  that  the  effects  of 
radiation,  convection  and  conduction  were  calculated  using  one  dimensional,  steady  state, 
inviscid  models.  The  calculations  from  the  MACH2  code  will  be  able  to  simultaneously 
solve  these  together,  and  a  better  representation  of  the  actual  system  should  be  obtained. 


In  addition  to  the  stagnation  pressure  ratio,  the  heat  transfer  for  convection, 
conduction  and  radiation  can  be  compared,  and  regimes  of  influence  determined. 
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Consider  conduction  heat  transfer  data  obtained  from  calculations  made  using  the  one 
dimensional  model  developed  in  chapter  3.5.  This  data  demonstrates  that  as  the  velocity 
increases  at  any  temperature  and  pressure  combination,  the  heat  transfer  decreases,  due 
mainly  to  tl  ie  redaced  transit  time  of  the  working  fluid.  Generally,  the  trends  indicate  that 
at  a  constant  pressure  of  1,000  atmosphere.s,  and  velocity  ranges  from  500  m/'sec  to 
13,000  m/sec,  the  following  range  of  heat  transfer  is  obtained: 


Table  4.9  Conductive  heat  transfer  data 


Ifnlet  tcinpci^turc  (K) 

OriiCSt  SiidJViiuiiCM 

8,000 

7.464x10’ to  1.816x10’ 

— 

10,000 

1.17x10"  to7.83xl0‘ 

23,000 

3.125x10*  to5.50I.xl0' 

34,000 

1  444x10*10  1.066x10* 

Under  the  same  conditions  mentioned  above,  the  and  Planck,  heal  transfer 

ranges  were  determined  and  are  shown  in  table  4. 10.  Table  4. 1 1  shows  the  data  for 
convection  using  the  same  conditions  as  conduction  and  radiation. 
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Table  4. 10  Radiative  heat  transfer  data 


Inlet  temperature  (K) 

.  ■  "1 

Rosseland  heat  transfer  (J/kg) 

Planck  heat  uansfer  (J/kg) 

8,000 

4.421.X10’ to  1.786.X10’ 

1.034x10“  to  4.717x10’ 

10,000 

1.403.xl0*to  I.252.X10" 

1.303x10“  to  1.219x10'' 

23,000 

2.120x10To  1.412.\10‘ 

8.244x10'°  to  5.521x10'' 

34,000 

1.605x10*10  1.287x10'' 

2.508x10"  to  2.017x10’ 

Table  4.11  Conveaion  heat  transfer  data 


Inlet  torapersture  (K) 

Range  of  heat  transfer  calculated  ( l/kg) 

8,000 

1.522x10'  to  4.007x10'' 

10,000 

2.139x10’  to  7.83x10° 

23,000 

8.924x10’ to  4.428X10'’ 

1  34,000 

1.419x10’ to  7.748x10° 

From  the  above  three  tables,  it  is  now  possible  to  discuss  regimes  of  influence  for  these 
heat  transport  processes.  To  do  this,  the  reg  mes  for  which  data  was  obtained  will  be 
outlined  using  constant  inlet  pressure  of  1,000  atmospheres. 

For  temperature  regions  below  10,000K,  conduction  dominates  the  flow  at  higher 


inlet  velocities,  but  near  lO.OOOK,  both  conduction  and  radiation  arc  noted  to  be  nearly 
equal  in  magnitude.  Convection  tends  to  maintain  at  a  level  of  heat  transfer  which  is  an 
order  of  magnitude  less  than  conduction.  It  is  also  apparent  that  convection  is  more 
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important  in  this  regime  than  either  Rosseland  or  Planck  radiation  when  velocities  in 
excess  of  Mach  .8  are  concerned.  This  region  is  characterized  by  inlet  temperatures  less 
than  1  eV,  heavy  conduction  throughout  the  barrel  region,  and  radiation  effects  rapidly 
diminished  as  the  flow  progresses  throughout  the  barrel. 

For  flow  in  regimes  of  operation  between  1  eV  and  2  eV,  the  flow  has  nearly  equal 
magnitudes  of  conduction,  convection  and  Rosseland  radiation,  although  the  Planck 
radiation  model  shows  several  orders  of  magnitude  higher  heat  transfer.  This  was 
discussed  in  section  4.7,  and  from  this  point,  the  working  fluid  is  absorption  dominated, 
thus  limiting  the  amount  of  radiation  seen  at  the  walls.  Thus  Rosseland  diffusion 
approximations  should  more  adequately  model  the  fluid  in  this  regime  of  operation,  and 


the  data  fcr  Rcsseland  rauiauGu  will  be  used  exclusively. 


The  enects  or  velocity  in  ims 


regime  are  not  as  apparent  as  in  the  below  1  eV  regime,  and  due  to  the  combined  effects 
of  all  three  heat  tiansfer  processes,  the  radiation  will  decrease  due  to  the  temperature 
drop  associated  with  the  heat  loss.  This  reduction  in  radiation  due  to  the  temperature 
loss  will  rapidly  shift  the  r  adiation  component  of  heat  transfer  to  the  regime  discussed 
previously  due  to  the  dependence  orT*.  Additionally,  the  effect  of  heat  loss  is  to  increase 
the  flow  velocity,  thus  further  reducing  the  radiation  expected.  The  corresponding 
increase  in  velocity  as  the  fluid  progresses  will  keep  the  conduction  and  convection  heat 
loss  within  the  same  order  of  magnitude. 

The  final  regime  studied  is  that  from  2  eV  to  3  eV.  In  this  region,  Rosseland  radiation 
dominates  at  inlet  velocities  below  Mach  .7,  and  conduction  is  seen  to  be  a  strong  heat 
transfer  mechanism.  Both  conduction  and  radiation  are  an  order  of  magnitude  greater 
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than  convection  at  low  velocities,  and  at  high  inlet  velocities,  radiation  and  convection  are 
of  the  same  order  of  magnitude.  Again,  the  temperature  decreases  associated  with  heat 
loss  will  tend  to  rapidly  diminish  the  radiative  heat  transfer,  but  radiation  should  remain  a 
large  contributing  factor  throughout  a  significant  portion  of  the  barrel  length.  The 
following  three  figures  graphically  illustrate  the  flow  regimes  discussed: 


Figure  4. 1 1  Representation  ot  transport  processes  in  the  operating  regime 
of  less  than  1  eV  (1,000  atm) 


Figure  4  .12  Representation  of  transport  processes  in  the  operating  regime 
from  1  e V  to  2  eV  (1 ,000  atm) 
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B-CnntectlM 


Figure  4.13  Representation  of  transport  processes  in  the  operating  regime 
from  2  eV  to  3  eV  (1,000  atm) 

From  this  analysis,  it  is  apparent  that  conduction  will  play  a  significant  role  in  the 
flow  of  the  working  fluid  in  the  barrel  region  of  the  ELGG.  Convection  also  plays  a 
major  role,  but  is  only  significant  at  hign  inlet  velocities.  Radiation  has  been  shown  to  be 
absorption  dominated  throughout  the  operating  region,  dominating  at  temperatures  in 
excess  of  2  eV,  and  rapidly  diminishing  as  distance  through  the  barrel  increases. 

To  complete  the  flow  regime  analysis,  the  problem  of  choked  flow  caused  by  viscous 
eflfects  in  addition  to  these  heat  transfer  processes  must  be  included.  As  noted  previously, 
the  viscous  effects  will  limit  the  inlet  flow  velocities  based  on  the  formation  of  a  normal 
shock  in  the  tube.  At  this  point,  the  flow  has  reached  the  maximum  mass  flow  rate,  and 
no  further  increase  can  be  tolerated.  From  figure  3.3  and  table  J.  1 ,  it  was  determined  that 
inlet  velocities  must  be  below  Mach  .63  or  above  Mach  2.36  with  a  fhction  factor  of  .003 
in  order  to  prevent  choking  of  the  flow.  Using  the  equation  of  continuity  (equation  3.19), 
the  maximum  mass  flow  rate  is  determined  to  be  6  596  kg/sec,  with  an  inlet  cross 


sectional  area  of  7.069x10'*  interior  surface  area  of  .094  nt"  and  interior  barrel  volume 
of  7.069x10'*  m^.  The  ELGG  Barrel  length  is  1.0  meters. 

This  changes  slightly  the  flow  regimes  discussed  in  that  much  of  the  flow  detailed 
cannot  actually  be  expected  to  be  seen  in  the  ELGG  barrel.  Since  the  viscous  effects  limit 
inlet  velocities  to  below  Mach  .63,  the  low  speed  discussion  remains  valid.  However, 
since  only  inlet  speeds  in  excess  of  Mach  2.36  are  allowed,  it  can  be  expected  that 
convection  will  be  a  significant  heat  transfer  process  if  supersonic  flow  within  the  barrel  is 
required.  This  will  also  tend  to  decrease  the  radiation  heal  loss  and  conduction  due  to  the 
decreased  working  fluid  transit  time  through  the  barrel. 

The  power  remaining  at  the  exit  of  the  ELGG  barrel  is  an  indicator  of  the  ability  for 
the  working  fluid  to  accelerate  a  projectile.  From  the  results  discussed  in  sections  4. 1 
through  4.5,  it  was  determined  that  the  conditions  favoring  the  largest  percentage  of 
power  remaining  were  high  temperature  and  high  velocities.  This  can  be  seen  in  tables 
4.3,  4.6,  and  G.  1.  This  works  out  well  since  the  Mach  numbers  will  be  much  higher  than 
the  choking  Mach  number  of  2.36.  At  inlet  speeds  of  13,000  m/sec,  the  percentage  of 
power  rcmainiiig  at  the  exit  for  eacli  of  the  transport  processes  can  be  compared. 
Conduction  shows  a  dependence  on  pressure,  ranging  trom  a  low  of  86.82%  at  10,000K 
and  850  atmospheres,  to  90.08%  at  10,000K  and  l,i50  atmospheres.  At  constant 
pressure  of  1,000  atmospheres,  the  range  is  from  83  .98%  at  8,000K  to  88.96%  at 
1 1,605K,  and  then  dropping  off  to  77.88%  at  34,OOOK.  Rosseland  radiation  shows  no 
pressure  dependence,  and  the  percentage  ranges  from  97.3%  to  99.7%,  hence  the  high 
speed  tending  to  limit  the  amount  of  radiation  losses  due  to  low  transit  limes.  Finally, 
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convection  does  ;:ot  exhibit  pressure  dependencies,  but  shows  a  range  from  83.92%  at 
34,000K  to  a  high  of  96,47%  at  8,000K.  Thus  convection  is  again  shown  to  be  more 
effective  at  higher  temperatures  and  velocities.  These  results  also  show  that  conduction 
remains  a  major  contributor  throughout  the  range  of  operating  conditions  of  the  ELGG. 

As  the  flow  progresses,  the  decreasing  temperature  can  act  to  change  the  y  of  the  fluid 
which  will  also  affect  the  limiting  Mach  numbers.  According  to  figure  D.4,  however, 
generally  the  decreasing  temperature  at  a  nearly  constant  pressure  of  1,000  atmospheres 
will  tend  to  decrease  the  y  which  is  actually  beneficial,  since  a  lower  y  has  been  shown  to 
decrease  the  limiting  Mach  number. 

A  final  point  for  discussion,  and  an  important  finding,  is  the  effect  recombination  has 
on  the  overaii  energy  of  the  v/orking  fluid.  When  the  working  fluid  is  heated  by  the 
electric  arc,  a  substantial  portion  of  the  fluid  may  be  ionized.  From  figure  4.2,  it  was 
shown  that  this  could  be  as  high  as  50  percent  at  3  eV.  In  addition  to  the  amount  the 
working  fluid  is  ionized,  it  was  shown  that  it  is  nearly  completely  dissociated  into  atomic 
hydrogen  at  temperatures  as  low  as  1  eV.  Considering  the  dramatic  decrease  in 
stagnation  temperature  at  the  exit,  it  is  necessary  to  consider  the  energy  liberated  by  the 
recombination  of  positively  charged  atomic  hydrogen  and  the  free  electrons.  It  was 
shown  that  the  working  fluid  is  not  frozen  as  it  progresses  down  the  barrel,  and  in  fact  the 
relaxation  rates  are  sufficient  for  equilibrium  to  be  reached  at  every  point  during  the  time 
the  fluid  travels  through  the  barrel.  Because  of  this,  and  the  fact  that  the  fluid  remains 
substantially  atomic  hydrogen  down  to  approximately  lO.OOOK,  the  following 
recombination  reaction  can  be  assumed  to  be  dominant. 
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H^  +  e  — energy  (4.5) 

The  energy  of  ionization  is  given  as  1 .3 1x10^  KJ/mol  for  this  process.  Using  the 
assumption  that  the  recombination  brings  the  electron  to  the  ground  state,  and  data  on 
electron  number  densities  provided  by  Patch”,  the  table  4. 12  was  produced  which 
demonstrates  the  liberated  energy  which  may  be  expected  by  this  reaction.  Calculations 
have  shown  that  stagnation  temperatures  at  low  inlet  velocities  may  drop  to  or  below 
5,200K.  Higher  velocities  and  temperatures  may  not  have  such  a  large  change,  but  will 
retain  orders  of  magnitudes  similar  to  those  shown  in  figure  4. 12.  Below  5,200K,  the 
electrons  are  nearly  100%  recombined  and  molecular  hydrogen  begins  to  dominate  the 
working  fluid 


Table  4, 12  Energy  liberated  during  recombination  (1,000  atm) 


Initial  velocity 

Inlet  stagnation 

Enlct  sUgnAtiOfl 

Exit  static 

number  density 

Energy  released 

temperature 

pressure 

fenqxiature 

clunge  (e  ) 

by  recomb  (J) 

1,500 

10,000 

1,000 

7,169 

5.37X10" 

8.257X10' 

8,500 

10,000 

1,000 

6.1 19 

1.29X10" 

I  OO/iVini 

4  ,  V  ■  V/ 

13,000 

10,000 

1,000 

2,404 

4  55X10"’ 

6.997X10'' 

1,500 

23,000 

1,000 

954 

5.437X10" 

8.361X10* 

8,500 

23,000 

1,000 

17,330 

1.968X10" 

3.026X10* 

13,000 

23,000 

1,000 

10,506 

4.713X10" 

7.247X10’ 

8,500 

34,000 

1,000 

22,920 

5.971X10“ 

8.913X10* 

13,000 

34,000 

1.000 

13,470 

3.756X10" 

5.776X10* 
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This  information  is  very  important  when  compared  to  the  magnitudes  of  energy  loss 
for  each  of  the  transport  processes.  In  order  to  facilitate  direct  comparisons,  the 
following  table  of  energy  loss  for  each  of  the  transport  processes  at  various  initial 
conditions  was  produced. 


Table  4. 13  Transport  energy  loss  at  1 ,000  atm  (*  4500  m/sec  for  convection) 


Initial  velocity 

(m/sec) 

Inlet  stagnation 

temperature 

(K) 

Mass  flow 

(kg/sec) 

Rosscland 

energy  loss(J) 

Planck  energy' 

loss  t  J) 

Conducuon 

/aiergy  loss  (J) 

Convection 

energy  loss 

(J) 

1,500 

10.000 

2.56 

6.856x10' 

7319x10’ 

6.74x10’ 

3.65x10’ 

1,500 

23,000 

1.12 

1.587.x  10’ 

6.171x10’ 

2.339x10’ 

6.679x10’ 

1,500 

34.000 

0/54 

.... 

6.780x10’ 

i.oeixio* 

£  '*ri£..%r*i 

5500* 

10.000 

7.26 

1.410x10’ 

5.665x10’ 

1.649x10’ 

1.837x10’ 

8,500 

23,000 

4.83 

1.770x10’ 

6,914x10* 

3.682x10* 

3.01x10’ 

8,500 

34,000 

2.72 

8.041x10’ 

1.257x10’ 

4.614x10* 

2.89x10’ 

Note  that  the  energy  liberated  by  recombination  plays  little  part  in  the  flow 
characteristics  and  heat  transfer  at  low  initial  temperatures.  However,  as  the  static 
temperature  exceeds  approximately  1 5,000K,  the  energy  from  recombination  approaches 
that  of  the  energy  lost  through  the  transport  processes.  This  is  especially  apparent  at 
higher  velocities  since  the  transport  processes  are  limited  by  the  duration  of  time  the  fluid 
remains  in  the  barrel.  Note  further  that  between  2  eV  and  3  eV,  the  energy  liberated  is  of 


equal  or  greater  magnitude  than  that  from  all  the  transport  processes  combined.  From 
previous  discussions,  Planck  radiation  is  not  considered. 
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What  this  shows  is  that  as  the  temperature  mcreases  into  the  2  eV  to  3  eV  range, 
energy  is  made  available  for  projectile  motion  by  recombination  of  the  free  electrons  with 
ionized  atomic  hydrogen.  This  is  a  significant  finding,  and  provides  strong  support  for 
the  use  of  high  temperature  plasmas  as  working  fluids  for  the  purpose  of  projectile 
motion.  Table  4. 14  illustrates  the  hydrogen  specific  internal  energy  as  well  as  the  degree 
of  ionization  and  dissociation.  From  this  data,  it  is  obvious  that  for  hydrogen,  ionization 
energy  makes  up  a  significant  portion  of  the  internal  energy  at  2  to  3  eV,  and  contributes 
to  the  working  fluid's  ability  to  provide  for  projectile  motion.  This  effea  is  negligible  for 
lovi-'er  inlet  stagnation  temperatures,  at  1,000  aunospheres.  Dissociation  energy  is  made 
available  to  the  working  fluid  only  during  the  recombination  of  atomic  hydrogen  into 
molecular  hydrogen,  but  only  at  temperatures  below  lO.OOOK.  At  higher  inlet 
temperatures,  the  recombination  was  found  to  be  negligible,  and  hence  added  little  to  the 
working  fluid  internal  energy. 


Table  4.14  Hydrogen  internal  energy  data  (8,500  m/sec,  1,000  atm) 


■ 

■ 

inlet  static 

«ntp(K) 

;xit  static 

;cn^(K) 

inlet  degree  of 

lonizalioo 

exit  degree  of 

icnizalion 

inlet  degree 

^fdisiOC 

;xit  liegrce 

ofdiaoc 

4.41x10’ 

2.926x10’ 

7.451 

6,119 

0.71 

■n 

wmm 

23.000 

5.43x10’ 

1078x10* 

20,450 

17,330 

.135 

.0872 

full) 

34,000 

4,63x10’ 

9.54x10’ 

31,450 

22.920 

.5263 

.172 

■ 

full) 

CHAPTERS 


ANALYSIS  OF  MACH2  COMPUTATIONS 


Several  full  model  computations  were  conducted  with  the  MACH2  code,  first 
isolating  each  process  for  verification  of  proper  operating  conditions,  and  then  running 
full  model  simulations  with  all  effects  on  simultaneously.  It  should  be  noted  that  these 
computations  were  performed  with  argon  due  to  the  limitations  faced  Ln  getting  hydrogen 
properties  properly  coded  into  the  program.  Nevertheless,  the  previous  ground  work  and 
results  are  applicable  for  trend  analysis  and  can  provide  methods  for  verifying  the  results 
obtained  fi-om  MACH2. 

The  computations  completed  used  a  calculation  mesh  of  10  icells  (x  direction)  and  10 
jcells  (y  direction)  with  iniet  conditions  of  stagnation  temperature  ranging  from  lO.OOOK 
to  24,000K,  static  densities  from  1 .858  kg/m^  to  43. 163  kg/m^  and  velocities  of  Mach  .5 
and  3.5.  These  inlet  Mach  numbers  were  chosen  based  on  the  theoretical  length  analysis 
in  chapter  3  to  pr  .-elude  chcfked  flow  conditions  within  the  barrel  region.  The 
computations  were  conducted  with  both  Planck  and  Rosseland  radiation  parameters  and 
included  the  effects  of  viscosity,  conduction,  and  chemical  reactions. 
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MACH2  generates  data  in  two  main  ways,  numerical  data  and  plots,  both  of  which 
were  used  for  this  analysis.  A  sample  set  of  data  is  included  in  appendix  L.  Two  cases 
are  presented,  both  using  Planck  radiation  with  inlet  Mach  numbers  of  3.5,  one  at 
lO.OOOK  and  the  other  at  23,OOOK.  Table  5. 1  presents  a  small  portion  of  the  overall  data 
generated  for  five  separate  cases. 


Table  5. 1  Full  model  results  for  MACH2  computations  (average) 


Rad  Model 

Inlet  stag  temp 

(K) 

Inlet  Mach  # 

Stag  Enthalpy 

change 

Ratio  of  stag 

pressures 

U  diameters  for 

developed  flow 

Exit  Mach  # 

Planck 

10,000 

5.0x10' 

-1.818x10* 

1  73 

5.04x10' 

3.38x10' 

Rosseland 

23,000 

3.5 

+3681x10’ 

0.57 

5.10x10' 

2.60x10' 

Planck 

23.000 

3.5 

+1.652x10’ 

0.78 

5.10x10' 

3.09x10' 

Planck 

10,000 

3.5 

-4.772x10* 

1.01 

3  67x10' 

2.92x10' 

Rosseland 

10,000 

3.5 

•4.435x10* 

0.83 

3.67x10' 

2.54x10' 

For  computations  at  inlet  stagnation  temperatures  of  lO.OOOK,  the  flow  develops  a  low 
temperature  region  near  the  wall  with  the  lowest  velocity  in  this  region  as  well.  This 
increases  gradually  for  both  quantities  towards  the  axis,  leaving  the  density  at  the  highest 
values  near  the  wall.  The  flow  field  demonstrates  small  scale  velocity  magnitudes 
directed  slightly  towards  the  wall  with  a  fan  wave  expansion  at  the  exit.  The  flow  shows 
that  boundary  layers  grow  rapidly  to  folly  developed  flow  within  .504  diameters  of  the 
barrel  inlet.  Figure  5. 1  illustrates  this  data.  Note  that  the  proportion  indicated  is 
misleading.  The  barrel  is  1.0  meters  long  and  has  radius  .015  meters  from  the  center  line. 
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Figure  5.1  Full  model  computations 


It  is  interesting  to  note  that  density  is  largest  near  the  wall  and  barrel  exit,  suggesting 
that  the  elective  tube  radius  is  reduced  in  the  flow  direction  due  to  density  variations. 
This  is  supported  by  velocity  data  showing  highest  velocity  near  the  axis  and  trailing  off 
gradually  to  within  2  cells  of  the  wall,  then  a  rapid  velocity  decrease  occurs.  According 
to  data  from  the  Mach2  runs,  axis  velocities  of  5.044xl0^in/sec  decreasing  to  4.263x10^ 
m/sec  within  2  cells  of  the  wall  are  present,  with  a  sudden  drop  to  2.006x10^  m/sec  at  the 


wall  interface. 
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Due  to  the  low  inlet  flow  speed,  sound  wave  transit  times  are  three  times  greater,  on 
average,  than  the  working  fluid  transit  time  down  the  length  of  the  barrel.  This  explains 
the  near  uniform  pressure  contours  throughout  the  barrel  region,  and  suggests  that  the 
flow  has  plenty  of  time  to  communicate  changing  conditions  throu'jhout  the  flow.  This 
holds  true  for  the  barrel  length  as  well  as  from  the  tube  walls  to  the  tube  axis. 

Viscous  dissipation  was  found  to  not  be  a  strong  factor  in  this  flow  due  to  the  low 
inlet  velocities,  as  well  as  the  low  inlet  temperature.  Chemical  recombination  of  the 
working  fluid  appears  to  play  no  significant  role  in  this  operation  regime  as  well.  Finally, 
as  predicted  by  the  trends  discussed  in  chapter  4,  there  is  a  stagnation  pressure  increase  at 
the  exit,  and  the  Mach  number  drops.  This  indicates  a  direct  correspondence  to  the  one 
dimensional  model  developed  in  chapter  3.  (see  table  5. 1  for  data) 

As  the  Mach  number  at  the  inlet  increases  to  3.5  with  the  same  inlet  stagnation 
temperature  of  10,000K,  the  situation  changes  markedly.  There  is  a  high  temperature 
region  created  near  the  barrel  walls  over  the  entire  length.  In  this  region,  velocity  is  at  the 
lowest  value  and  density  is  also  at  the  lowest  value.  This  is  in  direct  contrast  to  the 
subsonic  results  previously  discussed.  A  region  of  high  density  (6.8  kg/m^  verses  1.3 
kg/m’)  persists  near  the  tube  axis.  Vorticily  is  initially  widespread  throughout  the  interior 
of  the  barrel,  but  by  10,895  cycles,  it  is  no  longer  a  factor,  suggesting  that  the  flow  has 
become  fully  developed.  This  is  clearly  seen  in  the  plots  located  in  appendix  L. 

Considering  the  two  radiation  models  employed  in  the  calculations,  Planck  radiation 
results  in  a  higher  enthalpy  loss  thtin  the  Rosseland  model,  suggesting  that  the  flow  may 
be  absorption  dominated,  much  the  same  way  that  the  hydrogen  working  fluid  was  found 
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to  be.  The  change  in  enthalpy  is  much  higher  for  both  models  at  this  higher  inlet  velocity. 

The  flow  shows  an  average  exit  Mach  number  of  2.54  for  the  Rosseland  case  and  2.92 
for  the  Planck  case.  This  demonstrates  that  the  choice  of  Mach  3.5  was  an  effective 
measure  in  avoiding  choked  flow.  However,  the  Mach  number  decreases  are  not  as 
severe  as  predicted,  mainly  due  to  viscous  dissipation  in  the  working  fluid  near  the  walls. 

Since  the  fluid  is  supersonic  throughout  the  barrel,  the  pressure  waves  cannot  keep 
up  with  the  flow,  and  non-uniform  pressure  variations  were  demonstrated.  This  effect 
was  present  in  the  axial  direction,  but  was  not  severe  in  the  radial  direction  due  to  the 
small  tube  radius  and  relatively  short  wave  transition  time  in  the  radial  direction. 

At  inlet  conditions  of  Mach  3.5  and  stagnation  temperature  of  23,000K,  the  effects 
noted  for  the  supersonic  flow  at  the  lower  temperature  are  present,  but  show  much  higher 
magnitudes.  The  temperature  band  near  the  wall  has  higher  temperatures  and  is  thinner 
than  that  demonstrated  at  10,000K.  Boundary  layer  growth  is  very  rapid  (less  than 
8.95x10"*  sec)  in  the  flow  direction. 

The  supersonic  flows  do  not  show  the  density  region  near  the  exit  acting  to  reduce  the 
tube  radius  as  demonstrated  for  subsonic  flow.  In  all  cases  the  temperature,  density,  and 
pressure  approach  towards  steady  state  values  as  time  progresses,  and  rapid  turbulent 
boundary  layer  development  resulted  in  fully  developed  flow  within  one  diameter  into  the 
barrel  region.  Once  this  occurred,  the  velocity  proflle  did  approach  a  uniform  distribution 
with  strong  changes  near  the  wall.  This  is  also  in  accord  to  general  trends  predicted  by 
the  one  dimensional  models  in  chapter  3  and  analyzed  in  chapter  4. 
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For  argon,  the  degree  of  ionization  is  very  low  («  \%)  for  most  regimes  studied,  due 
mainly  to  the  high  inlet  stagnation  pressure  of  1,000  atmospheres,  and  the  high  ionization 
energy  for  argon.  However,  at  2  eV  and  inlet  Mach  number  of  0  5,  argon  i.s  30.86% 
ionized,  and  this  energy  may  be  recovered  for  projectile  motion.  Considering  table  5.2, 
except  for  the  high  temperature,  low  inlet  velocity  case,  there  is  a  low  degree  of 
ionization,  making  recombination  energy  negligible.  The  argon  is  apparently  a  weakly 
ionized  plasma  at  1,000  atmospheres  at  temperatures  below  2  eV  and  at  high  velocities. 


Table  5.2  Argon  internal  energy  data  ( 1 ,000  atm) 


inlet  Mach  U 

iuletatag 

tfOnp  (K) 

inlet  lie 

(Jaqs) 

exit  sic 

(J/Kg) 

static  inkt 

temp(K) 

siAiicexit 

temp  (K) 

inlet  deg  of 

ionizjuton 

exit  deg  of 

»oniz«;k>a 

type  rad 

3.5 

10,000 

8.901x10’ 

1.138x10* 

1,964.73 

2,756.19 

«l% 

«i% 

Planck 

■E 

10,000 

s.goixKH 

1.373x10* 

1,961.25 

3,509.35 

«\% 

«i% 

Ross 

3.5 

23,000 

1.684x10* 

2.105x10* 

4,502.74 

6,12628 

«1% 

«i% 

Ross 

0.5 

10,000 

3.152x10* 

2  115x10* 

9,186.52 

6,136.72 

1% 

«i% 

Planck 

3.5 

23,000 

1.684x10* 

1.885x10* 

4,502.74 

5.146.82 

«\y 

«i% 

Planck 

0.5 

23,000 

L 

2.146x10’ 

6,013.71 

30.86% 

<:l% 

Planck 

It  is  apparent  that  the  MACH2  program  is  an  important  tool  in  the  study  of  ELGG 
barrel  regions  and  shows  results  following  the  general  trends  of  simpler  models.  The 
development  of  high  temperature  bands  near  the  walls  is  not  predicted  by  the  models 
developed  in  chapter  3,  which  is  a  consequence  of  the  more  sophisticated  techniques 
employed  in  the  programming  of  MACH2. 


CHAPTEP  VI 


CONCLUSION 


This  thesis  has  considered  the  thermal  characteristics  and  heat  transfer  within  an 
ELGG  barrel  section.  The  ELGG  examined  consisted  of  a  1.0  meter  barrel  with  .03 
meter  diameter  separated  from  an  arc  chamber  or  convergent-divergent  expansion  nozzle 
by  a  diaphragm  designed  to  burst  at  a  specific  pressure.  The  regimes  of  operation  studied 
were  inlet  stagnation  temperatures  from  8,000K  to  34,OOOIC,  inlet  stagnation  pressures 
from  850  atmospheres  to  1,150  atmospheres,  and  inlet  velocities  from  500  m/sec  to 
16,000  m/sec.  These  operating  parameters  are  far  in  excess  of  those  encountered  in 
typical  artillery  applications,  and  this  thesis  represems  the  first  in-depth  discussion 


desifmtffi  tn  <;n«»rifi<-allv  aHHr 


race  tWg  K-ot  transfer  mechanisms  within  an  ELGG  barr el 


section  for  these  conditions. 

The  heat  loss  processes  discussed  were  thermal  conduction  from  the  working  fluid  to 
the  surrounding  barrel  walls,  forced  convection  within  the  fluid  due  to  the  highly 
energetic  working  fluid  plasma  entering  the  barrel,  and  both  Planck  and  Rosseland 
radiation  losses  from  the  working  fluid  to  the  barrel  walls.  Viscous  effects  were  also 
considered,  as  was  the  chemical  recombination  effects  at  high  flow  temperatures. 
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pressures  and  velocities.  In  order  to  study  these  processes,  several  one  dimensional, 
steady  state  models  were  constructed  with  the  assumption  of  a  calorically  perfect  working 
fluid. 

Hydrogen  gas  was  selected  as  the  working  fluid  for  the  ELGG  due  to  the  high  natural 
sound  speed  and  low  molecular  weight  of  2.016x10'’  kg/kmol  Additionally,  hydrogen 
has  a  high  ionization  potential  winch  was  presumed  prior  to  the  beginning  of  tliis  thesis  to 
aid  in  the  energy  transferred  to  the  projectile  by  allowing  the  fluid  to  travel  farther  down 
the  barrel  before  losing  energy  due  to  transport  processes. 

This  assumption  was  proven  to  be  correct  when  inlet  stagnation  temperatures 
exceeded  20,000K  and  at  velocities  greater  than  8,500  m/sec  when  1,000  atmospheres 


was  used  as  the  stagnation  inlet  pressure. 


At  these  conditions,  it  was  dcleuniiicd  tliai  the 


energy  liberated  by  the  recombination  of  the  electrons  and  the  ionized  hydrogen  atoms 
could  offset  or  even  exceed,  in  some  circumstances,  the  losses  due  to  conduction, 
convection  and  radiation.  The  working  fluid  was  found  to  be  ionized  to  as  much  as  50 
percent  at  inlet  stagnation  coiiditions  of  30,000K  and  1,000  atmospheres. 

The  possibilities  of  non-equilibrium  of  the  working  fluid  were  detailed,  and  it  was 
found  that  due  to  the  high  densities  associated  with  the  inlet  conditions  mentioned 
previously,  as  well  as  the  fast  relaxation  rate:-;  the  working  fluid  was  able  to  reach 
equilibrium  conditions  throughout  the  barrel  region.  Frozen  flow  losses  were  shown  to 
be  negligible  as  well. 

The  determination  was  made  that  the  working  fluid  remains  absorption  dominated 
throughout  the  barrel  region,  suggesting  that  the  Rosseland  diffusion  model  accurately 
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predicts  the  radiation  losses  as  opposed  to  the  Planck  methodology.  This  finding  showed 
that  only  when  temperatures  drop  to  approximately  1  eV  does  Planck  radiation  begin  to 
be  applicable,  but  at  this  point,  the  radiation  losses  became  insignificant. 

Three  distinct  regions  of  operation  within  the  ELGG  baixel  section  were  determined 
where  various  transport  processes  dominated.  In  the  region  of  operation  with  inlet 
conditions  of  1,000  atmospheres  and  temperatures  up  to  1  eV,  conduction  dominated  the 
system,  with  radiation  limited  to  a  short  duration  due  to  the  rapidly  dropping  temperature 
of  the  flow.  Convection  tended  to  play  a  significant  role  in  this  region  only  at  highei  inlet 
velocities. 

The  second  region  of  operation  extended  from  1  eV  to  2  eV,  and  calculations 
demoustratcu  that  ail  three  transpon  processes  are  of  nearly  equal  magnitude  in  this 
domain.  Radiation  is  able  to  play  a  more  significant  role  due  to  the  dependence  on  T*, 
but  it  still  tended  to  drop  in  importance  as  the  flow  progressed,  whereas  conduction  and 
convection  remained  significant  throughout  the  barrel. 

In  the  third  region,  radiation  dominates,  with  conduction  contributing  a  significant 
amount  of  heat  loss  as  well.  Convection  remains  at  an  order  of  magnitude  less  that 
conduction  and  radiation  in  this  region.  It  was  determined  that  conduction  is  the 
consistent  heat  loss  mechanism  throughout  all  three  regions,  while  convection  tends  to 
have  the  largest  effects  at  higher  velocities. 

Viscous  effects  were  determined  to  be  present  at  all  phases  of  the  flow  through  the 
ban  el,  and  the  main  effect  was  to  cause  choking  of  the  flow.  However,  development  of  a 
theoretical  length  model  which  is  designed  to  identify  operation  regions  which  prevent 
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choking  of  the  flow  was  developed.  It  was  determined  that  the  conditions  for  choking 
the  flow  were  determined  by  the  friction  factor,  the  ratio  of  specific  heats  of  the  working 
fluid,  the  inlet  Mach  number,  and  the  length  of  the  barrel  itself. 

The  applicability  of  the  computer  code  MACH2  was  tested  for  use  with  ELGG  barrel 
sections,  and  data  obtained  from  MACH2  was  analyzed  with  trends  and  data  developed 
in  this  thesis.  Due  to  limitations  in  the  coding  of  M/.CH2,  hydrogen  data  was 
unavailable,  and  argon  was  used  as  the  working  fluid.  Despite  this  difficulty,  trends 
determined  in  the  hydrogen  models  and  analysis  provided  important  guidelines  in  the 
understanding  of  the  calculations  produced  by  MACH2. 

The  results  of  MACH2  computations  show  that  at  supersonic  flow  conditions  within 
the  ELGG  barrel,  a  band  of  high  temperature  forms  near  the  walls,  suggesting  that 
viscous  effects  produce  heating  of  the  fluid,  and  this  effect  is  seen  in  the  decrease  in 
enthalpy  loss  at  the  exit.  Subsonic  regions  showed  a  region  of  high  density  forming  at  the 
wall/exit  interface,  which  tended  to  actually  reduce  the  exit  area  of  the  flow,  and  in  effect, 
acted  as  a  sort  of  nozzle.  Velocity  data  produced  in  each  of  the  computational  ams  was 


able  to  verifv  t.his  effect.  Supersonic  flows,  hov/ever,  did  not  shov/  this  r 


density,  mainly  due  to  the  high  temperature  band  along  the  walls. 

All  flows,  however,  were  shown  to  develop  small  scale  directed  flow  from  the  axis  to 
the  walls,  and  a  fan  type  expansion  at  the  wall/exit  interface.  The  strength  of  this  directed 
flow  was  determined  to  be  directly  dependent  on  the  inlet  Mach  number  and  stagnation 
temperature.  Additionally,  flow  temperature,  pressure,  and  density  tended  towards  some 
steady  slate  condition  as  time  and  computational  cycles  progressed. 
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MACH2  proved  to  be  a  useful,  if  complicated,  tool  in  the  study  of  the  ELGG  thermal 
characteristics,  and  provided  insights  into  actual  system  operation  which  would  otherwise 
be  prohibitively  difficult,  if  not  impossible.  The  ability  to  study  transport  processes  either 
separately  or  simultaneously  is  an  import, ant  pari  of  the  code,  and  aids  greatly  to  its 
flexibility. 

Comparing  the  results  from  both  hydrogen  and  argon  chemical  reactions  (tables  4. 14 
and  5.2),  it  became  apparent  that  argon  makes  a  poor  choice  for  a  propellant  in  an 
ELGG.  Hydrogen  at  temperatures  from  2  to  3  eV  is  able  to  provide  significant  energy  to 
the  working  fluid  in  the  recombination  process  of  free  electrons  and  ionized  atomic 
hydrogen.  This  liberated  energy  is  thus  available  for  the  process  of  projectile  motion,  and 
allows  the  energy  to  be  "transported"  through  the  barrel  without  being  lost  as  conduction 
or  radiation  energy.  The  high  ionization  potential  of  argon  compared  to  hydrogen 
(15.81  eV  to  13.527  eV),  combined  with  the  high  pressures  in  the  ELGG,  prevents  this 
from  occurring  except  at  inlet  temperatures  in  excess  of  2  eV  and  low  Mach  numbers. 

The  results  of  this  thesis  show  that  the  ELGG  is  a  plausible  idea  which  can  take 
advantage  of  high  energy  plasmas  for  the  purpose  of  projectile  motion.  The 
determination  of  regimes  of  transport  dominance  is  essential  to  the  future  production  of 
ELGG  devices,  as  is  the  finding  that  chemical  recombination  can  be  an  effective  way  to 
minimize  transport  losses  at  high  temperatures  using  hydrogen  as  the  working  fluid. 
Finally,  by  carefully  tailoring  the  inlet  flow  conditions,  supersonic  flow  can  be  maintained 
throughout  the  barrel  region,  significantly  enhancing  the  possibilities  of  very  high  speed 
projectile  motion. 


APPENDIX  A 


Degree  of  Hydrogen  Ionization  Calculations 


These  calculations  wwere  perfonned  with  the  use  of  MATHCAD  4.0.  The 
methodology  utilizes  the  Saha  equation  with  hydrogen  temperatures  from  l.OOOK  to 
35,OOOK  and  pressures  of  100, 1,000,  and  10,000  atmospheres.  Also  included  in  this 
appendix  is  a  graph  showing  the  ionization  potential  lowering  as  a  fijnction  of 
temperature. 

(This  appendix  is  referenced  by  chapter  3.1,  page  23  and  chapter  4.1,  page  65) 
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APPENDIX  B 


SESAME  Plots 

This  information  is  obtained  from  the  Los  Alamos  New  Opacity  Library  and 
SESAME  Opacity  Data. 

(This  appendix  is  utilized  for  chapter  3.7,  page  52  and  chapter  4.5,  page  89 
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Figure  H.  I  SESAME  rosseland  opacity  for  hydrogen 
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Figure  B.2  SESAME  opacity  comparison  for  hydrogen 


APPENDIX  C 


Opacity  and  Blackbody  Radiation  Plots 


These  plots  were  created  with  MATHCAD  4.0  using  data  provided  by  N.  L.  Krascella 
and  R.  W.  Patch.  Other  Plots  courtesy  of  NASA-Langley. 


(This  appendix  is  referenced  by  chapter  3.7,  page  52,  and  chapter  4.5,  page  89) 
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Figure  C.3  Comparison  of  hydrogen  mean  free  paths  based  on  data 
by  Patch^  (Data  at  1,000  and  100  atmospheres) 
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APPENDIX  D 


Hydrogen  Composition  and  Data 

Hydrogen  species  composition  and  number  density  plots.  Plots  are  for  hydrogen 
number  densities  at  pressures  of  10  to  1,000  atmospheres  and  temperatures  from  300K  to 
40,000K.  Species  included  are  Hj,  Hj*,  H,  H*,  e',  H‘,  Hj*,  and  H*.  Data  obtained 
from  Patch^  and  INSPI. 

(This  appendix  is  referenced  by  chapter  2.2,  page  17,  chapter  3.1,  page  23, 
chapter  4.1,  page  65,  chapter  4.6,  page  93  and  chapter  6,  page  1 13) 
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Figure  D,  9  Hydrogen  composition  at  '  00  atmospheres 
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APPENDIX  E 


RosseUnd  DiiTusion  Approximation  Calculations 

These  calculations  are  performed  on  MATHCAD  4.0.  Data  is  for  hydrogen  gas  at  10, 
100,  500,  and  1,000  atmospheres  and  temperatures  from  1,000K  to  lO.OOOK.  Graphs  are 
comparisons  with  Rosseland  results  and  Blackbody  radiative  flux  predictions. 

(This  appendix  is  referenced  in  chapter  3.7,  page  52) 


Calcubitioa  oftbe  Radistivc  beat  transfer  using  Rosseland  Diffusion 
Approximation  at  1000  atmospheres.  Data  for  Rosseland  opacities  is  from 
Krascelia  and  Patch. 
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Calculate  the  radiative  heat  transfer  using  procccdui'c  in  Vmcienti  &  Kruger. 
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Figure  E.  1  Radiative  heat  transfer  to  wall  at  100  atmospheres.  QR  is  in  units 
of  Watts/meter^,  G  is  in  Watts/meter^,  and  D  is  in  millimeters. 
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Figure  E.2  Radiative  heat  transfer  to  wall  at  500  atmospheres,  QR  is  in  units 
of  Watts/meter,  G  is  in  Watts/meter^,  and  D  is  in  millimeters. 


APPENDIX  F 


CiUculation  results  for  conduction  model 


This  data  is  produced  with  Mathcad  4.0.  This  appendix  contains  one  sample  case  data 
and  computed  tabular  data. 


Variables  used: 

Molvvt  -  molecular  weight 

R  -  gas  constant 

y  -  ratio  of  specific  heats 

k  -  thermal  conducivity 

Mu  -  viscosity 

ITol  -  inlet  stagnation  temp 

IPo  1  -  inlet  stagnation  press 

VI  -  inlet  velocity 

Ipol  -  inlet  stagnation  density 

Aso  -  stagnation  sound  speed 

al  -  inlet  sound  speed 

Ml  -  inlet  Mach  number 

TI  -  inlet  static  temp 

PI  -  inlet  static  press 

time '  tranisit  time 

mass  -  mass  of  working  fluid 

pi  -  initial  density 

q  -  heat  transfer  per  unit  mass 

POlPl  -  ratio  of  press 

TOITI  -  ratio  of  temps 

Red  -  Reynolds  number 
hlO  -  inlet  enthalpy 

Le  -  Laminar  entrance  length 

(variables  with  a  2  in  them  indicate  exit  conditions) 

(this  appendix  is  referenced  by  chapter  3.5,  page  37,  chapter  3.6,  page  39,  chapter  4.2, 
page  75,  and  chapter  4.6,  page  93) 
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APPENDIX  G 


Calculation  results  for  radiation  models 

This  data  was  produced  with  Mathcad  4.0.  This  appendix  contains  sample  case  data 
and  computed  tabular  data  foi  Rosseland  diffusion  and  Planck  radiation  models. 


Variables  used: 

Molwt  -  molecular  weight 

R  -  gas  constant 

Y  -  ratio  of  specific  heats 

k  -  thermal  conducivity 

ivlu  -  viSCOsiiy 

iTo  i  -  inlet  stagnation  temp 

IPol  -  inlet  stagnation  press 

VI  -  inlet  velocity 

Ipol  -  inlet  stagnation  density 

Aso  -  stagnation  sound  speed 

al  '  inlet  sound  speed 

MI  -  inlet  Mach  number 

TI  inlet  static  temp 

PI  -  inlet  static  press 

time  -  tranisit  time 

mass  -  mass  of  working  fluid 

pi  -  initial  density 

q  -  heat  transfer  per  unit  mass 

POlPl  -  ratio  of  press 

TOITI  -  ratio  of  temps 

Red  -  Reynolds  number 

hlO  -  inlet  enthalpy 

Le  -  Laminar  entrance  length 

(variables  with  a  2  in  them  indicate  exit  conditions) 

(This  appendix  is  referenced  by  chapter  3.7,  page  52,  chapter  4.5,  page  89  and 
chapter  4.6,  page  93) 


147 


Hydrogen  (CHlorieally  perfec*.)  Planck  Rndi-'iloii  only  modci 
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Table  G. 3  Mathcad  rosseland  radiation  data 


APPENDIX  H 


Calculation  results  for  viscous  effects  model 


This  data  wcis  produced  with  Mathcad  4.0  This  appendix  contains  one  sample  case 
data  and  computed  tabular  data. 


Variables  used: 

Molwt  -  molecular  weight 
Y  -  ratio  of  specific  heats 
Mu  -  viscosity 
IPol  -  inlet  stagnation  press 
Ipol  -  inlet  stagnation  density 
al  -  inlet  sound  speed 
T1  -  inlet  static  temp 
time  -  tranisit  time 
pi  -  initial  density 
TOITI  -  ratio  of  temps 
Red  -  Reynolds  number 
f  -  fnction  factor 


R  -  gas  constant 

k  -  thermal  conducivity 

iTol  -  inlet  stagnation  temp 

VI  -  inlet  velocity 

Aso  -  stagnation  sound  speed 

Ml  -  inlet  Mach  number 

PI  -  inlet  static  press 

mass  -  mass  of  working  fluid 

POIPl  -  ratio  of  press 

hlO  -  inlet  enthalpy 

Le  -  Laminar  entrance  length 

LeT  -  theoretical  length 


(variables  with  a  2  in  them  indicate  exit  conditions) 


(This  section  is  referenced  by  chapter  3.4,  page  31,  chapter  4,3,  page  79,  and 
chapter  4.6,  page  93) 


155 


1  lydrogcit  (Ciilorically  |Krfcct)  Viscous  clfcvls  only  model 


5 


rknMubr>  Ls)cr  Thw-iiica  ii  I  mo«(  ufcrwb 
Tiri<  r*]ikU  ^  ODirwicu 


APPENDIX  I 


Calculation  results  for  convection  model 

This  data  was  produced  with  Mathcad  4.0.  This  appendix  contains  one  sample 
case  data  and  computed  tabular  data. 

Variables  used; 

R  -  gas  constant 
k  -  thermal  conducivity 
ITol  -  inlet  stagnation  temp 
VI  -  inlet  velocity 
Aso  -  stagnation  sound  speed 
Ml  -  inlet  Mach  number 
PI  -  inlet  static  press 
mass  -  mass  of  working  fluid 
q  -  heat  trarisfei  per  unit  mass 
TOITI  -  ratio  of  temps 
Le  -  Laminar  entrance  length 
-  wall  shear  stress 

(variables  with  a  2  in  them  indicate  exit  conditions) 

(This  appendix  is  referenced  by  chapter  3.6,  page  39,  chapter  4.4,  page  84,  and 
chapter  4.6,  page  93) 


Molwt  -  molecular  weight 
y  -  ratio  of  specific  heats 
Mu  -  viscosity 
IPol  -  inlet  stagnation  press 
Ipol  -  inlet  stagnation  density 
al  -  inlet  sound  speed 
T1  -  inlet  static  temp 
time  -  tranisit  time 
pi  -  initial  density 
POlPl  -  ratio  of  press 
Red  -  Reynolds  number 
hlO  -  inlet  enthalpy 
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1  !ydrog.:n  (Caloricaliy  pcrfcci)  Coiivcclion  only  model 
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Table  1. 1  Mathcad  convection  data 


Convection  Data 

FFCl|j  *  -  10X1^  =  POXl^  =  P0XI\  =  PRXl^  =  M.\I^ 


FFC3;.  ‘  \?0;^  =  T0X\  =  P0X3,^  =  P0X3Z^  -•  PRXi.  =  • 


JTC  =  VX4.^  =  TOX-!^  =  l><)X4^  -  POX->:^  ^  PlCX-i^  -  (}X-I^  = 


APPENDIX  J 


Calculation  results  for  barrel  theoretical 
length 

This  data  was  produced  with  Mathcad  4.0. 

(Ihis  appendix  is  referenced  by  chapter  3.4,  page  27,  chapter  3.6,  page  39,  chapter 
4.3,  page  79,  chapter  4.6,  page  93,  and  chapter  .S,  page  107) 
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Calculatioa  of  Th/joredcal  Length  o^ELGG  Barrel  based  on  friction 

Gas  data  for  Hydrogen; 

r  - 1.41 


Tube  data: 
d  =  .03  L  =  1 .0 


To  construct  the  graph,  I  will  calculate  the  friction  in  the  barrel  for  Mach 
numbers  ranging  f^om  M  =  .5  to  M=  8  in  steps  of  .05. 


MACH  .5 
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Calculate  the  theoretical  length  for  various  friction  factors  from  .001  to  .005 
in  steps  sizes  of  .0005: 
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Figure  J.3  Th?orcti.G5U  length  calculations  with  y  1.51  (hydrogen) 
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Figure  J.4  Theoretical  length  calculations  witn  y  I  6 .  (hydrogen) 
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Figure  J.  1  Theoretical  length  calculations  with  y  1.31  (hydrogen) 
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Figure  J  2  Theo^'etical  length  calculations  with  y  1 .4 1  (hydrogen) 
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Table  J.l  Theoretical  length  calculations  from  Niathcad 
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APPENDIX  K 


MACH2  Test  Results 

Note:  the  plots  shown  illustrate  the  ELGG  barrel  section  with  a  length  of  1 .0  meter, 
and  a  width  measured  from  the  axis  center  line  to  the  wall,  a  distance  of  .0 1 5  meter.  The 
proportions  of  the  plots  make  it  appear  that  the  radius  is  neariy  the  same  size  as  the 
length,  which  is  not  true.  These  plots  are  for  illustration  purposes  only,  and  are  not 
drawn  to  scale.  The  flow  enters  from  the  right  and  exits  to  the  left  of  the  plots  shown. 

(This  appendix  is  referenced  by  chapter  3.8,  page  63) 
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Figure  K.9  MACH2  Pressure  calculation  Cycle  5000  (L=1.0m  r=.015m) 
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Figure  K.IO  MACH2  Density  calculation  Cycle  33210  (L=1.0m  r=.0I5m) 
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Figure  K.l  1  MACH2  Temperature  calculation  Cycle  33210  (L=l  ,0m  r=,015m) 
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Figure  K.l 2  MACH2  Pressure  calculation  Cycle  33210  (L=l  Om  r=  015m) 
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Figure  K.  1 3  M  ACH2  sample  run  *lata 


APPENDIX  L 


MACH2  calculation  results  and  charts 

This  appendix  contains  sample  nins  for  two  cases  utilizing  th'^  full  effects,  with  Planck 
radiation.  The  im2  used  for  each  case  precedes  the  plots  which  were  produced.  The  first 
case  is  for  inlet  conditions  of  stagnation  temperature  10,OOOK,  stagnation  pressure  of 
1,000  atmospheres,  and  inlet  Mach  number  of  3.5.  The  second  case  is  for  inlet  stagnation 
temperature  of  23,OOOK,  stagnation  pressure  of  1,000  atmospheres,  and  inlet  Mach 
number  of  3.5. 

Note:  the  plots  shown  illustrate  the  ELGG  barrel  section  with  a  length  of  1 .0  meter, 
and  a  width  measured  from  the  axis  center  line  to  the  wall,  a  distance  of  .015  meter.  The 
proportions  of  the  plots  make  it  appear  that  the  radius  is  nearly  the  same  size  as  the 
length,  which  is  not  true.  These  plots  are  for  illustration  purposes  only,  and  are  not 
drawn  to  scale.  The  flow  enters  from  the  right  and  exits  to  the  left  of  the  plots  shown. 
Figure  L.  1  shows  the  calculation  mesh  and  dimensions,  and  is  representative  for  both 
cases 

(This  appendix  is  referenced  by  chapter  5,  page  107  ana  chapter  6,  page  113) 
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MACH2  ini2  for  full  model  using  planck  radiation 
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flow 


Figure  L.  I  Calculation  mesh  for  MACH2  calculation  (full  model) 


Figure  L.2  MACH2  planck  full  model  velocity  Cycle  1000  (T^=10,000K,  M=3.5) 


Figure  L.3  MACH2  planck  full  model  velocity  Cycle  5000  (T„=10,000K,  M=3.5) 
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Figure  L.4  MACH2  planck  full  model  velocity  Cycle  10895  (T^=10,00OK,  M=3.5) 


Figure  L.5  MACH2  planck  full  model  vorticity  Cycle  1000  (Tg=10,000K,  M=3.5) 
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Figure  L.6  MACH2  planck  full  model  vorticity  Cycle  5000  (T„=10,000K,  M=3.5) 
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Cycle  -  10895 
ABS  Vorticitv 
-  =  3.7E-K)2  'a='2.9E-K)5 
B=5.8E+05  C-17E-H)5 
D=1.2E-K)6  E=l.5E+06 
+  =  1.7E-K)6 


Figure  L.7  MACH2  planck  full  model  votrticity  Cycle  10895  (T^=10,000K,  M==3.5) 


Figure  L. 8  MACH2  planck  full  model  temperature  Cycle  1000  (T„=10,OOOK,  M=3.5) 


Figure  L. 9  MACH2  planck  full  model  temperature  Cycle  5000  (T„=10,000K,  M=3.5) 


Figure  L.IO  MACH2  planck  full  model  temperature  Cycle  10895  (T^=10,000K,  M-3.5) 


T  =  8.948E-04 
Cycle  =  1000 
Prc3S’jre 

-  =  2.5E+06  A=3.1E'+06 
B=4.0E+06  C=5.1E-K)6 
D=6.5E  -C6  E=8.3E-+06 
+  =  1.1E-K)7 


Figure  L.  1 1  MACH2  planck  full  model  Pressure  Cycle  1000  (To=10,000K,  M=3.5) 


T  *  4.495E-03 
Cycle  =*  5000 
Pressure 

--2.6E+06  A=2.9E-K>6 
B=3.3E-K)6  C=3.6E-K)6 
D=«4.1E-H}6  E=4.5E-K)6 
+  =5.0c+06 


Figure  L.  12  MACH2  planck  full  model  Pressure  Cycle  5000  (T,=10,000K,  M=3.5) 


T  =»  8.948E-04 
Cycle  =  1000 
Density 

-=  1.3E-K)0  A=1,7E-K)0 
B=2.3E-KX)  C=3.1E-KX) 
D==4.2E+00  E=5.8E-K)0 
+  =  7.8E-HX) 


Figure  L.  14  MACH2  planck  fiiU  model  density  Cycle  1000  (T,=10,000K,  M=3.5 


T  =  4.495E-03 
Cycle  =  5000 
Density 

-=  1.3E-K)0  A=1.7E+00 
B=»2.3E-+O0  C=3.1E+Oi 
D=4.2E-K)0  E=5.8E-HX) 
+  =7.8E+00 


Figure  L.  15  MACH2  planck  fiiU  model  density  Cycle  5000  (T„=10,000K,  M=3,5) 
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T  -  9.800E-03 
Cycle  =  10895 
Density 

-=1.3E-K)C  A»1.7E-KX) 
B=2.2E-K)0  C=2.9E-K)0 
D=3.9E+00  E=5.1E-H)0 
+  -6.8E+00 


Figure  L.  16  MACH2  planck  fiill  model  density  Cycle  10895  (T^=10,OOOK,  M=3.5) 


MACH2  iin2  for  ftili  model  using  planck  radiation 


Inlet  conditions: 

To  =  23,000K 
V  *  4,381  m/sec 
P  =  1.858  kg/m**3 
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Figure  L.  1 7  MACH2  planck  fiili  model  velocity  Cycle  1000  (T„=23,OOOK,  M=3.5) 


Figure  L.  18  MACH2  planck  full  model  velocity  Cycle  5000  (T„=23,000K,  M=3.5) 


Figure  L.19  MACH2  planck  full  model  velocity  C^J.  (T3=23,OOOK,  M=3.5) 
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Figure  L,2I  MACH2  planck  full  model  vorticity  Cycle  5000  (T„=23,OOOK,  M=3.5) 
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Cycle  =  10895 
ABS  Voracity 
-=1.8E+02  A=4,5E-K)5 
8=9.1E-K)5  C=l.4E-K)6 
D=-1.8E-H)6  E’2.3E-K)6 
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Figure  L.22  MACH2  planck  full  model  vorticity  Cycle  10895  (T<,=23,000K,  M=3.5) 
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Figure  L,23  MACH2  planck  ftiU  model  temperature  Cycle  1000  (T„=23,OOOK,  M=3.S) 


Figure  L.24  MACli2  planck  lull  model  temperature  Cycle  5000  (To=23,OOOK,  M=3.5) 


T  =  9.800E-03 
Cyc!-;  s*  1  Qg95 
Tetripcrature 
--1.7E4)!  A=4.9E-01 
B-5.7E-01  C=6.6E-01 
D-7.7E-01  E=-9.0E-0! 
+«1.1E+00 


Figure  L.25  MACH2  planck  full  model  temperature  Cycle  10895  (7^=23, OOOK, 

M=3.5) 


Figure  L.26  MACH2  planck  fiill  model  pressure  Cycle  1000  (T„=23,OOOK,  M=3.5) 


Figure  L.31  MACH2  pl-inck  full  model  density  Cycle  10895  (To=23,000K,  M=3. 5) 
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APPENDIX 


B.  Questionnaire 
BACKGROUND  INFORMATION 


NAME _ RANK _ 

NUMBER  OF  YEARS  MILITARY  FLYING  (INCLUDE  UPT) _ 

QUALIFIED  AS  A  CIVILIAN  PILOT  Y/N  IF  YES  UST  RATING _ 

UST  AIRCRAFT  FLOWN  IN  REVERSE  CHRONOLOGICAL  ORDER 
FORMAT:  ACFT  -  YEARS  FLOWN  -  HOURS 

EXAMPLE: 

F-16  89-92  600HRS  A-10  86-89  550HRS  T-38  84-85  150HRS  T-37  84  75HRS 
T-41  83  20HRS 


LIST  LOW  LEVEL  TIME  (EXCLUDING  UPT) 
FORMAT:  ACFT  YEARS  FLOWN  HOURS 

EXAMPLE: 

C-130  91-92  40HRS 


UST  AIRCRAFT  YOU  ARE  PRESENTLY  CURRENT  AND  QUALIFIED  IN 


HOURS  FLOWN  LAST  30  DAYS _  60  DAYS _  90  DAYS 


UST  ANY  AIRCRAFT  IN  WHICH  YOU  HAVE  OR  PRESENTLY  HOLD  INSTRUCTOR  STATUS 


TOTAL  FUGHT  TIME  CIVILIAN  &  MILITARY 


HAVE  YOU  EVER  FLOWN  THE  MS-I  SIMULATOR?  Y/N 
IF  YES,  WHEN  AND  FOR  WHAT  PURPOSE? 


94 


